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16. ABSTRACT 

This study is an investigation of the metastable auroral emission. The neighboring OH contaminants and 

low intensity levels of the emission itself necessitated the evolution of an instrument capable of separating the emission from 
the contaminants and having a high sensitivity in the wavelength region of interest. A new type of scanning photometer was 
developed during the study, and its properties are discussed in this report. 

The evolution of the scanning photometer involved theoretical and laboratory examination of transmission 
characteristics of all-dielectric, Fabry-Perot type filters. The photometer scanning principle employed a spatially moving field 
stop in the aperture plane of the instrument, creating an angular exposure of the filter to the detector. The moving field 
stop was accomplished by rotating a disk containing a spiral slit behind a spatially fixed, sector shaped, masking disk. The 
spatial scanning scheme accomplished a nearly constant peak transmittance and bandpass for large tilt angles. This advantage, 
coupled with the throughput advantage of the large Fabry-Perot type filter elements, offered a unique tool for examination 
of low light level emissions in contaminated backgrounds. 

An instrument was built and tested in the laboratory for verification of the theoretical concept. Another instrument, 
employing a multi-spiral configuration, was then prepared with its instrument design oriented to the development of a 
ground based system for the measurement of C^f^P-^D) auroral emissions. 

The theoretical aspects of auroral electron interaction with atomic oxygen and the resultant 0*^(^P-^D) emissions 
were examined in conjunction with NjlNEG emissions. The fact that Nj and O are differentially distributed in the 
thermosphere, and that 0'*'(^P) is quenched, lead to an investigation of the ratio of these two emissions as an emission 
signature relevant to the precipitating auroral electrons. 

Ground based measurements of 0"‘'(^P-^D) auroral emission intensities were made in College, Alaska, (64.9°N, 

148°N) using the spatial scanning photometer (sector spectrophotometer). Simultaneous measurements of nJiNEGi,o 
emission intensity were made in the same field of view using a tilting photometer. Time histories of the ratio of these two 
emissions made in the magnetic zenith during auroral breakup periods are given. Theories of I 7319 /I 4278 previous 
investigators were presented. These previous theories are not in agreement with the ground based measurements. 

A rocket measurement of NjlNEG©,© and O'^(^P-^D) emission in aurora was examined in detail and was found to 
agree with the ground based measurements. Theoretical examination resulted in the deduction of the electron impact 
efficiency generating P) and also suggests a large source of 0‘^(* P) at low altitude. A possible source is charge exchange 
of S) with 01 (^P). These new theories are different from those previously deduced by other investigators. 
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CHAPTER I 
INTRODUCTION 


The research described in this report involves two separate, but related, studies. 
The first study is that of a unique spatial scanning photometer system. The high 
throughput qualities found in flat plate interference systems and the near-constant 
transmittance and bandpass characteristics found traditionally in grating spectrometer 
systems are convolved properties of this unique system. The second study involves the 
investigation of the geophysical role of 0'*'(*P-*D) emissions in electron aurora. The 
photometer system described in the first study is one of the tools used to investigate 
oxygen ion metastable emission. 

The introductory chapter continues with a brief description of some spatial 
scanning schemes reported by other investigators, previous work done on the 0"^(*P-*D) 
in aurora (measurements and theories) and a brief description of the contents of the 
other chapters text. 

1.1 PREVIOUS INVESTIGATIONS IN SPATIAL SCANNING OF ALL-DIELECTRIC 

FILTERS 

The interference elements of this study are limited to solid spaced, all-dielectric 
filters. The spatial scanning discussion in the text, and history, will be confined to the 
field stop configuration in the aperture plane and- the shaping and positioning of that 
field stop in allowing imaged fringes, or portions of, to pass beyond the aperture plane. 

The pin-hole, or circular field stop, is the most commonly used field stop in 
Fabry-Perot type systems. For solid spaced systems, the circular field stop passes the 
annular shaped fringe for a normal orientation of the interference element. Father and 
Reasoner (1969) discuss extensively the photometer characteristics offered by tilting the 
interference element and translating the fringe pattern across the field stop. The tilting 
photometer offers the advantage of changing the peak bandpass wavelength of the 
photometer system and thus allowing a measurement of the background radiation 
underlying the emission of interest. In tilting the interference element, a bandpass 
distortion, or “broadening”, occurs which must be compensated for in the interpretation 
of data. The field of view and pointing direction of the system remains fixed in the 
tilting photometer scheme, 

A second type of spatial filtering has been accomplished by generating annular 
shaped field stops whose radial distances and annular breadths vary to accomplish a given 
resolution. Such a system is described in Shepherd et. al. (1965). Each annulus offers a 
different look angle through the interference element giving a bandpass shift at a constant 
resolution. Such a system requires a pointing lens to orient the look angle to the same 
sky position. A large number of annular stops are required such that appropriate peak 
transmittance stops can be achieved for a particular interference element. 



1.2 EARLY MEASUREMENTS OF Oil AURORAL EMISSION 


Identification of the 0'^(^P-*D) emission in aurora was first reported by Omholt 
(1957). Omholt reported on several spectra obtained with a grating spectrograph from 
observations made at Yerkes Observatory in 1956. On two of the spectra a broad feature 
at X7325A was noted and believed to be the O'*" doublet. He noted that the doublet 
appeared most prominent in the ‘upper’ part of the aurora. He speculated that since the 
feature was not notable in the lower parts or in all spectra that the emission was 
probably quenched. 

Dufay (1959) reported several measurements of a feature at X7327A, which he 
identified as the O"^ doublet. These measurements were made in 1957. For a particular 
aurora of September 30, an estimate of some relative intensities was given. Much like 
Ohmolt’s measurements, these were primarily of an identifying nature and not of 
sufficient resolution to completely resolve the doublet. 

A measurement was reported by Wallace (1960) from Yerkes Observatory for a 
December, 1958 aurora. This measurement was used to substantiate the presence of the 
emission to support the identification of the 0'*'(^P-^S), 3727A doublet in a spectrograph 
he was analyzing of the February 10/11, 1957 great ‘Red Aurora’. In the December event 
reported (although no related emission spectra were reported) the 0"'‘(*P-*D) doublet was 
resolved and relative line strengths of the doublet were compared and verified with the 
Seaton and Osterbrock (1957) theoretical values. 

Vallance Jones (1960) did obtain spectrograms of the ‘Red Aurora’ of February 
10/11, 1957 at Saskatoon. The Jones’ spectra was measured with a high dispersion 
instrument with a spectral slit of ~ 5.4A. Jones identified relative intensities of many 
other atomic and molecular emissions including theNj IPGs ,3 at 7384.2A. Intensities of 
interest to this study were the 0‘*'(^P-^D) wavelength positions at 7319.4A, and 7330A, 
and NjlPGsjS at 7384.2A with relative intensities of 17.0, 11 . 8 , and 2.7, respectively. 
According to Jones (1971), no further auroral measurements of the doublet have been 
reported. This particular aurora has been extensively studied. lesooAss?? ratios were 
reported between 3 and 10 by Shuisakaya (1967) in the 60-70® geomagnetic latitude 
range. Jones (1971) summarized various auroral spectral classes, and this type of aurora is 
estimated to have a typical lower border height of 250 km. The Jones’ (1960) data will 
be examined more extensively in Chapter VIII. 

Recently, Walker (1975) reported a measurement made on an Atmospheric 
Explorer Satellite from which pertinent 0'*'(*P) continuity parameters were deduced. 
These results will be discussed in Chapters VII and VIII. 

1.3 EARLY THEORIES OF O+I^P-^D) AURORAL EMISSION 

The atomic transition parameters for O’*" are given by Seaton and Osterbrock 
(1957). Electron impact ionization cross sections for O'*" from ionization threshold to 
several hundred electron volts were theoretically determined by Seaton (1959). The total 


ionization cross section was verified by Fite and Brackmann (1959). According to Seaton, 
the probabilities of ionizing to a ^P, ^D, or **S state are .2, A, and .4, respectively. More 
recently, Kazaks et. al. (1972) calculated a breakdown of 3:5:4, respectively, for the ^P, 
*D, and '*S states. These values for ionization probability are for E > 300 ev, or the 
energy regions where the Born approximation describes the ionization cross section. For 
perspective, examination of the total ionization cross sections reviewed by Kieffer and 
Dunn (1966) suggests the total ionization cross section ratio of O to Nj is ~ .5 for E > 
300 ev. 


Chamberlain (1961) has postulated that quenching of 0‘*'(^P) must be large in 
light of the auroral measurements discussed above. Dalgamo and McElroy (1965), in 
treating dayglow analysis, estimated daytime intensities with and without a large 
deactivation rate ( 1 O'* ° cm^ sec'* ). They found resultant intensities consistent with their 
earlier estimates which required a large deactivation rate. Walker and Rees (1968) studied 
ionospheric electron densities and temperatures in aurora. In their analysis, the N 2 
deactivation rate of 0"'’(^P) was assumed to be 10'*° cm® sec'* . The high quenching rate 
is usually referenced to Dalgamo and McElroy (1965) (i.e., Jones, 1971; and Gerard, 
1970). Dalgamo and McElroy, however, did not give an argument for the high 
deactivation rate in their dayglow study. Since they first used the high rate, however, 
reference will be made to their work for the rate. 

Gerard (1970) made a study of the metastable oxygen ion and related optical 
emissions based on the Seaton (1959) cross sections and the Dalgamo and Mchl’oy 
(1965) quenching rate of 10'*° cm® sec'* Gerard related the volume emission rates of 
the oxygen metastable 0‘*'(®P-®D) and 0"*'(®D-'*S) to NJINEGq,© for a theoretical 
particle flux distribution of an IBC II type aurora. Gerard made comparisons of his 
computed results to spectra observed by Omholt (1957) and Wallace (1960) to verify his 
computed volume emission rates. Gerard’s interpretation of the spectra used for 
comparison verified his analysis. 

1.4 A BRIEF SUMMARY OF THIS WORK 

A study of 0'*’(®P-®D) emission was undertaken by this author to further clarify 
0+(2p) metastable chemistry in aurora. The study involved the investigation of both 
ground based and rocket measurements of 0"'’(®P-®D) emissions in aurora. 

The ground based measurements were made with a unique photometer system 
which has application to investigations of other atmospheric glows. The first half of this 
report (Chapters II through V) describes this system. 

Chapter II is aimed at giving the reader a fundamental understanding of the 
theoretical aspects of the spatial scanning photometer. Transmittance and bandpass 
characteristics are derived and described for the instrument in terms of transmission 
half-widths rather than analytic descriptions of the entire transmission function. The 
individual broadening functions and the convoluted instrument function are discussed. 
The achievement of a nearly constant instrument function for all slit posifions is probably 
the most significant advantage of this scheme over other spatial scanning schemes. 
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The all-dielectric filter properties were investigated in the laboratory, independent 
of the photometer, to get an understanding of their properties as they apply to the 
instrument. The particular emphasis in the filter study was on the angular characteristics. 
These findings are the subject of Chapter III. Characteristics were found to be widely 
variable from filter to filter. Theoretical aspects of these elements are discussed in 
Appendix A. 

A single spiral photometer was fabricated and studied in the laboratory to test the 
theoretical hypothesis of Chapter II. This study is reported in Chapter IV. Several 
objective configurations and filters were employed in the system and all are discussed in 
this chapter. A Fortran IV program is listed in Appendix B which generates the spatial 
slit functions described herein. 

Finally, a ground based spectrophotometer was built with the instrument 
characteristics oriented to the measurement of 0"'‘(*P-^D) in aurora. A multi-spiral 
scheme was employed in the ground based system. The main advantage over the single 
spiral device is an increased throughput for a given resolution. This ground based system 
is described briefly in Chapter V. A thorough understanding of Chapters II and IV is 
necessary for the reader to fully understand what is presented in this chapter. 

The last half of this report describes a study of the auroral 0"^(^P-^D) emission. 
The emission is weak and is surrounded by emissions of the Nj IPG band system. Ground 
based measurements are further plagued with OH emissions emanating at wavelengths 
close to the emission wavelengths of 0'*'(^P-^D). 

The theory of 0‘*’(*P-^D) is presented in Chapter IV. The continuity equation for 
the 0‘*'(^P) state is derived, with flux divergence considered negligible. The probability of 
C^C^P) was related to the probability of NjlNEG band emissions for convenience in 
analyzing the measurements of the subsequent chapters. Appendices C and D provide 
detailed 0"*'(*P) transition probabilities and production by electron impact information 
respectively. 

In Chapter VIII, ground based observations of 0’'‘(^P-^D) and N^lNEGi^o 
emissions are presented and analyzed. Theories of 0‘'’(^P) for aurora are presented in 
terms of ratios of predicted I 7 3 1 9 and I4 2 7 s • Appendix E is dedicated to the details of 
reducing the ground based measurements. Appendix F describes the volume emission 
assumptions used for the NjlNEGi^o emissions as based on previous investigations. 

Chapter VIII reports the detailed analysis of a rocket measurement of the 
0 +( 2 p- 2 £)) and NJINEGojO emissions. The efficiency of producing 0'*'(^P) by electron 
impact and quenching extremes were deduced from the data. The resultant values for 
these parameters are radically different from previous investigator’s findings. Because of 
this, the data reduction and error analysis are belabored in Appendix G to lend 
credibility to the findings. The N 2 IPG band system used in this analysis is listed in 
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Appendix H in the Fortran IV language. Several ‘New Models’ are discussed in the 
Chapter which better describe these measurements than the models of previous 
investigators. 

Chapter IX summarizes the primary conclusions arrived at in this study. The 
unique characteristics of the sector spectrophotometer system as well as auroral 
O+C^P-^D) investigations are discussed. 

At the beginning of this study, two goals were visioned. The first was to explore a 
spatial scanning scheme for Fabry-Perot type interference generators. The desired 
characteristics of the scanning scheme included exploring the extent of wavelength range 
available in a single interference filter. Also desirable in the spatial scanning were a 
constant (wavelength independent) throughput, constant bandpass, and near continuous 
wavelength selectability or tuning of the instrument. In essence, an instrument having the 
scanning flexibility of a grating spectrometer but having the light gathering power of a 
Fabry-Perot aperture was desired. The second goal was to apply that instrument to 
measuring low light level emission, in aurora, that was of geophysical interest but difficult 
to examine with conventional techniques. The goal was to measure from the ground the 
©■•■(^P-^D) emission intensity at X7319A to discern the radiative loss rate of the doublet. 
The continuity of production and loss of O'^C^P) in aurora was expected to be 
enlightened by such a measurement. A rocket measurement of the emission in aurora was 
added later in the study. The combined measurements afforded a unique opportunity to 
postulate some new concepts of 0"'‘(^P) continuity in aurora. 
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CHAPTER il 

SECTOR SPECTROPHOTOMETER INSTRUMENT THEORY 


2.1 INTRODUCTION 

The intent is to design a spectrophotometer with a high throughput and spectral 
resolution for faint light sources in the near infrared and visible wavelength region. The 
system ideally should have a constant resolution and throughput over the scanning range 
of the instrument. Ideally, the system would be a continuously scanning device which 
would permit spectral analysis and selectable exposure interrupts for time integrations at 
wavelengths of interest. Such a photometer is favorable for studying atomic and 
molec ular emissions of chemical reactions originating from the earth’s upper atmosphere. 

A photometer was developed by Bather and Reasoner (1969) which employs a 
tilting principle to scan in wavelength. In their system, the aperture plane lies at the focal 
distance of the objective lens, which is fixed behind an interference filter. By tilting the 
filter, the fringe pattern passes over the aperture and the detector is exposed to a 
changing wavelength. In their photometer. Bather and Reasoner employed a dielectric 
interference filter to generate a fringe pattern. Filter characteristics were studied in the 
tilting filter photometer and filter properties suited for particular optical characteristics 
were investigated. 

The interference fringe pattern generated by plane reflective surfaces of the 
Fabry-Perot type is characterized by a spatially narrowing fringe with increased radial 
distance from the objective axis in the focal plane (Bom and Wolf, 1970). The dielectric 
interference filter generates an interference pattern but contains a solid spacer versus the 
gaseous type normally used with Fabry-Perot plates. The filters are also equipped with 
blocking elements to suppress the transmission of any wavelengths outside a specified 
band, that is, all orders are suppressed except one. Given a monochromatic source, a 
single fringe will be generated with a dielectric filter and the characteristics of this fringe 
can be described very similarly to the gas spaced Fabry-Perot. The spatial breadth of the 
fringe will depend on the radial distance from the objective axis intersect in the focal 
plane of the objective. A monochromatic source of a shorter wavelength would generate 
its fringe at a larger radial distance than a longer wavelength, due to the fixed spacer 
index and spacing. 

By tilting the filter, the fringe pattern moves in the focal plane and a fringe 
generated by some monochromatic source can be made to pass over a slit in the focal 
plane of the objective. Bather and Reasoner (1969) employed this principle in their 
tilting-filter photometer. In their instmment a spatially fixed, circular aperture was used. 
Since the aperture was fixed, as the fringe pattern passed over the field stop, a greater 
breadth for short wavelength fringes was allowed to pass through the aperture, resulting 
in an instrument broadening, varying with tilt angle, and consequently, wavelength. To 
design a mechanical scanning instmment with a constant resolution, one must accoimt for 
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the spatially narrowing fringes generated at radially larger distances (i.e., at shorter 
wavelengths) from the objective axis intersect in the aperture plane. These effects have 
been compensated for in the sector spectrophotometer explored in this report. 

2.2 INSTRUMENT THEORY (THE. IDEAL INSTRUMENT) 

2.2.1 APERTURE THEORY, WAVELENGTH AND THROUGHPUT 

The relationship between the angle of tilt of a filter and the normal axis to a 
filter is reviewed in Appendix A. Equation A- 10 of Appendix A can be put in 
abbreviated terms, i.e., 


AX _ 1 

Xo n*^ 2 


( 1 ) 


where 


= - — — (P and Q are defined in Appendix A) 

IX* — referred to as the effective index 
d — tilt angle (in radians) . 


The slope of the wavelength shift versus 0* (0 in degrees) will be defined as a, where 


_ XqTT^ 
2 m *2 ( 180)2 


It therefore follows that: 


AX = «02 (0 in degrees) 


(3) 


or 


AX = a02 (0 in radians) 


(4) 
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where 


X 
2m** 


( 5 ) 


A schematic of a filter, objective configuration employed by most photometers, is 
shown in Figure 1. This particular configuration will image an infinitive source at the 
aperture plane where the field of view is described by the aperture area and the objective 
aperture plane separation. The following definitions apply: 


At — area of interference filter 
A 2 — area of the aperture (field stop) 
f — focal length of the objective lens 
ri — radius of interference filter 

r — radial distance of the aperture slit from the objective axis intersect with 
the aperture plane 

6 — previously defined . 


The throughput of such a configuration is given by Born and Wolf (1970) as: 


L = Aw 


( 6 ) 


where L is throughput, A is area of cone of intercept, and w is the solid angle of cone of 
intercept. For a displaced aperture, as indicated in Figure 1 , the throughput becomes 


L = cos'* d . (7) 


In equation (7), note that for a given input aperture size (Ai) and focal length (f), A 2 or 
the aperture area must increase for L to remain constant. The most commonly used fixed 
wavelength photometers employ a circular pinhole for A 2 and fix 0 = 0“, 
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SOURCE 



Figure 1 . Optical schematic and symbol definitions of a 
source imaged on an aperture plane through an 
interference. 



Examination of an arc-aperture or annulus sector, as depicted in Figure 2, will aid 
in evolving the principle of a unique system having considerable advantages over fixed 
field stops. In Figure 2 is shown a sector of a circular shaped slit with Sr breadth at r 
distance from the objective axis intersect with the aperture plane subtending an angle 
500- The area of such a sector is given by: 


Aj = r8rS0(, 


( 8 ) 


If one could move such a slit from r = r^jj^ to some r maximum, but decrease Sr as the 

slit moves to account for the changing breadth of projected fringes, a constant throughput 
with scan could be approached. That is, since 


r = f tan 0 


(9) 


where, f - objective lens focal length. Equation (7) can be rewritten as 


L = 



( 10 ) 


where 


cos d 


p;- and sin d 

(r* +f*)/* 


r 

(r* + f®)‘/* 


( 11 ) 


A spatial scheme to move the slit in the 50® sector from rjjjjj, to may be 

envisioned as follows. One can approximate the arc aperture with a “spiral” on a disk 
spiralling from r = to Visualize a second disk with a “pie” sector cut from it, 

overlaying the spiral disk. Now, by rotating the spiral disk overlayed by the spatially fixed 
disk, a spiral chord which moves from r = rj^in* *■ “ ^max ^ ®xposed in the sector “pie”. 

In order to scan in wavelength in uniform time, the relationship 


r = f tan 0 



( 12 ) 


1. # 0 since Sr must be finite demanding in turn that r be finite. 


10 




Figure 2. An aperture schematic for a radially shaped slit of 6r 
breadth as viewed normal to the aperture plane through a 
sector with angle 8<I>q . 
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must be examined. Consider turning the disk at a uniform rate from 0 = 0° to 
requiring 


(p — j3(X " Xq) 


( 13 ) 


or 


^ M = constant 

5t 5t 


(14) 


where 


_ ^max 
^^max 


constant 


(15) 


It follows then that 



But, 


(16) 


Vax - 2 7rn 

where n = number of turns in the spiral from r = rj^jjjj to 
rotations per scan.^ 

Thus, for a constant wavelength scan, it can be written that 


(17) 

or the number of 



(18) 


2. It is possible to have n greater than one if some mechanism is employed to mask the 
unwanted spiral chords in the sector exposing only a single chord to the detector for a given 
spiral position. 
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Another way of expressing Equation (18) may be derived by manipulating Equation (4) 
to read 


d 



(19) 


or 


max 




max 

a 


( 20 ) 


resulting in 


r = f tan 


(^max ) 


( 21 ) 


Expressions (18) and (21) both describe the criterion for r with respect to (p for a 
constant wavelength scan rate in a fixed rotation rate spiral system. Note in Equation 
(21) that for (p = 0„,ax> 


*^max ^ ®max 


( 22 ) 


2.2.2 APERTURE THEORY, BROADENING CHARACTERISTICS 

Since a 8<Pq arc is exposed in the sector, there exists an integral number of 5p 
increments in that arc sector, i.e., 

S4> = |36X [A derivative of Equation (13)] . (23) 

For the S<po of the sector then, 

5(^0 ~ PS\q (24) 
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where 5X^ is the broadening introduced by approximating the circular fringe pattern by a 

q>iral chord. In a sense, the instrument broadening dependence on tilt mentioned earlier 
has been minimized, but a curvature broadening as described above has been introduced. 
Equation (24) can be rewritten to read: 


^ -is ^^max 
5^0 0 2 Trn 


or 


s . ( 26)3 

AXj„ax 2 7m 


This is a square function broadening with the curvature flnesse now written as 



2 Trn 
5^ 


(27) 


The broadening introduced by a finite aperture (breadth broadening-SXg) may be 
computed by first approximating Equation (1) as 


X - Xo = AX = a sah* 0 


(28) 


The breadth broadening, written in terms of 0 follows from Equation (28) as 


5Xg = 2a sin d cos 0 S0 . (29) 

In order to express Equation (29) in terms of r. Equation (9) is differentiated to read: 

S6 » 5°^ fir (30) 

3. AXjjj{^(/6X(^ is analagous to finesse, when speaking of a Fabry-Perot system (Bom and 
Wolf, 1970). In this instrument, finesse (N) will be defined as AX^m^/fiX or in this case, the 
curvature finesse is = AXp ^ay /^^C* 

14 



which may be substituted into Equation (29), along with the definitions of sin 0 and cos 
0, to yield: 


6Xb 


2 arf^ 5 t 
(r^ +f2)2 


(31) 


Using a finesse definition analogous to the curvature broadening, the breadth finesse can 
be written as; 


5Xg 2 arf* 6r 


(32) 


2.2.3 TOTAL INSTRUMENT FUNCTION 

The ideal instrument is defined as one which satisfies a constant aperture 
broadening function, a constant instrumental throughput, and a constant wavelength 
change for a fixed spiral rotation rate, i.e., 

8\^ = (5X^+6X|)‘^* = constant , (33) 

L = ^ r5r60o ^ = <=°*^stant , (34) 


or 


rSrS^o 

(r2 + {2)2 


constant 


AtJ e- ^ 
f*Ai f2 


(34a) 


where Uq = 03 at 0 - 0° , and 


r 


f tan e^ax 


2 irn. 


where 8<j>l&t - constant. 


(35) 
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From the previous relationships of and 5Xg, Equation (33) can be rewritten 
as 

" ■ 06) 

However, from Equation (34a) 



r5r 


coo(r^ + P)^ 
f*60o 


(37) 


permitting a rewrite of Equation (36), i.e.. 




(38) 


It is shown, therefore, that with a constant throughput requirement, the breadth 
broadening is inversely proportional to Thus, for constant (i.e., a “pie” shaped 
section, 6^o) the aperture function is constant [Equation (38) J. 

The aperture is then constructed by selecting an instrument finesse. 


XT ^^max 


(39) 

The aperture finesse becomes 



1 1 ^^max 

NX Nf 6X^1 

^^max 

(40)^ 


4. 6X|| and 5X^ are filter characteristics contributing to the instrument function and are 
discussed in Appendix A. 
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By setting the two aperture finesses equal, the aperture finesse becomes: 


Nb - Nc - v^TRa 


( 41 ) 


The sector angle can then be computed from Equation (27) depending on n selected and 
the spiral may be computed from Equations (31) and (35). 

Another way of expressing Equation (37) is given by: 


6r = Sr, 


max 


( r^ + f^ \ *‘max 
fmax + / r 


(42) 


where 


Sr = 




2 aT™„„ f* 
* max 


(43) 


Expression (43) serves as a convenient form for computing the spiral breadth as well as 
does Equation (31). The radial function r is computed as a function of <j> [Equation 
(35)] to satisfy the constant wavelength scan rate criteria. 

Broadening for the total instrument (instrument function) is due to the 
convolution of each broadening component in the instrument. The instrument function 
for a filter-aperture is made up by the convolution of these functions, i.e.,® 


SXi = (SX|j + 6Xp‘/* 


(44) 


where 


SXj — instrument function, 
SXp — filter broadening, and 
SX^ — aperture broadening . 


5. An approximation. The aperture function is a square function broadening. The 
approximation involved in this convolution impUes the filter broadening is a well behaved 
function, i.e., square, triangular, or Gaussian. 
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A relationship that is constant for a constant throughput has been derived. 

Note, however, that for Acu to be constant with 6, the cos d dependence of A must be 
compensated by the inverse in cu. cu is related to the cone-half an^e (^) by 

' 03 = . (45) 

Note from Equation (A-13) that 

or GKw) , 

or, rewriting Equation (44), using Equation (A-13), it follows that 

1/2 

8Xi = (8Xf| + 6X5,+ 6X| + 6Xy 


(46) 


(47)« 


The instrument apertures will then allow 


6 X 





(48) 


where - cone 1/2 angle at 0 t 0°. In terms of r and f, Equation (48) becomes 


6X4 


V ‘/a 

P Vo(r"+f^) 
1-Q/ f 


(49) 


For SXj to be truly constant, then 6 X^ must be kept small relative to 8 X|j and 5X^. 

Therefore, once the aperture finesses and parameters have been computed, one should 
compute the solid angle of the instrument (^) and check 5X^ (49) for cone angle effect 

on the total instrument function. Note that f controls r for a given A 3 and may be 
altered without affecting the aperture finesse. 


6 . Note that broadening due to polarizing has been ignored. 
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2.3 INSTRUMENT THEORY (MODIFIED) 


The apertures, buUt for the spectrophotometers used in this study, employed a 
slight modification to the theory. 

It was stated after Equation (44) that for a truly constant throughput, the cos d 
change in A (A| cos 6) results in an inverse cos 6 effect in (i.e., co/cos 0). If co changes, 
the cone half angle or field of view of the instrument is allowed to change. This is not a 
desirable trait for an instrument to be used in the aurora. 

Equation (8) was modified such that 


L = Aw = A,w , (50) 

which results in 

w = ^ cos* 9 = constant , (51) 


or, in terms of r and f. 


CO = = Wo 

(r» + f*)' 


(52) 


By substituting Wo [Equation (52)] into Equation (38), Equation (36) would become 



and Equation (31) would become 


6Xb 


2 arf r 
(r* +f>)*/* 


(54) 
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In other words, by allowing L to change with cos 6, a constant w and field of 
view are achieved. The sacrifice is for a 6Xg varying as cos d [Note Equation (54) differs 

from (31) by a cos 8 factor] . A constant field of view also removes the cos 6 effect from 
the cone half angle broadening function [Equation (48)] . 
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CHAPTER Ml 

FILTER CHARACTERISTICS - LABORATORY ANALYSIS 


The purpose of this study was to analyze, in the laboratory, a number of available 
filters and examine their characteristics for potential instrument use. The filter 
broadening and transmission characteristics as a function of tilt angle were of primary 
interest. The transmission peaks of filters used were limited to the 5000-7500A 
wavelength range. Both cryolite and zinc sulfide spacers are used in this region. Below ~ 
4500A, only cryolite is used due to dispersion in zinc sulfide (Lissberger, 1959). 

The cone angle, temperature, dispersion, and spacer index effects have been 
studied by numerous authors (as referenced in Appendix A). The effective index and its 
use has also been explored by Lissberger and Wilcox (1959), Pidgeon and Smith (1964), 
Bather and Reasoner (1969), and others. 

A monochromator arrangement, as depicted in Figure 3, was used to examine 
these filter characteristics. A collimated beam from the monochromator was passed 
through the filter and then a capture lens (not shown) imaged the beam on a detector. 
The filter was removed and the source intensity was determined and stored in a 
computer. The filter was then replaced and the source was again scanned for transmission 
versus wavelength computations. 

The effective 0^ slope versus A\ is shown in Figure 4 for a high index (\q = 
6296A) filter. Note that up to 0 = 30°, the slope is relatively constant as has been 
verified by others. The transmission bandpass (5X|() is shown in Figure 5. It was noted 

that at 36°, the transmission slope showed two peaks with a very distorted shape as if 
the beam had polarized. Note a gradual increase in bandpass with tilt angle. This 
particular filter was not examined as to whether the gradual broadening was due to 
polarizing or some other effect. 

Another filter (low index) was examined and the more rapid wavelength shift 
versus 0^ is shown in Figure 6. Broadening with tilt is more extreme than the high index 
filter discussed previously, as displayed in Figure 7. 

It was noted in the filter section that the low-index filters should theoretically 
polarize at lower angles for a given order filter. The order of these filters is unknown but 
the index trend is apparent. What is not explainable is that the broadening is apparent at 
very low angles of tilt, and follows a near 0^ effect rather than the very abrupt polarizing 
effect derived in the filter theory (Appendix A) with respect to the polarizing 
broadening. 

This effect is furthermore born out in the \o = 6308A (low index) filter shown in 
Figure 8. In this case, the transmission profile was examined at 0° and 10° tilt, with and 
without a polarizer. Even at 0° tilt to the beam, this transmission curve shows an obvious 
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F^re 3. Schematic of the laboratory arrangement used for filter 

transmission studies. 




WAVELENGTH SHIFT (A) 



N) 


Figure 4. Wavelength versus 0* as deduced from the laboratory 
arrangement shown in Figure 3 for a Xq = 6300A filter. 
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Figure 5. Halfwidth (full width at half maximum transmission) versus 
wavelength measured at incremental tilt angles for the 
Xo ® 6300A filter shown in Figure 4. 
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Figure 6. Wavelength shift versus 0* for a high index filter 
(Xo = 5896A). 
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Figure 8. Transmission shape of a low index filter (Xo = 6308A) with extreme 
polarizing characteristics. It is noted that even at 0 ° tilt, a distortion 
in the transmission profile exists which is removed with a 

polarizer. 
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distortion of the transmission profile, due to a polarized splitting. At 10®, the polarizing 
effect is very notable. One also notes that the polarizer did not account for all of the 
broadening. It can only be speculated that this non-polarized broadening, as observed 
with a collimated beam, might be due to roughness of film deposition, resulting in a 
changin g reflectivity of a given film with angle. This is only speculation but users of the 
very narrow band filters should be aware that the effect exists. 

In summary then, there are two effects not well defined. The polarized effect 
(SXp) does not adhere to the simplified theory examined in Appendix A. The trend that 

the high index filters (Zj^S spacer) polarize at much higher tilt angles than the low index 

(cryolite) is verified, althougli as noted in the last case, the degree and angular 
dependence are not well defined by the theory presented in Appendix A. The 
non-polaiized broadening as a function of tilt is not explained. The investigation shows 
approximately 2 to 3 Angstroms gradual broadening prior to a tilt, corresponding to 
transmission distortion of a non-Gaussian nature. The rate of this change seems to follow 
a 0* relationship. 
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CHAPTER IV 

A SINGLE SPIRAL SECTOR SPECTROPHOTOMETER 


The previous two chapters have outlined the theoretical aspects of a spatial 
scanning photometer, and theoretical and laboratory investigations of some all-dielectric 
filters have been discussed. The next step, then, was to build a fairly simple instrument, 
primarily for laboratory investigation of the spatial scanning technique. This chapter 
describes the concept, fabrication, and testing of the first sector spectrophotometer. 

4.1 BASIC CONCEPTS AND DETECTOR 

It was decided to specify the initial spectrophotometer as having a single spiral 
(n=l) and individual finesses of 20 (i.e., Ng=N(;;=20). Equation (27), of Chapter II, now 
sets 500 at 18°.’ 

The radial limit of r (rjjjax^ taxes the requirement for the optics in the 
capture-detection portion of the instrument. The increased with increased 

photomultiplier detection surface, with nearness of the detector to the aperture, and with 
decreasing angle (0) of acceptance through the aperture, which depends on f and r. In 
this case, an EMI 9558 photomultiplier tube was used, offering a circular detection 
diameter of 44 mm. The tube was cooled to -20°C to reduce the dark count requiring a 
minimum physical separation distance of 9 cm between the photomultiplier tube face and 
the capture lens because of the cooler constraints. The determining factor remaining in 
ascertaining was the capture lens, the lens diameter, and the angular dispersion of 

photons through the aperture to be focused on the phototube detector. A further 
restriction on of course, was that it must not exceed the small angle approximation 

in required in the theory section. 

4.2 CAPTURE AND OBJECTIVE OPTICS 

The capture lens imaged the objective lens on the photomultiplier (since the 
source is diffuse). Figure 9 depicts a sketch of the lenses and detector where Sj is fixed 
at 9 cm as previously mentioned. 

It was decided to restrict the objective lenses to reasonably thin lenses and to 
shelf items. Objective lenses of 166 mm, 106 mm, 88 mm, and 57 mm were selected, and 
all objective combinations were used with the same aperture, capture lens, and detector 
combination. With Sj (Figure 9) fixed, a capture lens focal length of near 45 mm was 
approximated as being an optimum focal length for the objective selected by inspection 


7. An aperture configuration is specified by the number of degrees a spiral rotates in going 
from 0 = 0 to 0 = ^max> where 0niax ~ ^Trn. For n = 1 , 0max ^®™ 

the 2ir or single spiral spectrophotometer. An n = 3 system will be discussed later which will 
be classified as a multi-spiral spectrophotometer. 
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Figure 9. Schematic of the capture optics used in the sector system. 



of the lens maker’s formula (Jenkins and White, 1957). This is only a crude 
approximation, since the thin lens approximation is not totally justified. On a bench, a 
number of capture lenses were tried with the objectives chosen. A non-focused image on 
the detector resulted in a wander of the image on the tube face as the slit was moved in 
the aperture plane. A lens with a focal length of 52 mm was chosen for the capture lens 
and an of 25 mm (~ 1 inch) was selected as a result of the bench study. The 

capture lens was a plano-convex lens with the flat side positioned up against the aperture 
to minimize angular loss with various positions of r of the slit. The limit was such 

that all positions of r from to would permit capture on a mock tube face at 

the detector position in the bench study. Noticeable wander of the image across the tube 
face did occur for the 57 mm objective, but the remaining objectives yielded no 
detectable wander of the image. The ^niax’s s®l®cted Flax’s from 8° for 

the 166 mm objective to 23“ for the 57 mm objective. The optics chosen are summarized 
in Figure 10. 

Note in Figure 10 that the capture optics are offset from the objective axis. The 
scan involves moving the aperture radially away from the objective axis, which necessarily 
centers the capture lens over the sector. This is best put into perspective by examining 
Figure 1 1 , which shows the capture lens and optics as viewed from the detector. Since 
the object (the objective lens) is at an angle to the capture lens axis, the detector must 
be offset from the capture lens to capture the image. 

An additional point on the “semi-focused” image on the detector should be made. 
The intent is to count the total number of photons passing through the aperture. They 
can be counted as long as they hit the detector surface of the tube, which ideally would 
be equally efficient over all areas of its surface. Most photomultipliers decrease in 
efficiency to the outer extremes of the tube face. For a wandering image, efficiency 
variations in the photomultiplier must necessarily be accounted for. 

4.3 INTERFERENCE FILTER 

The filter used in the spectrophotometer was a Xq = 6300A dielectric, which had 
drifted towards the blue to 6296A (peak transmission at 0 = 0“). It was felt that the 
filter could still function for nightglow spectroscopy in looking at the (9,3) hydroxyl 
band at X6356A. 

The filter was first studied, independent of the spectrophotometer system, to 
determine its characteristics. The 6300A filter described in the previous chapter was used 
for this purpose. This filter was a high index (Spacer - Zj^S) all-dielectric with oc = .26. 

The above filter and its role in the photometer system are discussed in this 
chapter. Before proceeding, examination of the effect of the index of refraction of the 
spacer in the system may be enlightening to other potential users of the technique. 
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Figure 10. Schematic of the single spiral sector spectrophotometer optical 
arrangement. The objective focal lengths used in the 
spectrometer are indicated at the top of the figure. 








Should a low index (cryolite) spacer be employed with the same 8X|| as above, an 

effective index of 1 .34 would yield an a of .5 (a factor of two larger slope). For a given 
objective, the scan range would have been doubled. However, as noted in Equation (31), 
SXfi = G(a) requiring a larger f for a given 6Xg. For faint source work, a slower scan rate 

would be necessary to match the count rates detected for the spectral range of the high 
index filters, since now the throughput would be less. In addition, the angle at which a 
filter polarizes, as noted earlier, occurs at lower tilt angles for the low index filters. 

4.4 APERTURE AND SECTOR 

The sector, or fringe mask, is best fixed very close to the spiral disk to 
accommodate the planar concept derived in the theory. The aperture must be rotated 
behind the sector in the spectrometer arrangement (Figure 12). 

The sector was made by overlapping two semicircles. The two semicircles 
comprising the sector disk and the aperture disk were generated by computer plots, 
etched on rubylithe masters, photographically reduced, and chemically milled on 3 
mil-304-rolled stainless. The sector and spiral disks were then mounted into a block so 
that the sector was fixed and the sprial could be rotated with a d.c. motor. 

It should be noted in Figure 12 that the aperture spiral can be in a position 
which exposes to the detector, the two ends of the spiral. The portion of a rotation for 
which this exists is 18° or 1/20 of a rotation and will henceforth be termed as the 
transition period. This transition takes place when returning to the high wavelength 
portion of the scanning cycle to begin a new scan. 

The final assembly is depicted in Figure 13. The optical portions were enclosed in 
a box and the detector cooler was butted into a press fit o-ring in the rear plate of the 
instrument, 

4.5 CALIBRATION TESTS 

4.5.1 LABORATORY CALIBRATIONS 

The test set-up used in generating the intensity and transmission profiles is 
depicted in Figure 14. Note that in these tests, one is observing what the instrument, as a 
whole, is doing to the shape of near-delta function (~ lA) input provided by the 
monochromator. In simpler terms, the transmission profiles are a convolution of the 
integrated effect of the interference filter, aperture, and capture optics on the source. 

Figure 15 depicts the transmission curves with the 166 mm, 106 mm, and 88 mm 
objective configurations.® The scatter in the transmissions is largely due to the inability 


8. The 57 mm configuration required a spacer in the instrument that prohibited the 
removal of the filter in the setup without disturbing the source. True transmission profiles 
could not be calculated as a result of this physical restriction. 
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Figure 1 2. Sector and aperture wheel configuration for an n = 1 spiral. 
The sector disk is fixed and the spiral rotated. 
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Figure 13. Sketch of the single spiral spectrophotometer. The 
aperture-sector wheels were mounted in a block to which 
the motor was attached. The capture lens was mounted 
in another block. All components were then fastened 
to a metal plate and boxed in a light tight 
environment. 
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Figure 14. Schematic of the laboratory arrangement used for the study of the 
single spiral spectrophotometer. 
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Figure 15. Transmission profiles versus wavelength for three objective 
configurations used in the single spiral system. The 57 mm 
objective configuration is not shown since true transmission 
profiles were not feasible due to physical 
limitations of the instrument. 
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to replace the source precisely, since it had to be moved between the filter-out and 
filter-in scans from which given transmission curves were generated. 


The trend of the peak transmissions is to decrease with decreasing wavelength or 
higher 0. For example, in the extreme case, the 88 mm objective configuration showed a 
transmission change from .45 at the high wavelength end to .33 at the low wavelength 
end. This loss is partially attributed to the same loss observed in studying the filter alone. 
Whether the photons are entering the filter specularly or diffusely, the filter itself 
generates a forward lobe in the transmission. The reflective loss at the filter is greater for 
higher angles (or lower wavelengths). 

The halfwidth of the transmission curves for the configurations are summarized in 
Figure 16. Note that some broadening occurred with increased angle for the 88 mm and 
106 mm configurations, but a narrowing with wavelength exists with the 57 mm 
objective. The explanation, as to the broadening slopes for the various configurations, is 
attributed to the fact that there is a combined filter and aperture effect. As pointed out 
earlier, there exists a slight angular dependence of filter broadening possibly due to 
irregularities of the interference surfaces in the filter itself. For the longer focal length 
objectives, the instrument function is primarily due to the filter. Figure 5 of the previous 
chapter displays the angular broadening of this particular filter. In the 57 mm 
configuration, a dominating aperture effect is observed. In the ideal aperture [Equation 
(36)] , the breadth broadening expands as (r^ + f^ )^ /r and in the modified slit used here 
[Equation (53)], the breadth expands as (r^ + f^) ^ /r. The aperture then is not a 
perfectly unchanging resolver of the spectral line. The modified sht is narrower at a given 
r than the ideal sht to maintain a constant field of view discussed in Chapter II. The 
modified sht then should correspond to a narrowing bandpass as is observed in the 57 
mm configuration. For a given objective configuration, that is for a given filter and 
objective distance, the spiral aperture could be modified by opening or closing the sht 
produced by the sector to yield a constant broadening by the instrument over the 
spectral range. 

Degradation of transmission over the spectral range of each configuration was 
questionable with the monochromator arrangement, due to the specular nature of the 
source. A diffuse source was placed in front of each arrangement and scans of the 
continuum were made. (The source was a tungsten lamp behind multiple layers of opal 
glass). The transmission loss over the spectral range, as summarized in Table 1, was 
computed from this arrangement rather than from the monochromator. 

A summary of the instrument under the various configurations is shown in Table 
1. The counts/sec/rayleigh are for the peak transmission at the high wavelength position 
of each objective configuration. A Carbon 14 source, which had recently been calibrated 
by NBS, was used for this calibration. 
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Figure 16. Summaries of halfwidth versus wavelength of the objective configurations 
with a 6300A filter. The 57 mm halfwidth is deduced from intensity 
profiles as was done for the other three configurations. 


TABLE 1. SUMMARY OF SPECTROMETER CHARACTERISTICS WITH 6300A 

FILTER (a =.26) 


Objective 
Focal Length 

Spectral 

Range 

AX,/, 

®max 

Transmission 
Loss Over SR 
(approx) 

Aw (<an* 
Steradian) 

^^/Ray 

166 mm 

17.5A 

4.6A 

00 

6% 

.0073 

— 

106 mm 

43A 

6.3A 

13® 

24% 

.018 

16.5 

88 mm 

62A 

1.5k 

15.6° 

30% 

.025 

32.5 

57 mm 

140A 

20.0A 

23° 

50% 

.056 

















A second filter was examined for the purpose of a nightglow observation to be 
discussed next. The filter was a X6305A filter, whose characteristics are summarized in 
Table 2. The reason this filter was used was to capture the 01, X6300A line emission, as 
well as, the (9,3) hydroxyl band for comparison. The a was very near that of the Xo = 
6300A filter, A calibration plot of wavelength versus wheel position (in degrees) is 
plotted in Figure 17. The beginning of a scan is generated by a microswitch trigger on 
the scanning aperture wheel. Note the 1 8° of transition discussed earlier and the starting 
wavelength of 6300A. The starting wavelength is less than the normal filter transmission 
maximum, since the minimum 0 is and not zero degrees. 

4.5.2 A NIGHTGLOW OBSERVATION 

Figure 18 depicts a section of nightglow data taken from Broadfoot and Kendall 
(1968). Note in the range scanned are the 01, X6300A emission and the Pi , 2 , Pa,a» Qj. 
and Qi lines of the OH (9,3) band. 

The data in Figure 19 were taken with the single spiral sector spectrophotometer 
on August 11, 1971, at Ann Arbor, Michigan. The intensity of the 6300A line was 
measured to be 50 rayleighs and 15 rayleighs for the P 2 ,a line of the OH (9,3) band 
system. The scan is coming up on the Q 2 , Qi lines just prior to transition, but the peak 
is not reached prior to transition. This observation was made with the 106 mm 
configuration as summarized in Table 2. This particular scan was made with a seven 
minute scan rate at 1 1 ;00 p.m. E.S.T. on the date mentioned. The sky was very clear and 
the moon had not yet risen when the observation was made. It is noted that the Pjj 2 
line at X6287.6A was not evident on this night in the observation but is very prominent 
in the Broadfoot and Kendall data shown in Figure 28. It is not the purpose of this 
report to examine the OH in detail. Had a filter comprising a low index spacer of the 
same resolution been employed, the remainder of the band system, including the R1-R2 
branch, could have been scanned. 
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TABLE 2. SUMMARY OF SPECTROMETER CHARACTERISTICS WITH 6305A 

FILTER (a =.24) 
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Figure 18. Nightglow data observed by Broadfoot and Kendall (1968) with the 
range of the single spiral spectrophotometer superimposed. 
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CHAPTER V 

A MULTI-SPIRAL SECTOR SPECTROPHOTOMETER FOR MEASUREMENT 
OF AURORAL O+C^P-^D) EMISSIONS 


5.1 INTRODUCTION 

A ground based instrument was then designed and built to accomplish a particular 
measurement [0‘^(^P-*D), 7319A] in the auroral zone. The theory of the emission will 
be discussed in the next chapter. 

The object was to build a spectrometer with a high throughput and sensitivity 
with two inch optics to examine the 0'''(*P-^D) emission and to scan and, thus, separate 
the nearby OH and Nj IPG background emissions from the signal. 

5.2 OPTICS 

A high index filter peaked at Xq = 7325 A with a 5X|| = 5.1 A was acquired for the 
instrument. 

The detector was an EMI 9658R tube, with a .30 mm diameter sensitive area. The 
tube was housed in a thermo-electric cooler to minimize dark count. Glass lenses were 
used in the system. The capture lens was selected after a bench study similar to that 
discussed for the single spiral instrument. 

5.3 APERTURES AND SECTOR 

For this machine, a 3n (6 tt) spiral was chosen to expand the throughput (i.e.. 
Aw = 3x larger for a given curvature finesse). A second spiral disk (shutter spiral), 
rotating at 1/3 the rate of the aperture spiral, was necessary to mask off two of the three 
spiral sectors exposed to the detector at a given aperture spiral position. In summary, the 
aperture seining configuration involves a fixed sector disk, a rotating n=3 spiral aperture 
disk, and a rotating n=l spiral mask disk. Figure 20 depicts the aperture schematic. The 
value for was set at 25 mm for detector capturabUity similar to the single spiral 

system. The masking spiral was generated by examining the boundary of the aperture at 
the sector extremes for a given angle and enlarging by a small increment to insure that 
no vignetting of the aperture mask would take place. 

The two rotating wheels were individually rotated by stepping motors. A single 
pulse translator was used with a divide by three circuit going to the mask spiral. 
Microswitches were used on the apertures as a scan fiducial and for a mask alignment 
pulse. The wheel alignment was electronically checked and realigned after each scan by 
synchronizing pulses generated by the microswitches. 
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5.4 GENERAL DESCRIPTION 


An instrument schematic is shown in Figure 21. The primary optical components 
are indicted. 

A generalized electronics flow chart is shown in Figure 22. The dashed lines 
indicate the preferred data acquisition scheme, however, a magnetic tape recorder was not 
used in gathering the data reported in this report. 

Six thousand motor steps resulted in one scan, making the scanner essentially 
continuous. The main advantage of stepping was position accuracy, digitizing 
convenience, and versatility such as computer control of a stepping sequence or exposure 
interrupts. 

The instrument was mounted on a Cassegrain telescope yoke mount and a 
synchronous motor control of zenith angle was used to point the instrument. 

The scanning scheme involves accepting light at various angles through the filter. 
A large fresnel lens, with its focal point being the filter, was attached to the front of the 
instrument in order to keep the instrument from scanning across the sky in the process 
of looking at various angles through the filter. 

5.5 CALIBRATIONS ^ 

A resolution summary is shown in Table 3. A range summary for other objectives 
is shown in Table 4. The transmission profiles as measured at various angles are shown in 
Figure 23. Laboratory calibrations were made by an arrangement similar to that described 
for the single spiral spectrophotometer described in Chapter IV. 
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1. - Fil 

2 . - Objective Lens 

3. - Rotating Spirals 

4. - Capture Lens 

5. - Detector 

6. - Spiral Stepping Motors 

Figure 21. An assembly schematic of the optical, scanning, and detecting 
components used in the multi-spiral sector spectrophotometer 
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Figure 22. A general electronic flow diagram outlining the main components of the 
field system. The dashed lines indicate the capability in data 
acquisition, although a tape recorder was not available 
for use for data reported in this study. 

































TABLE 3. SUMMARY OF RESOLUTION CHARACTERISTICS OF THE 
MULTI-SPIRAL SPECTROPHOTOMETER (X® = 7326A) 


Filter 

u* « 1.9 
Xo = 7326A 
6X||= 5.1A 
Instrument 

f (objective focal length) =127 mm 


®min 2 1/2 




Wavelength 

(A) 

Instrument Temperature 
(°C) 


22°C 

8°C 

Xo (Filter) 

7328A 

7326.1A 

^in 

7326.5A 

7324.7A 

^max 

7301 .3A 

1299 Ak 


6Xj (Theoretical) 
5Xj (Measured) • 










TABLE 4. SUMMARY OF RANGE CHARACTERISTICS OF THE 
MULTI-SPIRAL SPECTROPHOTOMETER 
(Xo ~ 7326A, u* = 1.9) 



^min 

^max 



Scan Range 

(A) 

(approximate) 

Objective 
Focal Length 
(mm) 

''mm 
(deg sq) 

''mm 
(deg sq) 

^^min 

(A) 

^^max 

(A) 

57 

5°40' 

22°0' 

Q Q 

149.4 

140 

32 

484 

y .y 

88 

3°36' 

16°p' 

4.0 

66.5 

62 

13 

256 

106 

2° 5 8' 

13° 20' 

2.7 

46.3 

43 

9 

177 

127 

2°30' 

11°12' 

1.9 

32.6 

30 

6 1/4 

125 

166 

1“55' 

8° 35' 

1 1 

1 Q 4 

18 

3 2/3 

73 1/2 

1 .1 

1 


NOTE: If Xo ~ 3660A, the last 3 columns are multiplied by 1/2. If u* ~ 1 .3 
(low index filter), the last three columns are multiplied by 2. 
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Figure 23. Transmission versus wavelength as measured for incremental aperture 
positions. The 0 profile- refers to the upper wavelength position of the spirals 
(i.e., 9 - Positions 1,2,3, and 4 refer to the number of motor 

steps from 0 (in thousands) the spiral was stepped for the indicated 
transmission curve. Beyond step 4000, the spiral was in transition 
and shuttered at 5000 steps and began a new scan at 
6000 (or 0) steps. 




CHAPTER VI 

THEORY OF 0+(*P-*D) EMISSION IN ELECTRON AURORA 


6.1 INTRODUCTION 

The study of aiuroral emissions from ground stations has been a fundamental and 
useful tool in analyzing many physical phenomena of the aurora. With increasing 
knowledge of ionization cross sections by electron impact, constituent distributions, and 
measurement capabilities, further refinements of the physical understanding of the aurora 
can be recognized from ground observations. 

The distribution and morphological characteristics of auroral primary electrons 
have been reviewed and summarized by Rees (1969). Rees summarized the measured 
electron energy and pitch angle distributions. In general, the energy distributions are 
representable by a power law function with a wide range of function coefficients 
representing the various measurements. The power law distribution takes the form 


G(Eo)«E?= 


(55) 


An exponential form is also used where the power law distribution fails and takes the 
form 


G(Eo) ocEoe-Eo/“ 


(56) 


Rees (1963) theoretically analyzed electron streams and deduced ionization rate height 
profiles for various primary electron energy distributions. Berger et. al. (1970) performed 
a detailed study of the energy disposition of auroral electrons in the atmosphere. More 
recently. Banks et. al. (1974) have performed a detailed analysis of auroral ionization 
calculated from a measured incident electron distribution which included the secondary 
distribution and its energy degradation with height. 

Height distributions of ionization rates can be deduced by examining emission 
rates arising from the fraction of ions generated in excited states. N^INEG emissions 
have been extensively examined for this purpose due to their prominence and simplified 
vibrational and rotational structure (Jones, 1971). Borst and Zipf (1970) have shown that 
the ratio NjlNEGo.o electron excitation cross section to the total ionization cross 
section is .07 for electron energies ranging from ~ 30ev to lOKev. Since collisional 
frequencies of the atmospheric constituents at auroral altitudes are much less than the 
emission lifetime of nJ1NEGo,o>® the emission rate is indicative of the electron 

9. The collisional frequency at 90 km is ~10^ sec"* . The transition probability is 1. • 10’ 
(Shemansky and Broadfoot, 1971) which corresponds to a lifetime of 9.1 -lO"® sec. 
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ionization rate for precipitating auroral electrons. The energy required for generating the 
from the N 2 <X‘ Sg) is ~ 19ev. Banks et. al. (1974) suggest that 60 to 70% of 

the auroral ions generated are due to the primaries in the upper altitudes of precipitation, 
and the remaining 30% are due to E < SOOev electrons. The percentage of primaries 
responsible for ionization decreases at lower altitudes and the percentage due to 
secondaries increases. 

Historically, the OI(‘D-^P) i.e., X6300A , emission has been used, along with 

NJINEG emission, as a photometric indication of auroral electron hardness or 
atmospheric range (Eather, 1969, Rees and Luckey, 1974). The oxygen emission is a high 
altitude phenomena relative to the nitrogen emission due to the differential constituent 
distribution of the atmosphere. Secondly, the OI(‘D-^P) emission is quenched at lower 
altitudes leading to a height dependence of the two emissions. 

The weakness of the above technique is two fold. First of all, the OI(‘D) has a 
110 sec lifetime, is not bound to the magnetic field, and diffuses or plumes in discrete 
forms. The allowed N^INEG emission does not diffuse prior to emission. The flux 
divergence effect of OI(‘ D) then can be significant. Interpretation of magnetic zenith or 
distant arc photometry measurements should consider divergence when using OI(‘ D) as a 
hardness or range indicator of precipitating electrons. A second weakness lies in the 
possible source functions for the 0(‘D) state. The 1.96ev threshold is significantly 
different from the ~ 19ev threshold required for N^INEG emissions generated from 
ground state molecules. As the downward flux of electrons changes in altitude or the 
initial distribution varies from one aurora to another, preferential generation of one 
emission over the other is possible. 

The ©■'■(^P-^D) emission is studied here to better understand its geophysical 
characteristics in aurora and to explore its use as an electron hardness indicator in lieu of 
the OI(‘D-*P) emission discussed above. The 0'*’(^P-*D) has several advantages over the 
OI(^D-^P). The lifetime, or reciprical emission probability, is ~ 5 secs versus 110 secs for 
OI(‘D). Consequently, it will not diffuse in discrete forms to the degree OI(‘D) will. 
Secondly, it is an ion and the displacement will primarily be controlled by magnetic field. 
Another advantage is that the 0"''(*P) has an excitation threshold ~ 18.5 ev from ground 
state 01. The excitation energy, consequently, is very near that required of NjlNEG 
emissions. There should not be the problem, therefore, of the emission being generated 
from a different portion of the electron energy distribution as is the case for OI(*D). 
Furthermore, the 0'*‘(^P) arises from neutral oxygen and is quenched as is OI(‘D-^P) 
emissions. 

Among the disadvantages in the 0‘*'(*P-*D) emission is its relatively low intensity 
compared to OI(‘D-^P) and NjlNEG emissions (Chamberlain, 1961). Further 
complication arises from the fact that it underlies a strong N 2 IPG band. When observing 
the emission from the ground, nearby mesospheric OH contaminants complicate the 
ability to resolve the 0"^(^P-^D) emission. 
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The theory of 0"*‘(*P-*D) emissions as related to NjlNEG emissions is derived in 
the remaining portions of this chapter. The derivation of the volume emission rate 
relationship of these two emissions includes an investigation of what historical 
investigations have found for the quenching rates, probabilities of emission, and pertinent 
excitation cross sections. New measurements and new found values for some of these 
rates are discussed in the two subsequent chapters. 

6.2 CONTINUITY OF 0+(*P) IN ELECTRON AURORA 

Since diffusion rates are small relative to the spontaneous lifetime below 300 km, 
the continuity equation will be represented as a simple production and loss. 

6.2.1 LOSS - RADIATIVE 

The spontaneous cascade of 0"''(*P) to lower O’*" states is given by 


0+(*P) ^ 0+(*D,'‘S) + ^ 


(57) 


where 


— Spontaneous emission probability, 
h — Planck’s constant, 
c — Speed of light in vacuum, 

X — Emitting wavelength . 


The 0'*‘(*P) has two substates (0'*'(*Pi/2) and O'^C^Ps/a)) whose statistical weights are 2 
and 4, respectively. The 0'*‘(^D) also has two substates (0‘*’(^D3/a) and 0'*’(^Ds /a)), 
whereas, O'^(^S) has one (0'*'(^S3/a)). The subscripts on the term values above refer to 
the J value, or total angular momentum quantum number associated with the particular 
substate. As a result, there are six possible loss channels from the ©"^(^P), three from the 
*Pi/a state, i.e., 

*P,/a-^*D3/a , *Ds/a , or '‘S^/a (58) 
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and three from ^Pa/a state, i.e., 


^Ps/j"^ * 1 ^ 3/2 , *Ds /2 , or ‘’S 3/2 


( 59 ) 


All six probabilities are finite for both electronic quadrapole and 'magnetic dipole 
transitions, Seaton and Osterbrock (1957) have provided the atomic parameters for 
0 +( 2 p) in Table 1-C of Appendix C. A diagram describing the spontaneous loss from 
©■'■(^P) is shown in Figure 1-C of Appendix C. The total radiative loss of 0'*‘(*P) can be 
written as: 


7?(0+*P)r = t7(7319A) + t7(7330A) + t?( 2470A) = ,219nIO+2pj (gQ)*®’" 


In this study, particular volume emission losses were experimentally measured in 
aurora and these fractional losses of 0'*’(^P) are related to the total 0"''(^P) population. 
The particular losses were determined in Appendix C. Specifically, the total 0'^(^P-*D) 
loss can be written as: 


t?( 7319-30A) = .171 n[0+2p] 


(61) 


(See Equation C-11), or 


T?(7319A) = .098 n [0+2 P] 


(62) 


(See Equation C-13). 


10 . r)i )j refers to the volume rate of change of the constituent state inside the 

parenthesis or the volume emission rate in the case of a wavelength. The subscript i refers 
to radiation loss rate (i=R), to quenching loss rate (i=Q), to a production rate only (i=P), 
or to a loss rate only (i=L). Absence of a subscript implies that all volume changes are 
considered. 

11 . The wavelengths 7319A, 7330A, and 2470A reflect closely spaced wavelength pairs 
since these pairs are not resolved spectroscopically in this study, i.e., 

T?(7319A) = r?(7318.6A) + t?(7319.4A), 
t?(7330A) = t?(7329.9A) + t 7(7330.7A), and 
r?(2470A) = t?(2470.3A) + t?(2470.4A) . 

Another combination employed in this study combines the two pairs combining all 2p-*D 
emissions, i.e., t?(7319-30A) = n(7319A) + t?(7330A). 
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6.2.2 LOSS RATES - QUENCHING 


The following reactions were considered as quenching losses for the ion 
metastable: 

0+(*P) + Nj T?*0+(*D,^S) + Nj (63)‘» 

( 7 jq^ 10"‘® , Dalgamo and McElroy, 1965.) 


and 


0+(*P) + e ^ 0+(*D,^S) + e (64)**’*^ 

( 7 g = 3 • lO"*, Seaton and Osterbrock, 1957.) 

12. Walker et. al. (1975) recently deduced a preliminary value of 7 jq^ 5 • 10"'® from 

analysis of 7319A measurements taken on Atmospheric Explorer C Satellite. Their data 
analysis was for dayglow passes of the satellite. 

13. Seaton and Osterbrock (1957) determined the deactivation coefficient for inelastic 
electron collisions in the form of the collisional strength parameter ^(n,n') (Hebb and 

Menzel, 1940, and Seaton, 1953). The deactivation rate takes the following form: 

Qnn' “ ’'^^(nn')Z^n *^n 
where 

Kjj = 2 Trm v„/h, 

m = electron mass, 

Vj^ = initial electron velocity, 

ojj^ = statistical weight of level n, 

h = Planck’s constant, 

J2(n n') ~ collisional strength parameter for states n, n'. 

The ratio ranges from .06 to ~ .15 for substate collision possibilities. A value 

of .148 for the ratio corresponds to the dominating ^Pi/a - transition. The Qn,n' 

expression for 0'*’(*P) reduces to: Qnn' = 1.3 • 10**^^ when applying the collisional 
strength factor, etc., provided by Seaton and Osterbrock (1957). For a range of Tg from 
300“K to ~ 3000“K, Qj^' takes the values of 2.2 to 5.5 • lO"* respectively. 

14. Henry et. al. (1969) deduce a value for 7 g of ~ 2 • 10"’ for Tg ~ 300“K. 
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Two loss mechanisms for O**" include 


O’*' + Nj NO’*’ + N , (65) 

(7N2 = 2 • 10-‘2 , Fite, 1969) 


and 


O’*’ + Oj ^2* Oj + O (66) 

( 70 j = 2 • 10-“ , Fite, 1969) 


These rates have been determined in the laboratory or deduced in the atmosphere where 
0’*’('*S) is the dominant ionic state. It has not been determined, to the author’s 
knowledge, whether these rates (65 and 66) are accelerated or decelerated for the 0’*’(^P) 
state. Without a reason to speculate them to be different, it will be assumed they cannot 
be less for 0’*’(^P). 

Recently, Walker et. al. (1975) also suggested a possible alternative to large Nj 
quenching (See footnote (12)), i.e.. 


0’*’(2P) + 0(3 P) T2 0’*’('*S,3D) + 0(3 P) . (67) 

( 70 == 2 • 10'*® , Walker et.al., 1975) 


It was suggested by Walker et. al. (1975) that if 7 q was as large as 2-10"‘ ®, then 7 j^^ 

5-10"“ rather than 5-10-‘® as indicated in footnote 11. At any rate, the day glow 
analysis of Walker et. al. suggest that large quenching of 0’*’(3p) by either Nj or O was 
required to account for the measured surface brightness profiles of 0’*’(*P-3D) emission 
measured on Atmospheric Explorer Satellite. 

The quenching loss for 0’*’(3p) can be written as 

t?(0-*-3P)q = S 7in[il n[0’*’3pj (68) 
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where 


7 j — collisional deactivation rate of P) by constituent i, 
n{i] — volume density of constituent i . 


The quenching possibilities were described above in expressions (63) through (68). 

The total continuity loss, therefore, can be written by adding the radiative loss 
(60) to the collisional loss (68), i.e., 

7?(0+*P)l = .219n[0+2p] + Z 7in[i] n[0+2p] (69) 

i 


6.2.3 PRODUCTION OF 0+(*P) BY ELECTRON IMPACT 

The next consideration will be the production of 0‘*‘(^P) by electron impact, i.e.. 


OI + e 0'*’(®P) + e + e 


(70) 


The relationship for producing 0‘*'(^P) by electron impact is derived as a function 
of the volume emission rate of N^INEGq,© Appendix D, i.e.. 


7j(0+"P)p = R„ (E) • t7(3914A) • 

^ ^0 niNol 


where 


(71) 


R„,(E) 




^(E) 


“Nr® 


(72) 


The cross section factor (Rj^^(E)) is discussed in Appendix D. The cross section 

ratios vary < 30% for E > 50ev and will be assumed constant for the interaction of 
auroral electrons in this study (See Figure 1-D). Seaton (1959) eatimated 
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ao+ 2 p(E)/Oo+^^^ to be .2 and Kazaks et. al. (1972) estimated a value of ,25. The total 

O to N 2 ionization cross section (a^^(E)/aj^^(E)) was estimated at .5 from the Kieffer 

and Dunn (1966) review of previous measurements (Appendix D, Figure 1-D). The ratio 
of the total N 2 ionization cross section to the cross section leading to NJB* SjJ in 

the zeroth vibrational level (<r^,+*) was found to be 14.1 by BorsJ and Zipf (1970). 

The state in the zeroth vibrational level emits the X3914A band system. 

Seaton’s (1959) value for ratios discussed above for the 

remaining cross sections yield a value of 1.41 for R^^(E). These historical values are 

theoretically and experimentally verified with the exception of o which 

has been determined theoretically. v-» 1* U 


Another frequently observed emission in aurora is that arising from the N^B* 2 J 

in the first vibrational level (the X4278A band system). The work of Shemansky and 
Broadfoot (1971) suggests the cross section for the zeroth vibrational state is ~ 3.3 times 
larger than for the first vibrational state of NjB*2j. For the cross section factor 

corresponding to the X4278A emission (R^^^ (E)), the above historical values give a value 
of 4.65 for Rp^(E). The production term (Equation 71) would now read 


7?(0+"P)p = R(;j(E) • t?(4278A) • 


(73) 


6.2,4 CONTINUITY SUMMARY OF 0+(*P) 

The continuity equation can be written in terms of cross section ratios, emission 
rates, and quenching rates by setting the production of 0'*'(*P) (Equation 71) equal to 
the loss (Equation 69), i.e.. 


R_ (E) • t? 3914A) • = .219n[0+*P] + Z 7 j n[i] -n[0+*P] . (74) 

^0 nLNj] i •* 


The ©■'■(^P) population can be related to the emission rates resulting from the 
excited population. The population has been related to the total *P-*D volume emissions 
(j 7(7319-30A)) and to the *Pi/ 2 , 3/2 ^^s /3 emissions (t?( 7319A)) in Appendix C, i.e.. 


t?( 7319-30A) = .171 n[0+2p] 


(75) 
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and 


TK7319A) = .098n[O+*Pl 


(76) 


Equations (74) and (7S) can be combined to eliminate the term and solved for the 
volume emission rate ratio, i.e.. 


17(7 3 19-30A) ^ ^ .g. n[OI] 

i 7(3914A) n[Na] 


.171 

.219+ 2 7in[i] 


(77) 


Similarly, Equations (74) and (76) can be combined to give 


r?(7319A) ^ g n[OI] 
t 7(3914A) ntNal 


.098 

.219+ ^ 7in[i] 
i 


(78) 


If the production is written in terms of the 427 8 A emission rate, expressions similar to 
(77) and (78) can be written considering the branching ratio discussion of Appendix D, 
i.e.. 


17(731 9-30A) ^ g nlOll I lAll 

i 7(4278A) n[Nal | .219 + J 7in[i] j 


(79) 


and 


t7(7319A) ^ g .g. nIOIl . .098 

i 7(4278A) ' n[Nal .219+ZTin(i] 


(80) 


To repeat, historically Rj,^(E) “ 1.41, and R^^(E) “ 4.65. Equation (80) will be used 

extensively in the following chapter in interpreting ground based surface brightness ratios 
of X7319A and X4278A emissions. In the following chapter, a rocket payload 
measurement of X7319-30A to X3914A intensity ratios will be interpreted with the aid of 
Equation (77) above. 


63 



I II I III I IIIIIIIIIIIIIIIIIII llllllllll 


CHAPTER VII 

THEORETICAL INTERPRETATION AND MEASUREMENT OF 
0+(2p-*D), 7319A and NtlNEGi,o, 4278A 
EMISSION INTENSITIES IN AURORA 
FROM GROUND STATIONS 


7.1 INTRODUCTION 

A number of measurements were made at College, Alaska, of the intensity of the 
X7319A doublet (I 7319 ) along with the N^lNEGi^o intensity (I 4278 ). Expression (80) 
of the last chapter relates the volume emission rates of these emissions. This expression 
will prove to be valuable in this chapter when interpreting the .ground based data. 

In the theoretical development of the previous chapter, it was assumed that the 
primary electrons were mainly responsible for the NJ ionization rate and resultant 
NJINEG emissions. In interpreting the total volume emitting layer, in a ground based 
measurement where one is oriented to measure along the flux tubes, a significant 
contribution to the total intensity comes from the lower scale height of the atmospheric 
range of the precipitating particles. Significant energy degradation is beginning to take 
place in the primary electron energy distribution at these altitudes. One of the intentions 
here is to ascertain where the electrons are ranging into the atmosphere. Consequently, 
the description of this production function shape is pertinent to this interpretation. 
Berger et. al. (1970) and Banks et. al. (1974) have calculated theoretical profiles of 
production in aurora for monoenergetic, isotropic electron beams in the atmosphere. The 
particle range indicative of these investigations and spectroscopic studies by Rees and 
Luckey (1974) will be utilized in the interpretation of these ground based measurements. 

Previous theories of 0 '*’(*P) quenching and branching will be presented and 
compared to the ground based observations. 

The refinement of the theory and appUed methods describing the C^f^P) and 
Nj INEG emission relationship in the context of interpreting ground based measurements 
is discussed near the end of the chapter. 

7.2 GROUND BASED OBSERVATIONS OF I 7319 AND RELATED EMISSIONS 
7.2.1 INTRODUCTION TO THE GROUND BASED MEASUREMENTS 

The multi-spiral spectrophotometer (discussed in Chapter V) was set up and 
operational between February 2, 1973, and February 24, 1973, at the University of 
Alaska, Geophysical Institute (Ester Dome). The instrument was mounted on the roof of 
Ester Dome Observatory in a magnetic meridian orientation with a range of ±60® from 
the zenith. 
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The University of Michigan Airglow Observatory was located adjacent to Ester 
Dome. The Michigan facility was equipped with a four channel, tilting fUter photometer 
which monitored the OI, 6300A, and 5577A, as well as H^, 6S63A, and NtlNEGi,o» 

4278A emissions. The orientation of the Michigan photometer and the sector 
spectrophotometer (set up for 7319A) were position aligned by imaging the moon 
through the instruments for calculated moon zenith positions. A time reference was 
frequently coordinated between the two systems. The cone half angle of the two systems 
was 2 ®. 

7.2.2 OBSERVATIONAL PROCEDURE AND RAW DATA SAMPLES 

The primary mode of operation was lo keep both instruments oriented to the 
magnetic zenith. Arc scans would have been useful provided photometer triangulation 
scans for arc location were made simultaneously. When weather, moon conditions, and 
aurora cooperated for such an observation, equipment did not. This type of scan would 
be desirable for future measurements. 

A strong band system (NjlPGsjs) induces a varying background in the 
0‘*‘(*P-*D) wavelength region (Chamberlain, 1961). The ground based observations are 
also contaminated by OH(8,3) band emissions originating from the upper mesospheric 
region. Specifically, the Pi ,j(7316.4A) and Pj,j(73()4.5A) fell within the scanning range 
of the instrument. 

A monochromator was placed inside the spectrophotometer dome and calibrations 
were run after each night’s sky data was gathered. A calibration of February 8 , 1973, is 
shown in Figure 24. Spectrophotometer scans of the monochromator source for 7319A 
and the two OH positions are shown in the figure. 

The hydroxyl was always present. Figure 25 shows a typical series of scans with 
aurora absent (i.e., determined by lack of I 4278 emission for this same time period). 
Note in Figure 25 that the instrument did not totally separate or resolve the OH at 
7316.4A from the 0 "*'(*P-*D) at 7319A. The data reduction scheme used to discern the 
0 +( 2 p. 2 o) will follow in the data reduction section. 

A typical 7319A enhancement with an overhead auroral surge is shown in Figure 
26. Scan one of this figure typifies a hydroxyl only type scan. In the succeeding scans, 
the peak was shifted to a higher wavelength than indicated in spectra one. This indicates 
the 7319A was enhanced and not the OH. The scan rate in this case was ~ 1 per minute. 

7.2.3 GROUND BASED REDUCED DATA SAMPLES 

The deduction of the I 7319 and 17318.4 (OH) intensities from the instrument 
scan involved an algebraic deconvolution of the effective instrument transmission at the 
two wavelength positions of emission. This procedure was necessary since the instrument 
did not totally resolve the X7319A emission from the X7316.4. The details involving the 
analytic technique for reducing the respective intensity values are given in Appendix E. 
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Figure 24. Spectrophotometer scans of a monochromator source tuned to wavelength 
positions 7319A, 7316.4A and 7304.5A on Feb. 8, 1973, after the Feb. 7/8 
sky data scans at Ester Dome, Alaska. 



Degrees 


Feb. 6/7, 1973 
Mognetic Zenith 
A • Aperture Position 



Figure 25. Spectrophotometer scans of the sky on Feb. 6/7 at Ester Dome, Alaska. No 
aurora was evident and these scans typify OH presence in the night sky. The main 
feature corresponds to OH P, j at 7316.4A. 



Figure 26. Spectrophotometer scans of the sky (magnetic zenith) on Feb. 7/8, 1973, 
at Ester Dome. This series of scans typify the passage of an auroral surge passing 
through the magnetic zenith with ~ 35 rayleighs 7319A present in scans 2 and 3. 
The flag over each spectra indicates the 7310A(BG), 7316.4A, and 7319A 
wavelength positions by vertical hashmarks. Note that the peak shifts 
towards 7319A wavelength with increased signal amplitude. 
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Reduced intensities of 7316.4A and 7319A for a two hour period during an 
auroral breakup on February 22/23, 1973, are shown in Figures 27 and 28 respectively. 
The scatter in the OH Pj ,2 data is ~ 8 rayleighs with a mean of ~ 14 rayleighs for the 
two hour period. The scatter increases during periods of enhanced aurora and resultant 
7319A bursts depicted in Figure 28. 

Figure 29 shows a plot of the time history of the intensities of X4278A and 
X6300A for the same time period as the previous two figures. Figure 30 represents the 
ratio of the two emissions displayed in Figure 29. 

A cursory examination of the 7319A and 6300A data shows the 6300A emission 
enhanced over nightglow intensities between events but the 7319A data is not. The 
7319A data is much more bursting in nature. There is, however, a high event correlation 
between the 6300A, 7319A, and 4278A data. A closer look at Figure 30 suggests a ratio 
of I 6300 A 4278 value of .5 prior to 2325 A.S.T. and after 0025 A.S.T. with a ratio of 1. 
for the hour in between. These values correspond to those of Bather and Mende (1971) 
for the oval region. There also appears to be a slight ratio decrease trend in the period 
just before midnight. At eight minutes after the hour, note that a major event occurred 
in all data. 

Figure 31 further suggests the implication made earlier that I 7319 /I 6300 is very 
erratic. The 6300A data has a large continuous background, as well as the 4278A data, 
but 7319A does not and it is further reflected in the variation of I 7 3 1 9 / 1^300 ii^ this 
plot. 


Another case of 7319A and 4278A history plots is displayed for February 7/8, 
1973 in Figure 32. Here also noted is a high event correlation and again a major event 
very near midnight with a large amplitude in 7319A emission. Note the scales differ by 
an order of magnitude for the two emissions on the figure. A total of seventeen nights of 
data was taken with about half of them being fairly active. The February 7/8 and 
February 22/23, 1973, cases represent some of the more active periods. 

7.3 INTERPRETATION OF GROUND BASED I7319/I4278 DATA 

The approach used to interpret the ground based measurements of I 7 3 1 9 /I4 2 7 8 
ratios takes advantage of what is already known about the I 6300 /I 4278 ratio. 

Rees and Luckey (1974) have calculated the I 6300 /I 4278 ratios for particular 
monoenergetic, isotropic electron beams in the thermosphere. Berger et. al. (1970) and 
Banks et. al. (1974) have computed NJ ionization profiles for similar monoenergetic 
isotropic electron beams. The altitudes of maximum ionization for given energies from 
Berger et. al. (1970) and Banks et. al. (1974) are point plotted in Figure 33 and a hand 
drawn fit is drawn through the points. 
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Time lA.S.T.) 

Figure 28. 1 , 3 1 , versus time for a two hour period near midnight on Feb. 22/23, 1973. 
The data points have been straight line segmented to bring out the bursting 

nature of the emission. 
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Figure 31 . Ratio of 17319/16300 for 22 / 23 , 1973 . 4300 was time interpolated 
to correspond to the time of the I7 3 1 9 measuremetit for ratio calculations. 
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Figure 32. I 4278 I 7319 versus time for a two hour period near midnight 

on Feb. 7/8, 1973. The orientation of the instruments was the 

magnetic zenith. 
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Figure 33. Altitude of ionization maximum versus initial electron energy 
for monoenergetic, isotropic electron beams. The points are 
provided by Berger et. al. (1970) and Banks et. al. (1974). 

The curve through the points is a hand fit to the points 
on the tigure. 



Rees and Luckey (1974) values of I 6300 /I 4278 are plotted in Figure 34 (Curve 
A) where the energy to altitude conversion was made using the range curve plotted in 
Figure 33 which was discussed above. The interpretation then is that a given intensity 
ratio of l6 3 ooA 427 8 as measured along the magnetic field should have an ionization 
maximum at the corresponding altitude for a monoenergetic, isotropic electron beam. 

A plot of the I7 3 1 9 /I4 a 7 8 ratio can also be generated using the information 
discussed in the previous chapter. The intensity of the respective emissions can be 
expressed as the integrated volume emission rates, i.e., 


00 

I7319 = / T?(7319A)dz . ( 81 ) 

z 


Equation (80) of the pirevious chapter can now be expressed in an integral form as 


i”,» - *0.® / "ton 


.098 


.219+ ^ 7in[i] 


dz 


(82) 


The integral of Equation (82) can be further simplified with some reasonable 
assumptions. In the derivation of the Nj production by electron impact in Appendix D, 
the production to density relationship was related to an energy integral (Equation 2-D). 
Further assuming that 7 j(NJ)« tj( 4278A) (Borst and Zipf, 1970), expression (2-D) can be 
written as 


F(E)a +(E)dE 
nlNj] ‘'e ^2 


(83) 


where C is a proportionality constant. 

In the interpretation of the U 300 / 1437 8 intensities Rees and Luckey (1974) used 
exponential functions (discussed in the introduction to Chapter VI) in the description of 
F(E) for auroral electrons. It was assumed that the Rees and Luckey (1974) assumption 
is very near the shapes described by Berger et. al. (1970) and banks et. al. (1974) which 
were used as a source for the altitude of maximum ionization indicated by the Rees and 
Luckey (1974) intensity ratios. 
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ALTITUDE (Km) 



Figure 34 . Altitude versus Ix/Usts where X refers to 6300 A emission 

or 73 19 A emission. A theory for Issoo/Uavt ^wo theories 
for 17319/14 37 • are plotted on the figure and discussed in the 
text. Ground based measurements of 17319/1437$ <ure plotted 
for diffuse aurora and for **disctete events’*. 
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In Appendix F, a broad range of Nf production shapes is plotted and compared. 
It was deduced that the Nj proportional profile shape, as described in Appendix E, 
closely parallels the exponential shapes of the previous models, particularly for the very 
energetic beams ranging below 200 km. The Nj proportional shape assumes that the 
volume production of the emissions arising from the NJ excited states is constant to the 
range altitude of the electrons, and zero below. This assumption enables the ratio 
described in Equation (84) to be specified as constant over all altitudes in the electron 
precipitation regions, and furthermore, may be expressed in the integral form, i.e.. 


OO 

t?(4278A) ^ I r?(4278A)dz 

n[Na] «> 

/ n[Nil dz 
z 


(84) 


For an exponential atmosphere, (84) becomes 


t?(4278A) ^ 
n[Nj] 


142 7 8 
nlNj] • 


Z 


(85) 


where is the nitrogen scale height. 

Equations (83) and (85) can now be combined to give the expression for the 
^7319 /I 4 278 r^tiO 


^73 19 
^42 78 


Rg,(E) 

n[Nj] 


00 

/ n[OI] 
z 


.098 

.219 + E 7i n[i] 


dz 


( 86 ) 


Two theoretical profiles of I 7 3 1 9 /I 4 2 7 s (Curves B and C of Figure 34) were 
generated using Equation (86) above. They represent an ‘Old Theory’ (Curve B) and a 
‘Recent Theory’ (Curve C). 

The ‘Old Theory’ assumed values of = 1 • 10"‘° (Dalgamo and McElroy, 
1965), 7 q = 3 • 10‘® (Seaton and Osterbrock, 1957), and a p(E)/oQ+(E) of .2 after 

Seaton (1959) (i.e., R^^(E) = 4.65). The neutral composition arid temperature model of 
Jacchia (1971) for To„ = 1000°K was used for this curve and all theoretical curves which 
follow in this study. A value of 10® cm'^ was assumed for n[e] . 
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The ‘Recent Theory’ (Curve C of Figure 34 ) was generated identical to that above 
except for revised quencldng coefficients. Walker et. al. ( 1974 ) values of 7 q“ 2 ' 10 "‘®, 
and 7^^ = 5 ’ 10 '“, replace the Nj quenching of the ‘Old Model’ and the 7g = 2 d 0 "’ 

value of Henry et. al. (1969) replaces the ‘Old Theory’ above. 

The curves plotted in Figure 34 will assist in assessing measured ratios of 
Ij3oo/l43 7 8 ai'd I73i9/l4378> ®nd Luckey ( 1974 ) have carefully worked out the 
ratios of the classical auroral lines included in Curve A of Figure 34 . It was assumed that 
Curv'e A was correct and, for a measurement of l33oo/l4378> f^e altitude of maximum 
ionization was inferred from the height relationship of Figure 33 . As discussed earlier, the 
ground based measurements included I« 300/^43 7 8 measurements as well as 17319/14373. 
The measured 17319/14378 was then plotted on Figure 34 at the altitude which 
corresponded to that indicated by the I« 3 00/14378 discussed above. This was 

accomplished for two types of aurora discussed below. 

Time histories of I731 9/I4378 ratios for all data points are plotted in Figures 35 
and 36 for the two nights’ data sample previously discussed. ‘These data points include 
the ‘Diffuse Aurora’ as well as discrete events. A hand fit mean (solid curve) and 
extremes (dashed curve) on each of the two figures indicate the ratio tendencies for the 
indicated time periods. Figure 30 shows the time changes in the observed I6300/I4378 
ratio for the night of February 22 / 23 , 1973 . In Figure 30 , it is noted that the 
I3300/I4378 ratio remains nearly constant for approximately 15 minute time periods. 
I3300/I4378 ratios were selected from Figure 30 at ~ 15 minute intervals and the mean and 
range of the 17319/14378 were determined from Figure 35 for the same local time. The 
mean I7319/I4378 the range are indicated in Figure 34 by a point and range bar 
respectively, as described in the previous paragraph. The same was done. for the February 
7/8 which has also been plotted on Figure 34 . 

Discrete events are also plotted on Figure 34 (triangles). In these cases, the diffuse 
aurora component was removed from both the I4300 1437s measurements and the 

ratios were taken at their maximum value as that maximum occurred in the field of view 
of the instrument. The I7 3 1 9 showed little or no background as shown in Figure 28 . 

The altitude coverage of the ground based measurements was not large, as 
indicated in Figure 34 . The I63oo/l4 37 8 ratio ranged from ~ .5 to ~ 1.0 on both nights 
indicating that most of the aurora was ranging to the 120 to 130 km region. It is 
apparent that the theoretical predictions of the I7 3 1 9 /I4 2 7 s ratio are larger than the 
measurements. The predicted ratio of Curve C (Walker et. al., 1975 , quenching reference) 
suggests a value of .13 for I7319/I4278, whereas the measurements suggest a ratio of . 06 , 
a factor of 2 difference. The ‘Old Theory’ (Curve B) predicts ratios much larger than 
those measured at these altitudes, on the order of a factor of 5 difference. 

The picture is expanded in the next chapter by adding ratios measured at higher 
altitudes. Further theory and measurement inconsistencies will be discussed following the 
addition of this other data. 
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CHAPTER VIII 

Na1PGs,3, NjlNEGo,o, AND 0+(*P-*D) INTENSITY 
PROFILES IN AN AURORA AND DEDUCED 0*{^P) 
ELECTRON IMPACT PRODUCTION EFFICIENCY 
AND QUENCHING RATE \ 


8.1 INTRODUCTION 

Described in this chapter are the data, data analysis and geophysical interpretation 
of 0'*'(*P-^D), N^lNEGo,o» and NalPGs.a emission measurements made on a rocket 
payload in an active aurora. 

In the introduction to Chapter VI, some discussion was given on the difficulty in 
measuring the Q^C^P-^D) emission. The 7319A, 7330A emissions of the metastable ion 
are surrounded by NjlPG emissions. Chamberlain (1961) has published auroral ground 
based spectra of this wavelength region. The NjlPGs ,3 and NjlPGg ,4 bands are both 
overlapping the 0''‘(*P-^D) emission. The ground based spectra also shows nearby OH 
emissions which originate in the mesosphere. The OH contaminants do not need to be 
considered on a rocket payload located above the OH emission altitudes. 

The next section of this chapter presents the rocket photometer data in its 
measured form and also includes the electron density profile as deduced from a Langmuir 
probe on the payload. The three photometer channels shown were bandpass peaked at 
3914A, 7326A, and 7354A. The data analysis section (Appendix G) describes how the 
rocket photometer outputs of the previous section were used to deduce the amount of 
respective NjlNEGo,©. N 3 lPGs, 3 , and 0 "‘’(*P**D) emissions present in the field of view 
of the rocket photometers.The technique used in deducing the amoimt of 0'*'(® P-* D) and 
N 3 lPG 5,3 was somewhat complex. The amount of NjlPG passed by the 7326A and 
7354A channels required an N 3 IPG band model to synthesize the band shape. The 
N 3 lPGs ,3 and 0 '*'(*P-*D) emission intensities were then deconvolved from the 7326A 
and 7354A channel currents. 

One of the outcomes of the analysis was the intensity ratio of the N 2 lPGs ,3 and 
NJ1 NEGo,o l>and systems. This ratio is discussed next before proceeding into the main 
outcome of the measurement, the 0 '^(*P-^D)/N 2 lNEGo,o intensity ratio. 

The appropriate emission branching ratios were applied to the 0‘*'(*P-*D) and 
NJ1 NEGo,o measured emissions to convert to the I 7319 /I 4 J 78 ‘measurement’. This was 
done so that the data of the last chapter might be compared and plotted with this data. 
The rocket payload and groimd based I 7319 /I 4378 data were found in good agreement 
and the ‘Old Theory* and ‘Recent Theory’ discussed in the last chapter did not describe 
the measurement. 

Two approaches were taken to modify the theory such that appropriate 
quenching and cross section ratios would in turn describe the measurements. One 
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approach was to reduce the quenching previous investigators had applied to 0"'’(^P), as 
well as, the rate at which 0'*'(^P) was produced by electron impact. The second approach 
was to accept large quenching and postulate an additional low altitude source of 0'*'(*P). 
The chapter concludes with some limits on the manipulated parameters. 

8.2 NASA 4.329 AURORAL MEASUREMENTS 

On March 39, 1972, NASA 4.329 was launched into an active aurora from Fort 
Churchill, Canada (58.7°N, 93.8°W). The University of Michigan Airglow Payload was 
equipped with a turret photometer, particle detector, and Langmuir probe from which 
the following analysis was made, (Sharp and Hays, 1972, and Sharp and Hays, 1974). 

The particle detector scanned from .5ev to ~ SOOev electrons. Particles were first 
detected at 1 20 km. This altitude was used to discern that the payload was in the aurora 
above 120 km and prior (or below) this altitude, the emissions observed from the 
photometers were that of the aurora above the payload. 

The photometer channels (3914A, 7326A, and 7354A) were broadband (AXj/j en 
20 A to 30A) and the channel voltage profiles are plotted in Figures 37, 38, and 39 for 
the three channels of interest in this study. The circles and dots shown on Figures 38 and 
39 represent measurements using a wide angle and narrow angle field of view. The curve 
on the two figures represents a hand fit of the data. This was done to interpolate values 
at 2 km increments over the range of data as required for the data analysis. The 
sensitivities were determined with a C, 4 low brigluness source calibrated by the National 
Bureau of Standards (Appendix G). The detector dark current has been subtracted from 
these profiles prior to making these plots. 

Finally, the electron density, as deduced from the Langmuir probe, is plotted 
versus altitude in Figure 40. This data is required for the theoretical application of 
electron quenching in analyzing the continuity equation of O'^'C^P). 

8.3 REDUCED AURORAL EMISSION FROM NASA 4.329 

The data reduction procedure involved the laborious task of deconvolving the 
0 '*’(^P-^D) emission signal and the contribution from the neighboring N 2 IPG band 
emissions from the two photometer channels peaked at 7326A and 7354A. This 
procedure is described in Appendix G and the final results are plotted in Figure 41. The 
curve for the 0"*'(^P-^D) emission represents a hand fit curve through the data points. 

The rocket payload passed through an enchanced region beginning at ~ 140 km 
altitude. It is uncertain as to whether a new discrete event formed while the payload was 
ascending or whether the payload moved through an existing event. 




Figure 38. Altitude versus channel voltage for the 7326A channel for NASA 4.329. 
The circles and dots represent measurements using a wide angle and narrow 
angle field of view, respectively. 









Intensity 


Figure 41 . Altitude versus intensity of P- 

as measured on NASi' 






8.4 A DISCUSSION OF THE I1PG AND I1NEG INTENSITY RATIOS FROM NASA 

4.329 

Before going into the analysis of 0'*'(^P) continuity theory and deduction of the 
quenching and production efficiency coefficients, a brief discussion of the nitrogen 
emissions will follow. 

Plotted in Figure 42 is the intensity ratio observed on NASA 4.329 which is 
shown to change very little over the altitude range of measurement. The measured data 
(point plots) were hand fit with a straight line (dashed curve) which is drawn at a ratio 
of .34. The data deviated very little from this and the scatter seen is most probably 
measurement noise rather than actual ratio changes. 

A measurement summary of this ratio was discussed by Shemansky et. al. (1972) 
with respect to energy distributions of auroral electrons. The author speculates that the 
reason previous investigators have observed this ratio to change with altitude in aurora is 
due to misinterpretation of Nj IPG photometer data. The temperature dependence of the 
band shape and the manner in which band shape alters effective filter transmission must 
be considered for proper interpretation of this photometer measurement. 

Examination of electron impact ionization cross sections leading to the N^INEG 
emissions (as reviewed by Kieffer and Dunn, 1966, and Green et. al., 1973) show 
threshold energies of ~ 19ev with maximum ionization efficiency at energies ~ lOOev. 
The threshold for creating the state of N2 is ~ 7ev and the excitation cross section 

peaks at ~ 20ev (Green et. al., 1973). 

Opal et. al. (1971) have studied the energy distribution and angular dependence 
of secondary electrons generated by the impact of 100-200ev electrons for a number of 
gases including N2 . They found that the average energy of the secondaries for the impact 
energy range used was ~ 13ev and the distribution was fairly insensitive to the energy of 
the impact electrons (primaries). 

The number of secondaries is of course directly related to the number of electron 
impact ionizations on the ambient amtospheric constituents. The number of N^lNEGojO 
photons is directly proportional to the N2 ionization rate and the number of secondaries 
produced. Since the energy distribution of the secondaries doesn’t vary with the energy 
of the primaries and presuming the secondaries are the source of N2B^7Tg state, the 

N2IPG emission intensity should be directly proportional to NJ1 NEGo,o aS observed 
here. 

8.5 COMPARISON OF NASA 4.329 DATA WITH PREVIOUS THEORY 

The intensity measurements versus altitude provide another data source to add to 
the picture begun in the last chapter with the ground based measurements. In Figure 34 
of the previous chapter there were plots of two theoretical I7319/I 4278 versus altitude 
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Figure 42. Ratio of IlPG 5 , 3 /llNEGo,o versus altitude for NASA 4.329. 





and ground based data. NASA 4.329 measured 17319.30 and I 3914 values as a function 
of altitude. The appropriate branching ratios (discussed in the last chapter and Appendix 
C) were applied to NASA 4.329 data to interpret the I 7319 /I 4278 ‘measured’ ratio. The 
x’s in Figure 43 represent this new addition of measured values. This rocket payload 
measurement agrees quite well with the groimd based measurements but deviates 
significantly from the theoretical curves. Furthermore, the error analysis discussed in 
Appendix G suggests the error in the measurement is least at the higher altitudes where 
the theoretical ratio predictions are in large disagreement with the measurement. 

8.6 THEORETICAL EXPLANATIONS OF THE AURORAL I73 19/4278 

MEASUREMENTS 

The basic formulation for generating the theoretical ratio curves of Figure 43 was 
Equation (86) as developed in the previous chapter i.e.. 


^73 1 9 

I4278 


nINJ 


00 

/ n[OI] 
z 


.098 

.219+ 2 Tin[i] 


dz 


( 86 ) 


where the values for parameters previously discussed were taken from historical findings 
and a model atmosphere. To review. Curves B and C of Figure 43 were generated using 
values of 3-1 O'® for 7 g and MO'*® for for Curve B and 2-10'’ for 7 g, 2-10'*® for 

7 Qand 5-10'* * for 7 jq^ for Curve C. 

Fundamentally, there are two basic disagreements between ‘old’ and ‘recent’ 
theories and measurements depicted in Figure 43. The theoretical shapes of the ratios 
with altitude suggest a larger change or larger slope over the altitude range than the data. 
Examination of (86) suggests that the ratio changes as a result of two possibilities. If all 
7 j were zero, the integral would reduce tonlOI]^- Hq, thus setting a lower limit on the 

altitude dependence of the ratio to approximately the N[OI]/n[N 2 l ratio. Finite values 
for 7 j tend to increase the ratio change with altitude. 

Secondly the magnitude of the ratio predicted by theory is much larger than 
measured, particularly at the upper altitudes where the measurement is least in error 
(Table 2-G). Assuming that the error in the emission probabilities and model atmosphere 
parameter errors are small, the (E) parameter (cross section ratios) requires 

examination. Of the product ratios involved in this factor (Chapter 6), only one, 
is not both theoretically and experimentally verified. A departure from 

the Seaton (1959) ratio of .2 for this branching ratio will be necessary to theoretically 
predict the lower ratios of I 73 1 9 /I 4 J 78 observed in these measurements. 
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Figure 43 . Altitude versus 17319/14378- The plot includes two previous theories 
(discussed with Figure 34 ) two new theory curves and the ground based 
and rocket measurements of I7 3 1 9 /I4 2 7 g • The particulars of the 
figure are discussed in the text. 
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A pictorial evolution of several theoretical possibilities will begin with the NASA 
4,329 measured ratio of Figure 43 replotted in Figure 44. The solid line through the data 
(D) is a hand drawn fit of the data and the error flags represent the non-systematic error 
of I 7 3 1 9 described below. The error bar reflects the assumption that all the uncertainty 
in the ratio lies in the uncertainty in the I 7 3 1 9 measurement. The non-systematic errors 
were a result of the root summed square of the statistical error and the error component 
in the first four columns of Table 2-G, which exceeded those systematic errors at the top 
of the profile. The filter drift error, for instance, was not treated as a non-systematic 
error. It was assumed that if . the filter peak drifted between calibration and the 
measurement time, the filter stayed shifted for the measurement period and was 
considered systematic to the entire profile. The calibration error, however, varied from 
±36% at 120 km to ±24% at 160 km. If the measurement was a single measurement and 
not a deconvolution, these errors would be constant, or also systematic with altitude. The 
component of 12 % at 120 km for the calibration table then was also treated as a total 
shape error component or non-systematic error. 

It should be noted at this point, also, that all the theoretical synthesis from this 
point on involved a synthesis of the measured profile of I 731 9 /I 4278 from NASA 4.329. 
In the theoretical treatment, the measured electron density (Figure 40) was used rather 
than the MO"^^ cm'^ assumed in the last chapter as it applied to Equation 86 . 

8.6.1 THEORY 1, THEORY OF 'LOW QUENCHING' 

The first approach was to ignore early theories as to the quenching rate of 
©■^(^P), as well as the a (E)/a .(E) branching ratio. Instead, these parameters were 

Q-rz p QT 

iterated until a combination was found which would put theory and measurement in 
agreement. It was felt by the author that lower limits on quenching should exist for ‘well 
established’ or weU understood processes. It was assumed that 7 ^ could not be less than 

3- 10'® and that 7 ^^^ could not be less than 2-10‘‘*. The minimum was postulated since 

the additional internal energy associated with ©"^(^P) would not decelerate the rate of 
dissociative recombinate over the established 0'*"(‘*S) rate of 2 - 10'*^ (Chapter VI), 

The first step involves taking the theoretical curve B of Figure 43 and reducing 
the branching ratio until it was normalized to the data curve at 170 

km. This is shown as curve E on Figure 44 and is not within the confines of the error of 
the measurement. The branching ratio was reduced from the .2 value of Seaton to ,018 
for this normalization. 

The next step was to normalize the curve to the extreme ratio upper limit at the 
upper altitude (170 km) and to decrease the Nj quenching until the lower part of the 
curve was in the error limit of the data. This profile was generated by an Nj quenching 
rate of 1-10"* * and a branching ratio of .010 ( 7 g was assumed at 3-10'®), This is plotted 

as curve F in Figure 44. By using a value for 7 ^^^ of 240'*^, 7 g of 3‘10‘®, and a 

branching ratio of .008, the theoretical profile becomes almost identical to that of Curve 
D, the hand fit data curve. 
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ALTITUDE (Km) 


CURVE D, NASA 4.329, DATA & ERROR 

+ NASA 4.329 DATA POINTS 
THEORY. 

ASSUMED BRANCHING RATIO AND 
QUENCHING 

CURVE [ gp*2p/ °o* I TN2 I T> [ 


Figure 44. Altitude versus 17319/14278' Th® data from NASA 4.329 and two 
theoretical curves are plotted showing the effects of altering 
<^Q+2p(^y^Q+(^^ and tn of early theories. 






In summarizing, the above analysis suggests that the error limits on the data curve 
confine or define the following: 


Cj,«y(E)/a^+(E) - ,008 

> .005, < .10 
Te - 3 • 10-* 

7N, = 2 ■ 10->» 

< 1 • 10‘‘ ‘ 


Most probable 
Extremes 
Assumed 
Most probable 
Extreme 


The most probable values above were applied to Equation (86) to generate Curve N of 
Figure 43, which represents one possible combination of parameters in the theory that 
describes the data. 

The above combination is one possibility, but a recent analysis of data reported 
by Walker et. al. (1975) strongly suggests that large quenching is required of 0‘*‘(*P) to 
explain intensity measurements of 7319A emission on Atmospheric Explorer Satellite in 
dayglow. Examination of curve E on Figure 44 which represents large quenching by Nj 
does not fall in the error bars of the measurement. The only remaining possibility is that 
Equation (86) is not complete and that there possibly exists another production 
mechanism for 0‘*'(*P) that has not been considered. 

8.6.2 THEORY 2, LARGE QUENCHING AND N-^-f'S) PRODUCTION OF O'^(^P) 

Many other possible production schemes were considered, such as direct 
dissociation of the Oj by electron impact, but only one appears to be potentially large 
enough at low altitude to overcome large collisional deactivation. Torr (1975) suggested 
that 


N+(‘S) + 0(3P) -»• N + 0+(*P) 


(87) 


could possibly occur at a near kinetic rate. It was assumed here that for every N'^(‘ S) 
produced, an 0"‘‘(*P) was also produced. 

At what rate might N"'‘(‘S) be produced in aurora? Jones and Rees (1973) give a 
rate of N"*" production as related to the volume production rate as: 


»?(N+)p 


Nt 


7?(nJ) 


( 88 ) 
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where aj^^(E)/aj^^(E) was estimated at .25. Next, the volume production rate of NJ. is 
related to the X3914A volume emission rate by 


n(Nt)p 


"nJ.® 


»7(3914A) 


(89) 


(See chapter VI). To the author’s knowledge, the branching ratio for producing N‘*'(‘ S) 
to N"*" is unknown (i.e., g(E)/Oj^^(E)). The aj^^(E)/aj^^^(E) was discussed in 

Chapter VI. To further simplify. Equations (88), (89) and the N‘‘’(* S) to N"*" branching 
ratios may be combined to give: 




It was further assumed above that N‘*'(‘S) immediately reacted to produce an 0'*‘(^P). 
The new production source of 0'*’(*P) may be written, using historical values for the 
known branching ratios, as 


T?(0+*P)p = C^,^(E) • 77(3914A) 


(91) 


or as 


n(0+*P)p = C^^(E) • n(4278A) 


(92) 


where C^(E) is the cross section ratio product of Equation (90), where the subscript 
refers to the upper vibrational level of S“. 

Next is a cursory comparison of the combined production terms (i.e., the charge 
exchange process described above and the electron impact source treated earlier). The 
total production term may be obtained by summing the above source (Equation 92) with 
the electron impact source (Equation 73) to read 


Tj(0+2p)p = R^^(E) tj(4278A) + Cj,j(E)tj(4278) . (93) 
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In the model atmosphere used in this study, the nEOIl/nlNj] ratio is ~ 1 at 160 km. 

For the historical value of 4.65 for (E) a branching ratio of .4 for a . , _(E)/a .(E) 

* N *S Iv 

would make the two production sources equal. The other cross section ratios making up 
C^^(E) have been previously discussed. Since nIOIl/nINj] increases with altitude, the 

charge exchange mechanism becomes even more dominant at lower altitudes. 

Another way of writing (93) involves isolating the volume emission rate, i.e., 


.l(0+=P)p = ^R„_(E) + C„_(E)^ „(4278A) 


(94) 


where 


C„_(B) .. >1.5 


(95) 


The continuity equation of 0’*‘(*P) in Chapter VI can now be rewritten using 
Equation (94) above for the production equation, i.e., 

^R„^(E) + C^/E)^ t?( 4278A) = .219n[0+2p] + 2 7in[i] n[0+*P] . 


(96) 


This reduces to a volume emission ratio of the form 


t?( 7319A) _ L .p. n[OI] 
r?(4278A) \ ^ nENj] 
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.098 

.219+ 2 7in[i] 


(97) 


or in the integral form as: 
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= __.£i f n[oi] 
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(98) 
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Expressions (97) and (98) follow the identical derivations (80) and (86) did in Chapters 
VI and VII. 


Next, the quenching rates of Walker et. al. (1975) were assumed correct. First, it 
was assumed that 7q= 2‘10‘“* and = 5-10'‘^. It was also assumed that 7 g = 2 •10"'' 

(Henry et. al., 1969). The uncertainties in the two branching ratios, i.e., 
0 'o+ 2 p(E)/Oo+(®) in R(;j(E) and g(E)/Oj^^.(E) in C(^j(E) remain as variables. 

Plotted in Figure 45 is the hand drawn data curve and error bars as described for 
Figure 44. Also plotted are three theoretical curves using the quenching rates of Walker 
et. al. (1975) and Henry et. al. (1969) where the N‘‘'(‘S) branching ratio was assumed as 
.05, .10, and .20 for G®, G*®, and G*® respectively. For each assumed N'*'(‘S) branching 
ratio, an 0'*'(*P) branching ratio was computed such that the theoretical curve was 
normalized to the data at 170 km. The required 0‘*‘(*P) branching ratio for each curve is 
indicated on the curve table in the figure. It is noted that if curve G® were shifted to the 
right error limit at 160 km, the lower portion would fall at the left error limit at the 120 
km altitude. A branching ratio of .05 for g(E)/Oj^^(E) is therefore a minimum 

acceptable rate for this quenching combination according to these limits. It should be 
added that a branching ratio of .24 for N^(‘S) requires that the <^Q+ 2 p(E)/o^+(E) 

branching ratio is zero for 170 km synthesis of the measured I 7319 /I 4278 ratio. Thus, 
an upper limit of g(E)/o^^(E) of .24 is established for this quenching combination. 

Curve G*®, which closely parallels the data, has been plotted in Figure 43 for overall 
measurement comparison. 


Walker et. al. (1975) also suggested that a value for 7 j»^^ of 5-1 0'^® and zero OI 

quenching would describe their data. A similar analysis was accomplished for the above 
quenching combination as was done for the previous combination. The results are plotted 
in Figure 46. Analogous to the previous discussion, only Curve H^® can be shifted to faU 

in the error limits of the data, establishing a lower limit of o . , (E)/a .(E) at .20 for 

S 

this quenching combination. A value of .23 for a , (E)/a .(E) required 

N'*’* S N"^ 


aQ.,. 2 p(E)/ao+(E) to be zero at 170 km synthesis of the 17319/14278 ratio. The upper 


limit for a . , „(E)/o^ .(E) is therefore .23. 
N+*S N+ 


8.7 DISCUSSION OF THEORETICAL COMBINATIONS 


In the previous sections, discussion has surrounded three possible quenching 
combinations, all of which fall within the error limit of the measurements. In each case 
there are limits on two branching ratios which have not been previously verified by 
experimental measurement. The extremes and most probable values for these branching 
ratios are tabulated in Table 5. 
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ALTITUDE (Km) 


1 1 1 I f I I I I 

CURVE D, NASA 4.329 DATA & ERROR 

THEORY. 7 o = 2.10-10. WALKER et.al. 1975 

7M« = 5.10—11 


ASSUMED BRANCHING RATIOS 
«o+2„/oo+ 


Figure 45 . Altitude versus I7 3 1 9 /I4 2 7 8 • Same as Figure 44 except the 
theoretical curves show the effect of N'*’(* S) production efficiency 
for the “Large Quenching, OI” model. 





Figure 46. Same as Figure 45 except the theoretical curves show the effect 
of S) production efficiency for the “Large Quenching, N 2 ” model. 
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TABLE 5. SUMMARY OF THEORETICAL COMBINATIONS DESCRIBING 
0+(2p) CONTINUITY IN AURORA 



V-s® 

Max 

Most 

Probable 

Min 

<^0+2 p(E) 

Max 

Most 

Probable 

Min 

Theory 

V® 

V® 

Low Quenching 

.04 

.010 


0 

.008 

7e = 3 -10-® 

(assumed) 

0 

.005 

High Quenching, 01 






70= 

.24 


.030 

7N, = 

.24 

0 

7e = 2 • 10-’ 

Walker et. al. (1975) and 
Henry et. al. (1969) 

.05 

0 

High Quenching, Nj 






7n, = 5 10-‘» 

.23 


.003 

7e = 2 • 10-’ 

.23 

0 

Walker et. al. (1975) and 
Henry et. al. (1969) 

.20 

0 


102 




The low quenching case essentially did describe the data with a single production 
source by electron impact (See Curve N on Figure 43). A very low branching ratio for 
N"*'(‘S) produced an additional source of 0'*'(^P) at low altitude beyond the error limit. 
The limits of <^Q+ 2 p(E)/o^Q+(E) given -in Table 5 produce theoretical I 7319 /I 4278 ratios 

at the error extremes of the data. 

In the ‘high quenching, OF case, a range of .05 to .24 was deemed acceptable 
according to the earlier discussion surrounding Figure 45. It was noted that even the .20 
branching ratio curve (G^® of Figure 45) did not totally fit the data curve. If the charge 
exchange production mechanism completely dominated over electron impact, the 
theoretical curve would bend closer to the measured curve at the lower altitudes. 

The ‘high quenching, N 2 ’ theory required a large efficiency in a^.,(E)/a .(E) 

S 

to even skirt the error bars as discussed in relation to Figure 46. Even with all production 
due to the N'^(‘S) mechanism, the criteria of falling in the measurement error is not met. 
The author feels this mechanism is least likely. 

In summarizing the investigation then, some conclusions can be drawn. First, in 
all theories, the upper limit acceptable for the branching ratio 0^+2 p(E)/o^^(E) is .030. 

This is significantly lower than the values proposed by Seaton (1959) and Kazaks et. al. 
(1972) of .2 and .25 respectively. It has further been postulated that this ratio could be 
zero if large quenching does exist. Secondly, an upper limit of .24 has been established 

for a . , „(E)/o .(E) considering all three theories, with a lower limit of zero. Once the 
N"''* S 

quenching rates have been truly established, the bounds on the branching ratios can be 
narrowed. 

A puzzle remains. Why should the efficiency for producing 0"^(^P) by electron 
impact be at least a factor of 7 lower than the .2 value deduced by Seaton (1959)? This 
is especially puzzling since Walker et. al. (1975) were successful in describing the dayglow 
satellite measurement with a branching ratio comparable to Seaton’s. Lassetre (1969) 
suggests that differences in cross sections do exist for photon interaction versus electron 
impact in several molecular systems. It is difficult to understand why the difference here 
described is so large. 

8.8 A RE-EXAMINATION OF EARLY MEASUREMENTS 

The previous section contined discussion of several theoretical possibilities which 
describe the I 7319 /I 4278 measurements. It has been shown in this discussion that 
high quenching of 0'*‘(^P) as proposed by the Walker et. al. (1975) measurement can 
describe 0'*'(*P) continuity balance in aurora, provided the electron production efficiency 
for producing 0‘*'(*P) is reduced from the .20 value of Seaton (1959), and an additional 
low altitude source mechanism, such as N'*’(‘S) charge exchange, is introduced. Yet 
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unexplained is why Gerard (1971) concluded that the low altitude source was not present 
and that reduced electron impact was not necessary to explain early auroral observations. 
Many of the historical observations of these emissions have been accomplished from 
ground observatories. Historically, spectrometers from these ground stations viewed 
distant arcs and discerned altitudes of maximum emission of the respective emissions. To 
re-examine these ground based measurements it would be useful to plot a volume 
emission rate profile of the X7319A and X4278A emissions for a typical aurora as 
predicted by these theories. Since these early measurements viewed the aurora nearly 
perpendicular to the magnetic field rather than the magnetic zenith, it was assumed that 
the ratio measurements reflect volume emission ratios. 

Plotted in Figure 47 are two curves which approximate the relative volume 
emission — altitude relationship of the two emitting species, 0'^(*P) and The 

t}(4278A) profile was assumed as being the N} -proportional type discussed in Appendix 
F. That is, the i?(4278A) emission was defined proportional to above 120 km, and zero 
below. Equation (80) of Chapter VI was used to generate the t 7(7319A) profile with the 
most probable values of the ‘Low Quenching’ theory of Table 5 as input parameters. The 
‘High Quenching, OF theory with most probable values of Table 5 predicts a nearly 
identical profile for t?( 7319A). The important point to be stressed is the volume emission 
peak for ij(7319A) lies very low (~ 135 km in this case). A volume emission peak of ~ 
220 km is predicted by Gerard (1971) using the Seaton (1959) value of .2 for 
ao+ 2 p(E)/aQ^(E) and large Nj quenching with no additional 0‘‘‘(*P) source other than 

electron impact on 01. 

Plotted in Figure 48 are volume emission ratios — altitude relationships for the 
‘Old Theory’ (yjq^ = 1.10*‘®, 7 g = 3.10"®) verified by Gerard (1971) and a new theory. 

The new theory plotted is the low quenching model although, as discussed earlier, the 
‘High Quenching, OI’ (Most probable, Table 5) theory now predicts a near identical 
profile. The curves were generated using the reciprocal of Equation (78) with the input 
parameters reflected by the previously discussed quenching and u ^ (E)/a .(E) 
efficiencies. 

To culminate this study, reference is made to the early measurements and 
speculation as to why the high efficiency and quenching rate were verified by Gerard 
(1971) and others, especially in light of the Wallace (1960) data. 

Jones (1971) mentioned that the Oil Auroral lines were at least three times the 
intensity of the NjlPGs^s band for the Jones’ (1960) spectra of the February 10/11, 
1958, aurora. A relative intensity of 2.8 at X7384A was reported for the spectra 
measured with an instrument containing a 5. 2 A instrument function. 

A 5.2A triangular slit function peaked at X7384A was passed through the Nj IPG 
band model, previously described, with a rotational temperature of 700°K. This resulted 
in a determination of the percentage of the band intensity measured at the quoted 
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ALTITUDE (km) 



RELATIVE VOLUME EMISSION RATE (ARBITARY UNITS) 


Figure 47. Altitude versus relative volume emission rates of t}( 7319A) 
and t 7(4278A) assumed as a typical aurora. The t}(4278) is 
proportional to nlNj] to 120 km and zero below. 
t](7319A) was produced by a new model 
(Low Quenching, Table S). 
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,NEW THEORY. (LOW QUENCHING! 


■OLD THEORY" 


7} (3914AI/,, (7319M 

Figure 48. Altitude versus i}(3914A)/t7(7319A) for the “Old Theory 
(Oo+ 2 p(E)/Oo+(E) = .2 and = 10"*®) and a new model (Low 

Quenching, Table S). This figure is meant to aid in interpreting 
early ground based spectroscopic measurements. 





wavelength and amounts to ~ 5% of the total N2 IPGs ,3 band intensity. This would infer 
that the relative intensity of IPGs ,3 was 56 resulting in an ij(7319)/i7(7387) of .3. 
Jones (1971) suggested the doublet was three times the intensity of NjlPGs,3 for the 
same data but the above analysis disputes this. It was pointed out in the rocket data 
analysis earlier that IlNEGi,o/Ili*Gs,3 was ~ 1 at all altitudes, which would infer that 
»?(1 NEGi,o)/'* 7(1 PGs, 3) should also be one. This would mean that the interpretation of 
the Wallace (1960) data would correspond to a T?(7319)/r?(4278) of .3, or 
»7(3914 )/i 7(7319) of ~ 11. From examination of Figure 48, the ‘Old Theory’ would 
suggest a corresponding altitude of ~ 150 km. For the ‘New Theory’ of the same figure, 
the altitude indicated would be near 220 km. As pointed out earlier, Jones (1971) 
classified this type of aurora as having a lower border of ~ 250 km. 

The main argument given by Gerard (1971) for observational verification of his 
volume emission rate theory of Oil Auroral emissions was Omholt’s ratio of .1 for 
T?(7319-30)/i7(5577) as measured in the upper portions of a spectrograph taken in 1957. 
As reviewed by Jones (1971), the green line to NjlNEGo,© ratios are believed to be in 
the neighborhood of two. This inference leads to tj(3914)/tj( 73 19-30) of 5 of 
77(39 14)/t?( 73 19) ~ 8.3. Again, examination of Figure 48 suggests ~ 170 km for the ‘Old 
Theory’, and the ‘New Theory’ would correspond to ~ 235 km. It is difficult to 
speculate how high Omholt’s ‘upper’ portions of aurora were. Spectra exposed to aurora 
for long time periods is subject to exposure to many different types of auroral displays. 
The fact that the lower parts of the Omholt spectra showed ratios and intensities similar 
to several other spectra, where the doublet was not noted, could imply that these types 
of aurora took place, as well as higher altitude “red” type events during the exposure 
period of the particular plate. 


CHAPTER IX 

SUMMARY AND RECOMMENDATIONS FOR FUTURE WORK 


The previous eight chapters were dedicated to the description of research 
performed in two disciplines. Background and early investigations were siunmarized for 
both the spatial scanning photometer and auroral 0'*'(*P-*D) emissions in the first 
chapter. 

9.1 SUMMARY OF THE SECTOR SPECTROPHOTOMETER DEVELOPMENT 

Chapters II through V discussed the evolution of a unique spatial scanning 
photometer system. 

The second chapter gave the reader an understanding of the theoretical principles 
of the spatial scanning scheme. This was accomplished with a simplified half-width 
treatment rather than with analytic transmission functions. This chapter concluded that 
theoretically, at any rate, a rotating spiral slit could be designed into a photometer system 
so that the instrument would have many desirable characteristics over other photometers. 
It was deduced that a scanning aperture could be built which had a near constant 
throughput and bandpass for all slit positions. It was also recognized that the wavelength 
scan rate could be made proportional to the rotation rate of the spiral slit and that a 
continuous wavelength region was accessible for a continuous spiral. In essence, the 
evolution of an instrument having scanning characteristics of grating type spectrometers 
with the throughput advantages of flat plate interference elements was possible. 

A study of the interference elements both theoretically (Appendix A) and in the 
laboratory (Chapter II) was required. The purpose of this study was to determine to what 
angle filters could be pushed from normal (i.e., what wavelengths), to determine filter 
cone angle effects, and to explore general filter characteristics. The two effects limiting 
filters in this application were the polarizing effect (more prevalent in low index filters) 
and a slight low angle broadening of undetermined origin. In Chapter II, the cone angle 
effect of the filters did impact the characteristics of slit design in striving to achieve a 
constant instrumental throughput and broadening. All filters were found to differ. The 
limitations of these elements as they apply to the sector spectrophotometer were covered 
in this chapter and the noted Appendix. 

A single spiral photometer was built and tested in the laboratory to validate the 
theory described earlier. The development and laboratory analysis were discussed in 
Chapter IV. It worked! 

A more complex device which allowed for more light gathering power at a given 
resolution was described in Chapter V. This system was designed for ground based 
measurements of auroral 0’’'(^P-*D) emission. This instrument system included a 
multi-spiral aperture. Multi-spiral slits can provide two to three times the instrumental 
throughput over single spiral devices for a given bandpass resolution. 


Two disadvantages of the system include the large detector for image capture, 
which results in more detector noise than would be in, e.g., a tilting photometer system 
with a fixed field stop. The second disadvantage is that in spatially shifting the aperture 
and changing the look angle through the filter, the source is also spatially scanned (i.e., a 
wavelength scan describes an arc in the sky). A large lens imaged at the interference 
element corrects for this, but does cause some throughput loss. 

9.2 SUMMARY OF THE O+l^P-^D) IN ELECTRON AURORA 

The theory of 0‘*'(*P-*D) as related to Nj’lNEG emissions was the subject of 
Chapter IV. The analysis of ground based data and a rocket payload pertinent to 
0+(2p-2D) analysis were the subjects of Chapters VII and VIII respectively. 

In Chapter VI, pertinent cross sections and their ratios were analyzed. An 
expression of the ratio of the volume emission rates of 0'*'(^P) to N^INEG emissions 
expected in an electron bombarded atmosphere was developed. The ratio was found to be 
functionally dependent on collisional deactivation rates and populations of Nj and e in 
addition to the cross section ratios and transition probabilities. The source and values of 
the ratio dependent parameters were discussed from a historical standpoint. 

Ground based measurements of I 7319 /I 4278 ratio irr the magnetic zenith were 
presented and compared with previous theories. It was found that the measurements 
produced lower ratios for aurora than early theories predicted. 

Rocket payload data from an aurora were analyzed and added to the data base of 
I7 3 1 9 /I4 2 7 8 irr Chapter VII. The rocket analysis measurements were found to be in good 
agreement with the ground based measurements of the previous chapter and previous 
theories did not describe the measurement. New theories were presented that did describe 
the data which included a ‘low, quenching’ model and two ‘large quenching’ models. The 
‘large quenching’ models required a large, low altitude source for 0"''(^P). Charge 
exchange of 0(^P) with N‘*'(‘S) was a suggested source for the ‘large quenching’ models. 
The ‘large quenching’ models were required to coincide with the Walker et. al. (1975) 
dayglo w results and it was found that the ‘large quenching, 01’ model represented the 
auroral data better than the ‘large quenching, N 2 ’ model. In all the new theories, a large 
reduction in the direct electron production efficiency of 0"*'(^P) (i.e., a . j (E)/a .(E)) 
from the Seaton (1959) value of .2 was required. 

9.3 AREAS FOR FUTURE RESEARCH 

It would be of interest to make the auroral rocket measurement to a higher 
altitude and verify the findings of NASA 4.329. I believe that an N 2 IPG band emission 
would suffice for the nitrogen component to the ratios of interest. A narrower bandpass 
filter (something ~ 7-1 OA, peaked at X7319A) would reduce the background and give a 
better signal to noise ratio. Actually, a tilting filter might be ideal to give adjacent 
background information for an improved measurement in light of the ©"^(^P-^D) 
faintness. (The sector spectrophotometer could also perform the measurement, but the 
tilting filter would suffice since no OH would be present.) 
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Other measurements from the ground would be ideal, especially in the ‘soft’ 
auroral zone poleward of the oval. Bather and Mende (1971) have done morphology of 
6300A and 4278A in this region. Examination of their data and Figure 44 suggests 
intensities of 10 rayleighs of 17319 should be present according to the results found in 
this analysis. 

It would furthermore be interesting to look for the emission in SAR arcs. Even 
though the emission is less than it might be in one of the molecular b^d systems, the 
line nature of the emission might make it more attractive from an analysis standpoint. 

Further clarification of the production efficiency of the metastables over the 
ground state by electron impact may be feasible in the laboratory. Ideally, these cross 
sections should be known out to ~ lOkev or even larger. Simultaneous measurements of 
the 0 '*‘(^D-‘*S) and 0"'‘(*P-*D) doublets would be of interest for clarity as to the 
branching or metastable partitioning probability per ionization. 

The instrumentation studied was designed for a particular purpose and each 
wavelength region has its own peculiarities. Primarily, however, I believe that the greatest 
improvement required lies in the filter. The availability of narrow band filters, with good 
off axis transmission characteristics (minimum polarization etc.), is the main component 
of the system studied. The value of optimizing the spatial technique is dependent on the 
quality of this first element of the system. The high throughput advantage, however, 
vastly improves signal to noise characteristics offered by conventional grating 
spectrometers, which is a trait of the Fabry-Perot type interference generators. Where 
short wavelength scans are satisfactory and a high throughput is required, this system 
should be examined. 

The author ranks the ‘large quenching, OF model as being the favored theory of 
the three presented. The Walker (1975) data strongly supports large quenching of 0'*‘(^P) 
for dayglow continuity and certainly the auroral atmosphere differs little from the 
daytime atmosphere in terms of populations of quenching species. With large quenching, 
there must be a low altitude source of 0'*’(*P) in aurora in addition to direct electron 
bombardment on 01. The N'^(*S) mechanism is a very plausible one. The ‘high 
quenching, N 2 ’ model is less favored due to the fact that the Nj scale height leads to a 
marginally acceptable synthesis of the NASA 4.329 measurement. It is possible that 
temperature effects in the quenching rate (by Nj or OI) may be large. In addition to the 
random thermal rate of collisions, electric fields can produce large relative non-neutral 
motions, particularly at the higher altitudes. The higher altitudes could have an 
accelerated deactivation rate over the low altitudes. This effect would tend to nullify the 
arguments made on the OI versus Nj scale height effect on the altitude shape of the 
emission rates. These effects are not known in these measurements. 
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APPENDIX A 


FILTER THEORY 


Narrow bandpass filter (5X|| < lOA) technology is continuously being refined. A 

review of fundamental filter theory will be explored in this appendix to aid in 
understanding its utilization in the instrument design. The reader should realize that 
building filters is an art left to thin film specialists who put theory into practice. I will 
attempt to give some insight into theory, such that, these effects might be understood 
and compensated or anticipated by a user of the instrument investigated. The art of 
communicating to the manufacturer the filter requirements for a particular task requires 
some knowledge of filter theory and how and what particular optical properties are 
inherent in the dielectric filter. In addition to cone angle and dispersion effects that most 
users are aware of, off axis transmissions introduce polarizing properties that must be 
considered. 

It is the purpose of this appendix to summarize some filter theory so that some 
insight into optimizing specifications or trading off lesser undesirables of a filter for 
others, might aid the reader. The work of Lissberger and Wilcox (1959), Lissberger 
(1959), and Pidgeon and Smith (1964) offer two approaches to the theory. The approach 
of Pidgeon and Smith (1964) offers a rigorous matrix layer inversion scheme and 
transmission equations that give a laboratory specialist a tool in theoretically 
programming and building a filter. The earlier work of Lissberger and Wilcox (1959) and 
Lissberger (1959) follows a less rigorous approach, but one which yields nearly the same 
results. I will follow the Lissberger and Wilcox approach. 

The transmission of a light beam through two parallel reflective surfaces (such as 
a Fabry-Perot) is given by the well known Airy formula (Bom and Wolf, 1970) 


I(t) = 


I(i) 

1 + F sin^ 6/2 


(A-1) 


where I (t), I (i) = transmitted and initial beam intensity 


4R 

(1 - R)* 


(A-2) 


where R is reflectivity of the surfaces 
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5 = 2 7rm = 27r-?^ 
A 


(A-3) 


ju — spacer index of refraction 
t - spacer thickness 
X — wavelength . 

The half-width (full-width at half height) of a monochromatic beam can be expressed in 
terms of the above expression at half-width 

6 = 2 mw. ± ^ , (A-4) 


where (A-1) becomes 


1 = 1 

2 1 + F sin* e/4 


(A-5) 


or, for e small 


e 


4 


(A-6) 


The separation of neighboring orders (finesse) can therefore be written as 
F = ^ . 


(A-7) 


All dielectric filters are made with multiple layers of quarterwave layers of 
alternating high and low index for the reflective stack, and a multiple half-wave layer for 
the spacer. Theories have been derived for filter stacks similar to the simple two-layer 
Fabry-Perot reviewed above. Lissberger (1959), Lissberger and Wilcox (1959), and 
Pidgeon and Smith (1964) have made theoretical and experimental studies of the 
all-dielectric filter. A sketch of the theory involved and, in particular, examination of the 
major filter characteristics as they apply to the sector spectrophotometer will follow. 
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Figure 1-A depicts a schematic of a simple all-dielectric filter borrowed from 
Lissberger and Wilcox (1959). The dielectrics used in stacks are usually cryolite (hq - 

1 .34) and zinc sulfide (jijt * 2.34) which in turn are built on a glass substrate. Lissberger 

(1959) and Lissberger and Wilcox (1959) developed a general theory for the simple type 
niter shown in temis of the Fresnel coefficients of the boundaries and optical thicknesses 
of the layers. Pidgeon and Smith (1964) did a similar study with similar results. The 
subsequent review will follow the Lissberger (1959) and Lissberger and Wilcox (1959) 
and for brevity, this work will be abbreivated as LW and the Pidgeon and Smith (1964) 
as PS. 


The multilayer stacks (Figure 1-A), separated by a spacer layer of thickness t, can 
be represented by equivalent “single” layers separated by the same spacing. The 
equivalence relationship is ^own in the right half of Figure 1-A where the equivalent 

Fresnel amplitude (R^*^) of the air and glass stacks are related to initial amplitudes 

^l(OA,OG) and phase change within a particular stack (S'^*^). For details of the 
derivation, reference can be made to LW papers. The following are the results of the LW 
work in abbreviated form. 

The Airy formula takes the form: 


s ^max 
^ 1 + F sin» x/2 


where 


X = 6 + 6^ + fiA 

6 - previously defined 

5^,5^ — phase change in the glass and air stacks, respectively 
Tmax ~ transmission. 


(A-8) 


(A-9) 


By making a small angle assumption (cos d 1-6^ 12), excluding polarizing effects, and 
accounting for dispersion, the angular dependence of Tj^^x versus X becomes 


Xo 


|2(1-Q) 



(A-10) 
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Figure 1-A. Schematic of a dielectric filter as referenced in Lissbeiger and Wilcox (1959). 




where 


AX = X - Xo 

Xo — wavelength of for 0 = 0 (where d is external incident angle of 
light to filter) 


where for a high index (Zjj5) spacer 


P = -- 


n + 


n + 


lc_. 

Mz-Mc 


and 


Q. = 


n + 


Me 


Mc-Mz 




z 


5^ ^ MqMz ^ 
6X Mz-Mc Me 


5mc 

6X 


or for a low index (C) spacer 


p ^ I n + (p^Mz(Mz-Mc)) 
Mg I ^ ” Mc^^ 


Qr = 


McM; 


c'*z 


■6-) 


m| -mJ 


^2 5X 



(A-11) 


(A-12) 


(A-13) 


(A-14) 
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The half-width (for parallel - normally incident light) for a first order filter, as 
derived by PS, may be written as 


6 \ 


II 



(A-15) 


where 


Fk 


MkMr 


2x+l 


l(Mr-l) 


(A-16) 


Mr = Mh/M2 
Mjj — High Index 
Ml — Low Index 

X — Number of layers of reflective 1/4 wave layers 
M{^ — Index of spacer . 


The above relationship suggests that the narrow bandpass is achieved by many layers 
(high X) and to a certain extent, with high index spacers. 

The cone angle has a broadening effect on the transmission characteristics of a 
filter. According to LW, the bandpass relationship is as follows: 


6Xp - 




where 


6X^ - 



(A-18) 


(A-19) 
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P,Q,Xo — previously defined 

\p — cone 1 /2 angle (in radians) . 


At visible wavelengths, a 3° cone 1/2 angle corresponds to ~ 2A broadening. Therefore, 
for narrow bandpass applications (5X|j < lOA at visible wavelengths) care should be taken 

not to design an instument such that the cone angle effect alters the bandpass 
significantly. On the other hand, one normally wants as large a cone angle as possible 
where light gathering power is important, such as in nightglow. 

Dispersion has been discussed by LW. For ZjjS type filters, the dispersion is small 

and normally not of concern for \ < 4500A. For blue wavelengths, it becomes more 
significant. The effect is to make the tilting slope (P/l-Q) smaller than it would be if 
dispersion were absent. Dispersion of Zj^S is discussed in detail in Hall and Ferguson 

(1955). Dispersion due to cryolite is negligible in the visible wavelength region (LW). 

The effect of most concern, and observed in many filters by this author, was the 
polarizing effect with tilt. The effect was complex and was simply treated such that some 
indication of the filter characteristics primarily responsible for the phenomena might be 
identified. The simplified approach was only meant to give tendencies and not to be 
taken as explicitly defining the phenomena. As will be observed in the text of the 
laboratory analysis, some filters showed very large polarizing characteristics at very small 
angles (~ 2° tilt), much smaller angles than this theory would predict. Others polarized at 
larger angles than would be predicted. The following theory does suggest that the high 
index filters tend to polarize at larger tilt angles than the low index filters, and this 
general tendency was observed in the laboratory analysis discussed in the text. 

The approach was to redefine Equation (A- 10) in terms of AX|| and AXj^ where 
Pj, Qj^, P||, and Q|| (Equations (A-11, 14)) were redefined in terms of m"*" = M cos 0 and 

jull = n/cos 0 where 0 is the angle of refraction in a given layer. This was related to the 
external angle of incidence by Snell’s law of refraction. The separation of the two 
components can then be represented by 


6Xp — AXj_ - ^^11 


(A-20) 


Figure 2-A shows a plot of some theoretical separations (5Xp) based on this simplified 

approach. To minimize the polarizing tilting effect, a higher spacer order is indicated, as 
well as, a higher index spacer. The spacer index is the most marked effect but recall that 
the wavelength shift for a given tilt is almost twice as great for a low index filter. 

The transmission of one component versus another with tilt has not been treated 
here. However, the total bandpass integral did not change for several filters tilted to 0 > 
30°. This study should be carried further with a more refined theory, to account for the 
observed separations. A thorough knowledge of the particular filter stack and of the 
individual components with tilt for many different spacer indexes, orders, and stack 
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layers would be necessary to understand this effect. The study in Baumeister (1961), if 
extended to include the spacer stack, would be an ideal tool for analyzing particular 
requirements of this sort. The use of polarizers (although a loss of 60% in throughput 
results) is recommended in instrument design, where large tilt angles or large wavelength 
scans are desired. 

For other sources of interest, the author suggests referring to Baumeister and 
Jenkins (1957), Father and Reasoner (1969), MacLeod (1969), and Thelen (1966). 
Temperature effects are discussed by Baker and Yen (1967). 
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APPENDIX B 

FORTRAN IV PROGRAM USED IN GENERATING THE DATA ARRAY 
FOR THE APERTURE SPIRAL AND A MASKING SPIRAL 

INPUT: XNA = Breadth finesse 

XNC - Curvature finesse 

F = Focal length of objective (cm) 

RM = R Maximum (cm) 

DPHI = Angular increments desired in array (degrees) 

AMAX = \ • Xo , at R maximum (angstroms) 


OUTPUT: PHII = Angle of aperture point from r = 0 reference (degrees) 

RDRP, RDRM = Radial distance bounding aperture at given angle (cm) 

PHSM 1 , RSMM = Angle and radius of inner masking spiral boundary (degrees, cm) 
PHSPl , RSPP = Angle and radius of outer masking spiral boundary (degrees, cm) 




L 


1 

2 

3 

4 

5 

6 

7 

8 
'i 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
2d 
29 
3U 

31 

32 
:>3 

34 

35 


UNITS IN CM 

READ! 5tl)XNAtXNCiF»RMtN(OPHI «AMAX 

FORMAT (lXi-3.0tlXF4.lt 1XF5.2* IXF 5.2t IXI 2t lXF5.0t lXFo.2) 

KA0=57. 29578 

XI^=. 39370 

CM=2. 54000 

KM=RM*CM 

TM=ATAN(RM/F) 

TR1=TM*RAD 

DRM=F4TM*((l.0/CCS(TM))**2)/XNC 

PI=3. 14159 

A=RM*RM#PI 

UPHI =0PHI/RAU 

PHI = 4.0/RA0 

PHM=2.0«P1*N 


DELPHI=PHM/XNC 

PHMM=PHM+90./RAD 

GMtG=OhLPHl*TM*TM/XNA 

A01=A*0MEG 

ALPH=AMAX/(TM*TM) 

NN=1 

10 K=F*TANITM4(PH1/PHM)**.5I 

DR=URM*( ((R*R+F*F)**l.5)/( ( KM*RM+F*F ) 1. 5) )*(RM/R) 

RDRP=R+UK/2.0 

RDRM=R-OR/2.0 

T=ATAMR/F) 

A02=PI*RM*RM*UR*DELPHI*C0S(T )*CUS ( T I *S I N( T )/F 

ANA=T»CUS(T)*COS(T)*JR/(F*TM*TMI 

ANA=1 .0/ANA 

C THE FOLLOWING PERTAINS TU THE SHUTTER 
PFS = PHI/3.0 
PHSME=PhI-UELPHI/2.U 


PHSPE=PHI+ULLPHI /2.0 


IF(PHSME)63t63,64 
63 PHSME=0.0 
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36 


RSM=0.0 

37 


ksmm=o.o 

36 


RNh=0.0 

39 


GO TO 95 

40 

64 

B=F»TAN( TM*(PHSME/PHM }**.5) 

41 


RSM=F*TAN( TM« (PHSMt/PHMlt'^.B 1+ 1 ORM/2 .0 ) ♦( U Q*B<-F*F ) 1 . 5» / ( 1 

42 

1F*F 5) )*(RM/B) 

43 


R5MR=RSM+.1000 

44 


C^F^TANITM* ( (PHSME+2.0*PI+DELPHI)/PHM)«*.6) 

46 


RNE=C-( DRM/2.0)<“( ( ( C*C+F #F ) ♦♦! . 5) /( ( RM*RM+F ♦F J ♦♦ 1. 5) ) ♦ ( RM/C » 

4b 

95 

CCNTINUE 

47 


0=F*TAN( TM*(PHSPE/PriM)**.5) 

48 


RSP=F*TAN(TM*(PHSPE/PHM)**.6)-(0RM/2.0)*( ( (l)*D+F»F )**1.5) /( ( RM*RM+ 

49 


IF*F)**1. 5) }*(RM/D) 

60 


RSPP=RSP-.1000 

61 


PHSM=PHS-DELPHI/2.0 

62 


PHSP = PHS+OELPHI /2.0 

53 


CHEC=PHSPE-2.0*PI-DELPHI 

54 


IF(CHEC)98»96,99 

65 

98 

RNLT=0.0 

56 


GO TO 100 

5 7 

99 

E=F*TAN( TM*l (PHSPE-2.0*P I-OELPHl )/PHM)*»,5) 

56 


RNLT=E+( DRM/2.0i*( ( ( E *£+F*F ) 1 . 5) / ( ( RM*RM+F*F ) ♦♦! . 5 ) »*(RM/EJ 

59 

100 

CCNTINUE 

60 


PHS1 = PHS<‘RA0 

61 


PHSMl=PHSM*RAD 

62 


PHSP1=PHSP*RAU 

63 


PHI l = PHl*RAU 

64 


iF(PKlI.LE.360. )G0 TO 27 

65 


PHI I=PHII-(NN-1 1*360. 

66 

27 

CCNTINUE 

o7 


PIP=Phl/NN 

68 


PI2 = PI*2 .0 

69 


XRI T£ (6, 25) PHI I tRORP»RDRM,PHSMl ,RSMM»PHSP1 , RSPP 
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7Q 

25 

FCJRMAT(F6.2,F7.5f F7.5,F6.2,F7. 5,F6.2,F7.5) 

71 


IF(PI2.G£.PIPIG0 TO 29 

72 


NN=NN+1 

73 

29 

CONTINUE 

74 


IF(PHI-PHMM)ll,ll,12 

75 

11 

PHI = PHH-UPHI 

76 


GO TO 10 

11 

12 

WRITE (6, 17) TM1,URM,A01»ALPH 

18 

17 

FDKMAT(/1UXF5.2,2XF6.5,2XE10.3,2XE10.3) 

79 


WRITE (6,1) XNA,XNC,F,RP,N,CPHI,AHAX 

80 


END 

81 


END 


END OF FILE 



APPENDIX C 


SPONTANEOUS EMISSIONS FROM O^i^P) 


The 0'*'(^P) “nebulium lines” (Herzberg, 1944) have been studied and reported by 
astronomers due to their noted dominance in spectra made of cosmic nebulae. Studies of 
the 0’*‘(^P) and the forbidden emissions originating from the state provided astronomers 
with a radiative signature of the environment in which such emissions were originating. 

Seaton and Osterbrock (1957) calculated the effective transition probabilities, 
emission wavelengths, and line strength factors for spontaneous emission from the 0‘*‘(^P) 
and ©■'■(^D) states. Figure C-1 portrays the fortrat diagram for those magnetic dipole and 
electronic quadrapole transitions originating from the 0'*'(^P) state. The three equivalent 
electrons of the outer shell can be in a '*S, or ^P term configuration with total 
angular momentum quantum numbers of 3/2 for the S term, 3/2 and 5/2 for the D term 
and 1/2 and 3/2 for the P term. The spin coupling results in two ^P states (^Pj/j and 
^Pg/j), where the relative populations of the two states are given by their statistical 
wei^ts. Spontaneous cascade from the ^P states can occur through either a magnetic 
dipole or electronic quadrapole transition to the levels or the ‘*S level. 

Table 1-C depicts the atomic transition parameters for spontaneous emission from 
the ^Pi /2 and ^Pa /2 states of Oil. Note the transition probabilities represent the total 
transition probability including both the electronic quadrapole and magnetic dipole 
transitions. 

The total radiative loss rate of 0'*‘(^P) is given by the sum of the loss occurring 
from the ^Pj /2 state and the state, i.e.. 


7?(0+(^P)r = 77(0+^P,/2)r + 77(0+2 Pa /2>R 


(C-1) 


The loss rate of the ^Pi/z state becomes 


T?(0+2P,/2)r = n[0+2p,/2] • A2P^/^ (C-2) 

and for the ^Pz/z state 

7?(0+2P3/2)r = n[0+2P3/2] ‘ A 2 . (C-3) 
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DIAGRAM FOR O 0 METASTABLE, 2p TRANSITIONS 

(AFTER SEATON AND OSTERBROCK, 1957) 

STATISTICAL 

WEIGHT J TERM 



Figure 1-C. Diagram for Oil Metastable, transitions 
(after Seaton and Osterbrock, 1957). 
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TABLE 1-C. ATOMIC PARAMETERS FOR 0+(*P) TRANSITIONS 
(SEATON AND OSTERBROCK, 1957) 


Transition 
n -»• n' 

Wavelength 

(A) 

Ax* 

(sec‘* ) 

w** 

*Pi/j 

"Da /a 

7329.9 

.100 

2 


"Ds/a 

7318.6 

.061 

2 

*Pl/2 


2470.3 

i024 

2 

"P3/2 

"Ds/a 

7330.7 

.061 

4 

"Pa/a 

"Da/a 

7319.4 

.115 

4 

"Ps/a 

^Sj/a 

2470.4 

.060 

4 


^Includes summed electric quadrapole and magnetic dipole transition 
probabilities. 

♦♦Statistical weight, 2J + 1, of the upper state (n). 




The statistical weights of the respective levels relate the populations of the substates to 
the total 0"^(*P) population, i.e., 

n[0+*Pi/a] = I • n[0+*P] ' (C-4) 

6 

and 

ntO+^Pa/a] = | • n[0+*P] . (C-5) 

o 

Substituting (C-4) and (C-5) into (C-2) and (C-3) respectively, and substituting the 
respective transition probabilities from Table 1-C, the respective *P loss rates become 


T?(0+2p,/,)j^ = .062n[O+»P] 


(C-6) 


and 


^?(0+*P3/2)r = .157n[0+"P] 


(C-7) 


Further substitution of equations (C-6) and (C-7) into (C-1) gives the total radiative loss 
rate of ot(^P) as 


77(0+2p)j^ = ,219n[0+*P] 


(C-8) 


Considering the rocket payload measurements discussed in the text, it is desirable 
to represent the population of 0'*'(^P) in terms of the fraction of the 0‘*'(^P) radiating to 
the 0'*‘(^D) state. This total emission rate was measured on NASA 4.329 as discussed in 
the text. Equation (C-8) can be written as 


t?(7319-30A) + t?( 2470A) = .219n[0+*P] 


(C-9) 



where tj( 7319-30A) represents the radiative loss rate of *P to *D states and t7(2470A) the 
radiative loss rate of ’P to ^S. By substituting the emission probability and statistical 
weights for t7(2470A), 


t?( 7319-30A) = .219n[0+*P] - | • .024n[O+*P] - | • .060n[O+»PJ ,(C-10) 

o ft • 


or 


77(73 19-30A) = .171 n[0+*P] 


(C-11) 


The ground based measurements are directed at only the 7319A doublet. By 
expanding (C-ll) into the respective components, and representing the 7330A emission 
by the statistical weight and transition probability, (C-ll) becomes 

t 7(7319A) = .171 n[0+»P] - | • .100n[O+*P] - ~ • .061 nlO+»P] (C-12) 


or 


t?( 7319A) = .098n[O+*P] 


(C-13) 
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APPENDIX D 


PRODUCTION OF 0+(*P) BY ELECTRON IMPACT AND ITS 
RELATIONSHIP TO EXCITED N| EMISSIONS 


The volume production rate for an impact source such as energetic electrons is 
determined by the product equation 


T?(x')p = / F(E) ax'(E) n[x] dE 
E 


(D-1) 


where F(E) — energy dependent flux of streaming particles through the background gas, 
x; ajj'(E) — energy dependent cross section for producing an x' state from an impact on 

species x; n[x] - volume density of species x; and i?(x')p — production rate of species 

x'. The flux distribution (energetic and pitch angle distribution) of precipitating auroral 
electrons is difficult to measure but has recently been accomplished by several 
investigators for the low ev to tens of Kev range. The intent here was to relate the 
production rate to spectroscopic signatures since electron distribution measurements are 
not always available, particularly when making ground based measurements. It has been 
previously noted and referenced that the ionization rate of Nj is related to the 
NflNEG emission. The approach taken was to relate the ionization rate of 01 to the 
ionization rate of N 2 , which in turn is related to the NjlNEG emission rate. 

Equation (D-1) can be written for Nf and O'*' to read as, 


»?(N?)p = / F(E)a .(E)n[N2] dE 

c ^2 


(D-2) 


and 


Tj(0+)p = / F(E) a +(E) n[OI] dE 
E 


(D-3) 


Borst and Zipf (1970) have found in laboratory analysis that for E > 30ev electrons, 1 % 
of all generated emit the NtlNEGo.o radiation. The 7% probability is nearly 
independent of energy and permits the production rate of to be written as: 
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7j(Nt)p = — T7(NtlNEGo,o) 


(D-4) 


where is the cross section for producing state in the zeroth vibrational 

level. The identity now exists between Equations (D-2) and (D-4) that 


/ F(E)o„+(E)n[Nj] dE - t _. ij(3914A) . 

E “Nf® 


(D-5) 


Equation (D-3) can be rewritten to look very similar to Equation (D-5) by removing the 
non-energy dependent n[OI] to the left side of the equation and multiplying by n[Nj], 
i.e.. 


T?(0+)p = / F(E) a^+(E) n[Na ] dE 


(D-6) 


The only difference between the integral portions of Equations (D-5) and (D-6) lies in 
the cross section. 

The energy dependence of the total ionization cross sections of Nj and O are 
summarized by Kieffer and Dunn (1966). More recently, Kazaks et. al. (1972) have 
added theoretical deductions for the oxygen cross section. The energy dependence of the 
cross sections of the two constituents are very similar. The major discrepancy in the cross 
sections of concern lie in 01 from 20 to ~ 300ev. Broksenberg (1961) and Kazaks et. al. 
(1972) deduce a cross section of 30% higher at lOOev than Fite and Brackman (1959) 
and Rothe et. al. (1962). ForE > 300ev, little disagreement exists in the OI cross section. 
The assumption was made after careful examination of these cross sections (Figure 1-D) 
that 


V'® 

— = Constant (D-7) 

V® 
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where the constant was deduced as being .5 by this author. In other words, it was 
assumed that at all energies, the OI ionization cross section was one half that of the Nj 
cross section. Note in Figure 1-D that for E < 50ev, the ratio approaches one. It is the 
author’s deduction that the .5 ratio gives a good representation of the overall cross 
section difference for the bulk of the auroral primary electron energy range. 

Equations (D-6) and (D-7) can now be combined to give: 


= /F(E)%+(E)nlNjl dE 
a^+(E) n[OI] g NJ 


(D-8) 


An identity can now be written by eliminating the integral of equation (D-8) and (D-5), 
resulting in a relationship between the production rate of O'*’ ions to the NjlNEGo,© 
volume emission rate, i.e.. 


a (E) 

7j(0'*')p = 






• tj(3914A) 


n[m] 

n[Nj] 


(D-9) 


The next step is to determine the rate of 0'*'(^P) production specifically. The 
0'*'(^P) production is related to the O'*" production rate by the relationship 


t7(0'^"P)p = 


V® 




(D-IO) 


By substituting Equation (D-9) into (D-10), the 0'*’(*P) production rate becomes 


77(0'*'2P) 


<^0+2 p(^) t?(3914A) 


(D-11) 


The first three cross section ratios could be reduced to simplify Equation (D-11), but 
previous investigations have presented values for these particular ratios. The first ratio, 
o ^2 (E)/Oq+(E) is indicated by Seaton (1959) to take on a value of .2 for all energies. 
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Figure 1-D. Cross section versus energy of ionization for Nj and 01 after Keiffer 

and Dunn, 1966. 




Kazaks et. al. (1972) found a value of .25. The next two ratios of Equation (D-11) were 

discussed earlier. This author’s examination of a .(E)/a^ .(E) suggested a value of near 

O N2 


.5 satisfied the bulk of the energy range. Finally, Borst and Zipf (1970) give a value of 
14.1 for the third ratio of Equation (D-ll) as discussed earlier. In the text of this report, 
these ratios will be referred to but for simplicity of the total continuity treatment, the 
ratios will be defined as: 




V®) 






(D-12) 


where the subscript on Rj;(E) refers to the particular vibrational level of Nj*. This 
permits a rewrite of Equation (77) to read as: 


.KO«P)p = R„,(E),C39I4A) ^ 


(D-13) 


By assuming the ratios are not dependent on energy and inserting the values discussed 
above, one finds that R^^(E) = 1.41 if the historical values are retained. 


It is also referred to in the text that the production of 0’*'(^P) is related to the 
emissions arising from the first vibrational state of N^B^ S^, which leads to the X4278A 


band emission. The work of Shemansky and Broadfoot (1971) suggest that a 


N|* 


(E) is a 


factor of 3.3 lower for the first vibrational level than for the zeroth level. For this 
particular cr_,.^(E), the ratios will be given by R„ (E) and (D-13) will become 

NJK 'Jl 


7j(0+*P)p = R„j(E)t7(4278A) 


(D-14) 


Similarly, from historical values, Rj^^(E) takes a value of 4.65. 
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APPENDIX E 


DATA REDUCTION OF X7319A GROUND BASED INTENSITIES 


The intensity of 7319A and OH Pi ,3, (7316.4A), were determined in the 
following manner. From Figure 24 it was deduced that when the spectrophotometer was 
peaked at 7319A, the transmittance of 7316.4A (Tqj|) was 63% that of 7319A. When 

the instrument was peaked at 7316.4A, the transmittance of 7319A was 63% 

that of the OH (i.e., when the instrument is tuned such that = 7316.4A, then 

"^0+2 p ~ "^'max when tuned to - 7319A, Tqjj - .63 

The calibration factor (CF) was 16.15 counts per second per rayleigh. CF was 
determined with a low brightness radioactive source recently calibrated by the National 
Bureau of Standards. The instrumental integrated transmission function was accounted 
for in this caluclation. 

Instrument signal levels can be converted to intensity in the following manner. 
The equations for the two wavelength positions can be written. The equation for the 
7316.4A wavelength position can be written as; 


V7316.4 - VgG 
~ CF 




max*'!'3‘6-4 


(E-1) 


and for the 7319A wavelength position as: 


V7319 - VgG 
CF 


= T 


OH 


I 


7 3 1 6 .4 


"^max ■ '^73 19 


(E-2) 


where are the measured photomultiplier signals for the 7319A, 7316.4A, and 

73IOA (Background, BG) wavelength positions. The amount of background (V^q) was 

determined at 7310A, between two hydroxyl lines of concern. The two equations above 
can be solved for expressions of intensity, i.e.. 
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and 



(E-4) 


Recall from above that Tqh ” ^ri+ “ ^max» where T^ax determined from Figure 
23 of Chapter V. ® 

The analog charts were reduced on a Gerber digitizer where the voltage 
amplitudes were punched on cards for 7319A, 7316.4A, and 7310A wavelength positions 
of each spectral scan. 
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APPENDIX F 


Nj PRODUCTION FUNCTIONS FOR INTERPRETATION 
OF AURORAL EMISSION INTENSITIES 


The ionization rate produced by incoming electron steams has been studied by 
Rees (1963), reviewed by Rees (1969), computed by Berger et. al. (1970), and recently 
examined by Banks et. al. (1974). Power law and exponential shapes have been compared 
to measured N^INEG volume emission rates in attempts to describe the production 
versus altitude, for various electron energy distributions. 

The purpose of this study is not to analyze these electron energy distribution 
functions or flux rates. The purpose here is to give some indication of the altitude of 
maximum energy deposition and peak ionization rate. This is accomplished by viewing up 
the flux tubes along an auroral electron stream and measuring emission signatures of the 
perturbed atmosphere. 

To typify or envelope the extremes of the production rates (for N 2 ions), three 
altitude distributions have been assumed in this study. The peak ionization rate for each 
was fixed and the distribution above and below that altitude was described by their 
respective functions. The three functions are described below. 

Monoenergetic (M) 


Rees (1963) studied the production rate of NjlNEGo.o radiation as induced by 
isotropic streams (0°-80“) of monoenergetic electrons. The isotropic monoenergetic 
production rate then is that of Rees (1963), as inferred from his height versus normalized 
volume emission rate for the various energies that correspond to an altitude of peak 
ionization. Banks et. al. (1974) performed a study similar to Rees (1963) using refined 
electron scattering cross sections. Their results agreed with Rees (1963) for E > IKev 
primaries but lower energies peaked slightly hi^er in altitude. The Rees (1963) results 
are used here. 

Nj Proportional (P) 

The N 2 proportional shape assumes that the emission rate is proportional to the 
N 2 number density above the peak of maximum ionization and drops to zero at the 
peak. 

Chapman (C) 


Chapman (1931) described a theory of electron production by photoionization, 
which has classically been termed Chapman’s production function. This function and its 
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many forms have been reviewed in many texts including Rishbeth and Garriot (1969), 
Chapter III. The form used here is for normal incidence to the atmosphere and is given 
by 


P * Po exp {l - h - e‘^ > , (F-1) 

where 

h - reduced height = / -J “ 


and 


z • altitude (km), 

Zo - altitude of maximum production, 

H - scale height = kT/mg, 

Pq - maximum production rate at peak altitude . 


Figures 1-F and 2-F show altitude profiles of the production rates versus altitude 
for the three functions described for peak production at 120 km and 200 km, 
respectively. The peaks are normalized to one and the intention is to convey the 
distribution of the rate above and below the peak for the respective functions. 

When viewing an aurora from below with these types of production profiles, note 
that the peak of the Chapman profile corresponds approximately with the median point 
on the profile. For the monoenergetic and Nj proportional types, the median is 
somewhat above the peak production point. At 200 km (Figure 2-F), the three profiles 
are shown to diverge to larger extremes in their drop off rates from the production peak 
altitude. 
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Figure 1-F. Altitude versus relative production rates with peak production 
at 120 km for a monoenergetic, Nj proportional, and Chapman 
type production functions. 
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Figure 2-F. Same as Figure 1-F, ej 
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APPENDIX G 


DATA REDUCTION OF NASA 4.329 PHOTOMETER SIGNALS 
FOR INFERRED 0+(* P-2 D), NjlNEGo,©. AND 
N3lPGs,3 EMISSION INTENSITIES 


G.1 INTRODUCTION 

The NJINEG system is a fairly narrow band system. The effective transmission of 
the 3914A filter was determined to be 30%. This effective transmission was determined 
graphically using a 500°K laboratory spectra of the NJ1 NEGo,o band. No attempt was 
made to compensate for a changing rotational distribution, which might result from 
changing thermospheric temperatures over the measurement range. 

Figure 1-G depicts the reduced intensity of the NtlNEGo.o emission (solid 
curve). Note on the figure that the altitude at which energetic electrons were encountered 
is at ~ 1 20 km as indicated by the particle detector. Below this altitude, the photometers 
are viewing aurora above which is either weakening during the payload ascent, or the 
payload is passing into a weaker form during flight. At 150 km, an enhancement is 
encountered which nearly doubles the emission intensity. The rilNj ] • or effective 

columnar density profile shape is shown in the dashed curve (no scale indicated). The 
purpose of this curve is to suggest that this particular observation departs significantly 
from an ideal, steady state situation. 

G.2 0+(2p-2D) AND Nj1PGs,3 DATA ANALYSIS 

The introduction to Chapter VIII discussed briefly the difficulty in discerning the 
0+(2p.2D) emission intensity which is surrounded by N 2 IPG emission. This is 
particularly complex in broad band photometry where the spectra shape is not resolved. 
It was for this reason two channels (7326A and 7354A) were dedicated to the 
measurement. The data analysis requires a deconvolving of the emissions to discern their 
sky brightness. 

G.3 EFFECTIVE CHANNEL TRANSMISSIONS OF O+l^P-^D) 

The filter transmission profiles are shown in Figures 2-G and 3-G for the 7326A 
and 7354A channels where the solid curves represent the measured transmission profiles. 
Examination of the ©’’'(^P-^D) wavelength positions at 7319A and 7330A suggest a 78% 
transmission for both wavelengths through the 7326A filter. The 73 54A filter allows 10% 
transmission of 7330A and 3% of 7319A emissions. 

The effective transmission of ©'•'(^P-^D) (i.e., T^+) was determined by weighting 

the transmission of the individual volume emitting probabilities. The effective 
transmission relationship for a given filter can be written as; 


140 



Figxire 1-G. Altitude versus intensity of N^lNEGo.o is indicated by the solid line. 
Particle altitude is noted at 1 20 km. Dashed line is columnar Nj density 
(no scale indicated) for purpose of showing the changing aurora. 
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+ 


(G-1) 



A 7 3 1 9 .nr 
A7319 

A 7 3 19-3 0 


^7 33 0 
A73 1 9-3 0 


• T 7330 


9 


where has been previously defined (Appendix C) and represents the filter 

transmission at the line position of the particular emission. The effective transmission for 
the 7326A channel was .78 and for the 7354A channel was .’068. 

G.4 EFFECTIVE CHANNEL TRANSMISSIONS OF Nj IPG 

The effective N 2 IPG transmission for each channel was determined by modeling 
the Nj IPG band system and finding the transmission of each rotational line of the band 
through a given filter. The NjlPG band model description, verification with laboratory 
spectra, and Fortran IV listing are in Appendix H. For each of the twenty-seven substates 
of the Nj IPG band, twenty-nine rotational lines were computed and their effective 
transmission found. This process was performed for temperatures between 200°K and 
1000°K, in 50°K increments to account for substate and rotational distribution changes 
effecting the band shape over the temperature range anticipated in the atmosphere. 

Due to the number of transmission calculations involved, it was analytically 
expedient to fit the two filter transmission measurements and represent them analytically 
(the dashed curves of Figures 2-G and 3-G). The analytic representation permitted a 
transmission calculation for each line position versus a table look up and interpolated 
value from the measured profile. 

To give perspective to the reader, the N 2 lPGs ,3 band is plotted in Figure 4-G for 
rotational distributions corresponding to 2dO°K and 700°K. The 0'*’(^P-^D) wavelength 
positions and the peak bandpass locations of the photometer channels are also given. 

The N 2 lPG 6,4 band (7274A) also provides a background to the filters but to a 
lesser degree. Shemansky and Jones (1968) discuss relative band intensities and 
vibrational populations. Upon examining their discussion, the V=6 level was assumed to 
be populated 55% that of the V=5 level. This is further examined in the relative cross 
sections listed in Shemansky et. al. (1972). Considering the transition probabilities of 
N 2 lPGs ,3 and N 3 lPGe ,4 bands (Shemansky and Broadfoot, 1971), the intensity of the 
(6,4) was estimated at 16% that of the (5,3) band. The vibrational temperature range 
encountered in the atmosphere was not taken into account in this ratio, i.e., 
I1PG6,4/I1 PGs, 3 = .16 = Constant. 

Plotted in Figures 5-G and 6 -G are the effective transmissions of the pertinent 
N 2 IPG bands for the 7 3 26 A and 7354A channels respectively. The temperature effect of 
the band on the effective transmission is most notable for the T < 600°K portion of the 
7326A channel. Although the N 2 lPG 6,4 transmission is ~ .3 the N 2 lPGs ,3 transmission 
for the 7326A channel, the N 2 lPGe,4 band strength is only ~ .16 that of the N 2 IPGs , 3 - 
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Figure 4-G. A 700°K and 200°K synthesis of N 2 IPG 5 ,3 with a 3 A triangular slit 
function. The rocket filter transmission peaks and the 0'*'(^P-^D) 
emission line positions are noted in the figure. 
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Figure 5-G. Effective transmission of N 2 IPGs ,3 versus temperature. The point plots 
are transmissions calculated for a temperature with the model using the 
analytic function plotted in Figure 2-G (7326A channel). The dashed 
Unes are linear legs used to represent the calculated values and values 
in between calculated points. 
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Figure 6-G. Same as in Figure 5-G, except for the analytic function plotted 
in Figure 3-G (7354A channel). 
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G.5 REDUCTION OF 0+(2p-*D) AND NjlPGs.a INTENSITIES 

This section will describe the deconvolution of the two emissions from the two 
(7326A and 7354A) channel signals. 

The assumption was made that the only signal contribution to the two channels 
was due to the photon contribution resulting from the 0'*'(*P-*D) and NjlPG emissions 
and the photomultiplier dark current. The dark current was measured on the payload and 
had been removed from the signal prior to plotting Figures 37 through 39 in Chapter 
VIII. The remaining signal may be equated to the sky brightness of the emissions 
(Intensity, i.e., IO‘'‘(*P-^D) for 0‘*'(^P-*D) emission and IIPG for Nj IPG band system in 
rayleighs), i.e., 


JV 

SF 


^0+(2p-^D) ^0+2 P ^ 


(G-2) 


where 


V — the signal measured in a given channel (volts), 

SF — sensitivity factor of the given channel, volts rayleigh"’ ° A‘* 
(Table 1-G), 


To+2p’TlPG 


— effective transmission of the system (subscript) for a given 
channel. 


It was recognized in the previous subsection (G.4) that the first positive 
contribution was actually made by both the (5,3) and (6,4) bands. Equation (G-2) can 
then be rewritten as: 


_y 

SF 


^OVP-^D) ^O+^P ^ ^ IPGs ,3 


IPG 


S ,3 


+ I 


IPG, 


6 ,4 


IPG 


6 ,4 


(G-3) 


where the subscripts refer to the particular NjlPG band of a channel. It was further 
stipulated in the last section that the I1PG6,4 band is weaker than the I IPGs ,3. Since a 
channel was not available to monitor its contribution to the other channels, it was 
assumed that 
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TABLE 1-G. CALIBRATION SENSITIVITIES FOR NASA 4.329 
ROCKET PHOTOMETER 


Channel 

Sensitivity* 


(Volts per Rayliegh) 

3914A 

.24 

7326A 

.027 

7354A 

.025 


♦Calibration error ±.003 


4PG, 


6 ,4 


= .16 I 


IPG 


5,3 


(G-4) 


Equation (G-3) now can be written as: 


SF 


I 


0+(2p-2D) 


• Tq+^p + IiPGs.3 ^'^1PGs,3 


+ 


.16 T 


IPG 


6,4 


(G-5) 


A cursory examination of the Nj IPG effective transmissions of Figures 5-G and 6-G 
suggest Tjpq ^ < .3 Tjpq^ ^ or that the second term in the parenthesis of Equation 

(G-5) is 5% the first term. 


Equation (G-5) can be written for each channel (7326A and 7354A). The 
respective channel signals (Figures 38 and 39), sensitivity factors (Table 1-G), effective 
0'*'(^ P-^D) transmissions (Section G.3), and effective Nj IPG transmissions (Section G.4 
and Figures 5-G and 6-G) provide the necessary input for Equation (G-5) for both 
channels. The two channel equations (Equation G-5) can then be solved for the two 
emissions algebraically to yield: 


4PG 


5 ,3 






SF 


B . xB 


IPG 


+ .16 T^, 


5 ,3 


IPG, 


:„)i- 


0+2 P 
-B 


) 


IPG 


S,3 


'0+2p 


+ .16 T 

pB 

0+2 p 


IPG, 


J 


(G-6) 
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and 


I 


1 


0+(^P-*D) 


tA 

^0+»P 


[5X-’iPG...( 


T® + 
*1PGs,3"^ 




.)] 


(G-7) 


where superscripts A and B refer to the noted values of the 7326A and 7354A channels 
respectively. 

The thermospheric temperature, as indicated by model atmospheres, changes by 
several hundred degrees over the altitude range of this measurement range. In Section 
G.4, the sensitivity of the effective transmission of the Nj IPG band system to 
temperature in both photometer channels was discussed. In order to account for the 
temperature effect, it was assumed that the bulk of the Nj IPG emission was originating 
from a half scale height above the payload for a given measurement altitude. It was 
further assumed that the effective temperature dictating the distribution of the N 2 IPG 
emission was that of the Jacchia (1971), = 1000°K, model atmosphere. The nitrogen 

scale height was determined for a given measurement altitude, divided by two, and added 
to the measurement altitude to provide the reference temperature from the model 
atmosphere. The possible error induced in this method is discussed in the error analysis. 
All other inputs to Equations (G-5) and (G-6) were constant or were values at the 
measurement altitude. 

The results of the data analysis of NjlNEGo,© are plotted in Figure 41 of 
Chapter VII. The data is plotted only for the portion of the ascent in which the payload 
was in the aurora and not when it was being viewed from below. The solid curve drawn 
through the P-^ D) data points is a hand drawn fit to the data. 

G.6 ERROR ANALYSIS OF TH^^ 0+(2p-*D) PROFILE 

A number of potential errors in the resultant C^C^P-^D) profile arise from the 
uncertainty in the quantity inputs of Equations (G-5) and (G-6) outlined in the previous 
section. The method used in treating these errors was to best estimate the potential 
deviation of the input parameter. The deviated value for that parameter was then input 
into Equations (G-6) and (G-7) with all other input quantities remaining fixed. A new 
Q+( 2 p. 2 j)) intensity profile was then computed, and the percent difference the resultant 
profile departed from the original 0"''(^P-*D) profile was computed. 

One perturbation was made by increasing the model atmosphere temperature 
profile by 100“K, reflecting the rotational distribution or band shape effect of the 
Nj IPG system. Additional perturbations were introduced by shifting the peak 
transmission of the photometer filters by 5 A, adding an additional 1% transmission to 
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each photometer for N 3 IPG, and applying a ±12% uncertainty in the sensitivity factor. 
The percent deviations the perturbed profiles deviated from the original for each of these 
are listed in the first four columns of Table 2-G. Note that the model temperature 
assumption and the current leak error lead only to possible lower values of 0 "'‘(^P-^D) 
emission than deduced. It was felt by the author that the exospheric temperature profile 
would not be less than the model assumption used in an auroral situation. Also note that 
the entire temperature profile was increased by 100°K for this error perturbation, rather 
than just the exospheric temperature. Ambient temperatures will tend to shift 
transmission peaks on filters. This reflects the second error perturbation. The third error, 
the current leak, accounts for the possibility that other emissions not treated here could 
be leaking through the filter wings and contributing to the signal. The sensitivity error 
results from the calibration source uncertainty (Table 1-G). 

The statistical error or random error due to shot noise, counting statistics, 
transmitter noise, etc., is given in the last column of Table 2-G. As stated earlier, the data 
were reduced using the smoothed curve through the data points of Figures 37, 38, and 
39. The data were also reduced using the actual measured values of Figures 37 through 
39. The statistical error or scatter was determined as the deviation of the scattered input 
results to the measured 0‘*"(^P-^D) profile of Figure 41, where the curve in Figure 41 
represents a hand fit smoothing of the results. The root sum square of all errors is 
displayed in the last column of Table 2-G. 



TABLE 2-G. PERCENT ERRORS DEDUCED FOR 0+(2p-2D) EMISSION DATA 
FROM NASA 4.329 ROCKET PHOTOMETER 


Z(km) 

Model Temperature 
(T,T+ 100°K) 

Filter Drift 
±5A 

Current Leak 
TiPG = TiPG + .01 

Calibration 

±12% 

Statistical 

Errors 

V2(Error)2 

120 

-30% 

±7% 

-18% 

±36% 

±25% 

+44%, -56% 

140 

-12% 

±7% 

-12% 

±30% 

±20% 

+36%, -40% 

160 

-9% 

±7% 

-9% 

±24% 

±20% 

+32%, - 34% 
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APPENDIX H 
Nj IPGj j BAND MODEL 

The N 2 IPG bands are very broad and contain twenty-seven branches of comparable 
intensity. The energy equations for the B^TTg state are given by Budo (1935, 1936) and 

for the state by Schlapp (1927). The line strength factors as theoretically 

determined by Budo (1937) for Hund’s case (a) and (b) coupling were used. The 
molecular constants used were those summarized by Shemansky and Jones (1968). 

The model was verified by running a 3 A triangular slit function through the (1,0) 
band and comparing it against the data of Vallance Jones given by the Shemansky and 
Carleton (1969). The synthesis is shown in Figure 1-C where the B^TTg rotation 

equilibrium is depicted against a near 400°K spectra. Some slight mismatch is evident but 
considered insignificant for purposes of this analysis. Nj IPG spectra in aurora has been 
discussed by Shemansky and Jones (1968). Their spectra suggests the N 3 IPG is in the 
ground state distribution (‘ Sg). Transition probabilities for first positive, as measured by 

Shemansky and Broadfoot (1971), are too large to permit rotational relaxation to the 
B^TTg state prior to emission at collision frequencies corresponding to altitudes above 90 

km. A Fortran IV program listing of the model follows Figure 1-H. 


153 


Ol 



Figure 1-H. Dark line is 400° K spectra of Vallance Jone’s as shown in Shemansky 
and Carleton (1969). The dashed line is a synthesis of Nj IPGi ,o as modeled 
utilizing a 3 A triangular slit function. 


tev- 


FORTRAN S3URC: LIST 


1 DIMENSION FGI3CI 

Z DIMENSION ARU(30)»AQll(33),APl;(33)>ARI.2O0)tA912t 33I>AP12( 3?). A’. 

H3(30»,A013(30),APl3I30ltAR2U30»fAQ2H 301 ,AP2U 331, AR22I i0i,t,Q2ii 
i3j|, AP?2( 50, AP23I30I,AQ23I30» ,AP23I 30 1 , AR3 U 33 ) , AC,3U 30), AP31 I 30) 
3,AR32I -O) ,'.i32I30) , AP32 ( 30 ) , AR33T iO),AC.33( ^0),AP33( 30) 

3 01 MENS I UN MJK 30 ,FU2< 30) ,FU3(30) ,FLH30),FL2I33» ,EL3( 30),Pll« 3C ) , 

13111 3C),Kli(3( ),P22(30) ,022(30),R22(33) ,P3 3( 5-),03 3(30),R33I30,?1 
2200 , .12 (■>.)) ,R 12(30) ,P13 150) ,03 ( 30) ,iU3l iO),F2l( 30), v;2l( 30, R21< 
33C ),P3l(3C),031(30) ,R31 1 39) , P32 ( 3u ) , U32 ( 33 ) , R32 ( 3D ) ,P2 3( 30 ) , W2 3( 3 . 
A),'<23( 5U),XJl(3O,XU2()0),XlJ3( J0) 

A DIStNSC^'N XPll (2J),XwU(30),XiUl(30),XP22( 30),X322( 30),XR22(3C),xP 

135( 30 ,XCi33( U- ) ,X<33(;53) , XP13( ’0 ) ,XCil3C iC) , XRI ii 30 ) ,XPI2( iO ),X(j‘. 2( 
2 3.. ) ,X.U2( tO),XP23(3.':) ,XU23(3D) ,XR23(30 ) , XP3 1 ( iO ) , X ..3U 30 ) , XP.3I ( 50) 
3,X->32( J0),Xf.;52(.'0),XR32(3O,XP21(33),X0<l( 30),XR21( 30 
S l)l“i.-NSI JN XP( 3,3, 30 ,X.^(3, >,3: ),X;U3,3, 5 5) , VP ( 3 , 3, 3.0 ^ yQ( 3 , 3, .1 ), Y 

U( 5,3,3., ),PP( 5,3, it/) , O'. (3,.', 30), RRt 3, 3, 30) 
b OlKsNSlUN ?.P(3,3,3..),SJC3,3,3.'),S:U3, 5,30) 

C C 'INSr-'..NT;> f .Ik rue UPt'ER ANO LO-nEk STATcS kcRi. GIV; N bv Slll-MANSKY 

C CONSTANTS f.'jP. THE UPPtR A.NU LOWcR STATcS WERE GIVEN BY SH.:ManSKY 

C AND jnriESC I >l.t5,PLANE T. SPACE SCI . ,PP. 1 1 15) 

7 R = KjO . 

10 SL=2.9'(/7f>M-10 

11 bX = 1 ♦ 3..U 3 5t - 1 6 

12 M=o.2234E-27 

13 E=l"»SL/rtK 

iA bLiJ=l. 6374b 

13 0ei=1.4545 

lb AU=. 017)4 

17 AL’. 0179b 

.It, ,;al=-.i)03 

2a GAH=-b.44E-05 

11 Uc.J = 3.bC,6-06 

13 btru=3.Tt-0e 

14 <i co iriNub 

b READ UPPER AND LOWER VIBRATIONAL LEVELS AND TEMPERATURE 

25 READ(5,17IVU,VL,T 

.0 17 E3.<MAT(E3.u,F3.0,lX,F5.0) 

aT ,.<lTt (6,200 

•,U 200 FOkf-'AmHl ) 

51 (T=T 

;2 IF(IT.:;G’.)99)G0 TO 99 

G COMPUTE WAVcNUKbERS OF VIBRATION 

j 5 TEtl=5962l.‘31 

3b TEL’.SOIOS. »6 

^7 WU=I735.92 

40 WXU=l4.39N 

41 WYU=-.v04Ol 

*•2 •/I2u*.0f'0513 

43 WL =146(-'. bl 

44 WXl*l3.b5l 

45 wYLs.C3o25 

46 W/L=.Ot, 17:> 

4 7 WlU = Tei,* V« )♦( VU*.5)-WXU«‘(VU«.5)**2.»WVU*(Vj4.5)**3.-WZU'-( VU+.5)**4, 

‘ w a =TE1 »M 5( V1.*.5)-WXL*( VI 1.5)*«2.«WVL*(VL».5)**3.-W/L»( VL«-.S)«"<=4. 
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II ■■■ll■lll■■ I 


51 OWN»WNU-W.'<L 

52 WNA31./(DnN*1.0E*08) 

53 WRITE (6,7) VUtVLtTtMNA 

54 7 FORKAT(lX,4H8AN0,F3.0fF3.0,10X,5HrEHP*,F5.D«10Xtl3MBANO HEA 0 >,F 6 .J 

1) 

C COMPUTE Vlu distribution OF GROUND STATc 

55 TPI=6.2631 

56 El-E/TPI 

57 WG=235').61 

60 WXG=14.456 

61 WYG=. 00751 

CGY>MG*( VU«.5)-WXG«( VU«-*5 )**2.4-HVG«(VU«'.S) **3. 

6i SG*0.0 

64 DO 2 1=1,10 

65 V=I 

66 X=V*.5 

67 GV = WG*X-WXr.*(X**2. )♦WYG•(X*♦3. I 

70 G*£XP(-GV*f 1/T) 

/ I SG = SG 

72 2 CONTINUE 

,W G'< = lOOO.'<'£XP(-GGV*tl/T)/SC 

75 WRITE(6,73)GN 

76 75 F0RMAT(9X, 17H VIBRATION I NTEN« , 1 1 0. 3 1 

C COMPUTE ROT energies IN GROUND STATE 

C A 1 IN COL. 20 WILL GIVE ROT DISTRIBUTION OF cXCITEO STATE 

77 BVU = BEU-AU*( VU<-.5) 

10;> BVL=B£L-AL*( VL + .5) 

ICi ALL = -l.335«-.007*VL 

102 YU=25.6B+.253*VU^.u028*(VU**2.I 

lOi DVU=DEU*B: TU*( VU+.5) 

104 aec= 2 .uio 

1.5 A£G=.01i,7 

*(,.6 OVG = ( 5.8-.001*( VU<-.5))*1.0E-06 

1j 7 BVr.=BEG-AEG*( VU>.5) 

110 DO 19 J=l,30 

111 A=J-1 

112 EG( J) =bVG*A*( A*l. )-DVS*A*A*( A+1.)*«2. 

113 19 CONTINUE 

115 WR1TE(6,13)(FG(J),J=1,1S) 

122 13 F3RMAT(iX,15E7.0) 

123 WR1TE(6,13)(F6(J), 0=16,30) 

C COMPUTE ROT ENERGIES OF LOWER STATE OF TRANSIT1.0N 

C THE ROT LEVELS OF TRIP SIG STATE ARE GIVEN BY SCHLAPP( 1937,PHYS. 

C REV.) 

130 DO 3 K=2, 30 

131 0=K-2 

i?2 ALr=ACS(ALL) 

133 SK1=( ( (2.«D»3. )**2. )*(3VL**2.l+( ALT**2. )-(2.*ALL*8VL) )**.5 

i04 SR11=ABS(SR1) 

155 FL1(K)=BVL*D*(04-1. H-(2.»0*3. )*bVL-ALL-SRll*GAM*(D^l. ) 

13o 3 CONTlNUt 

lAO FLim=C.O 

14 . DO 4 K=l,30 

i‘-,: L) = K-1 

!•- 5 EL. (K. ) =bVL*J*'(f>+i. ) 
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i 


144 

14b 

147 

iiiO 

ISi 

lt'3 

it>4 

ir»5 

i oO 
loi 

L 

C 

c 

lb 
i :o 
*64* 
16b 
i bb 
ib7 
i ro 
l/i 
17c 

1 

;7b 

.7o 

177 

£ >■ 1 

^'4 

b 

^06 

2 10 
ill 
2^. 2 
2 13 
4 14 
^-15 
217 
2^0 


221 

d22 

2^3 
2<-4 
22 5 

22b 

22 7 
. :nj 
2 :,! 


4 CONTINUE 
00 5 K^l, 50 
D = K 

ALT»ABS( ALU 
PIP = ABS(2*»1)-IU 

SR2 = ( (P1P*^«2. )*(bVL*«2.)*(ALT**2. )-C2.»A\L»ki/Ln***5 

5K22 = Al^S( SR2) 

FL3(K)=faVL*J»(D+i. )-(2.*D-l. ) ♦B\ZL-ALL^SR22-CiAM*0 
b CONTINUr, 

tJLi=bVL«»2. ♦ALT**2.-2.*ALU'BVL 
RR« = AbS( 5ii ) 

FL3( 1 )=LVL-ALL-(RKR^*.5)-r,AM 

COMPUr- RJT GNFRGILS UF UPPER STATE OF TRANSITION 

THl rot LiVELS OF TRIP PI STATE ARE GIVEN BY BUOO( 1V35» 1936t Z . 

PHYSIK) 

00 U J = l»'30 
U = J-l 

Yi=YU*lYU-4, )+4./3, 

Y2=VU4(YU-l. )-4*/9. 

YY = ( Y2-2.^U4(l)^l. n/(Y1^4.<‘D*(D^l. ) ) 

YZ = ( Yl+4.^0*(D>1. ) ) 

YY7 = (Al>S( Y2))*«.5 

FUK J )-liV;J4(L.«( D+U )-YYZ- <2. ♦YY/3.) ) 

U CGMiNut 

ihj 10 J^i , 3 ; 

0:- J-l 

Vl=4,/ 

Y?=-4./V, 

YY:=(Y2-2.<=U^<0 + 1. ) )/(yi>4.*D^(n^l,n 
FU. ( J ) =BVU4( ) ♦•(4.«YY/3. I ) 

1 . CONTI'iUE 

UJ B J=l,30 
0 = J-1 

YI=YU*( YU-4. 

Y2=^YUt(YU-l. 1-4. /9. 

YY = { Y2-2. <'0^' ( 0^1. ) » / (Yl4-4.+0*CD+l. ) ) 

Y2 = ( Y I ^4. *(:*( U + 1. ) ) 

YY2=( ABS( Y2) )**.5 

FU M J ) =BVU#(U* ( 0+1 . ) <-YYZ- {2.4YY/3. ) ) 

e cu-niNuc 

L)WNL.= 0WN 
L)U'i2 = OWN 
0 W N 3 = 0 ^ 

CJU^PUTt LINE STRENGTHS FOR 27 BRANCHES 

LINE STRENGTH FACTORS ARE GIVFM BY BUOO (1 7 37 » Z . PHY S I K ) 

L>0 70 K = 2t30 
B = K-l 

Ui=(YU*( YO-4. )+4.*0**2. )*».5 
U 'i = < YU-f ( YU-4. )+4.<=(0*l.)*<‘2.)*».5 

Cl = J# ( . ) ♦YU*! YU-4. )*2.*(2.*n+l. ) ♦ ( U- 1 . ) *0* ( 0+ 1 . ) 

C2=YU*< YU-4. ) ♦4. *U* ( D*1 . ) 

C 5 = ( o-l. ) 4 (0+^ . ) *VU*( YU-4. J *2. ♦( 2. *0*1 . )*D* (0*1. )♦( 0*2. ) 

7 I P = YIJ* ! YU-4. ) +4 . 

/ I I =Ai; :S( / 1 IM 

C=<nr.( yu-2 • )4 ( ( / I T ) ♦•' . 5) *-YI)4 ( Y'J-«». ) >4 . 
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233 

234 
^35 
2 36 
^37 

4.^0 

241 

2H3 

<^44 

<:45 

^60 
2:^1. 
c -»2 

2t -> 
^*>5 

c ‘J 6 

r‘-.7 

^:o0 

2o 1 
<:o i 

i o4 
463 
460 

467 

4?0 

4?i 

2^2 

2Y-k 

4/4 

476 
4/6 

477 

300 

j>: 1 

3:2 

30 3 
36 3 

3 J 6 

■507 

.-.10 

m 

^12 

ai3 

3 1 '* 
;l'j 
5i 0 


C1*C1+CRU 

C3=C3<-CRU 

0IL=U**2. 

AUIHK)=( ( (BIL+O-I. )*Ul*< YU-2. >*2.*0*(0*0-l. » )**2. )/(0«Ctl 
Au22(K. ) = ( R.<=(2.*0+1. )*(0*0«-D-l,)**2.)/«D*«0*l. >*C2I 
AJ33< K) = ( J (D*0»U-1. »*U3-C YU-2. )*( U>2 .)) **2 .)/((&♦ 1 .) *C 
13) 

AJ12(K) = ((2.*D + 1. )•( <13*0+0-1. )«( YU-2. »*un**2.)/(0*»D+l.)*CU 
Arjl 3 (K)=( ( ( 0* 0+0-1. ) *Ul+( YU-2. )-2.*(0-l. »*(0+i. »**^. »**2. )/( (U+1. ) 

l*Cl) 

A021<K)=<2.*t0*( YU-2. )-2. )**2. »/(U*C2» 

31 1 = ( ( 0*0+ l)-l. )♦{ YU-2. )-U3) 
dI2=ABS<0I 1) 

Aa32(K) = ( <2.*0+l. )*(BI2**2.) )/(0*(U+l. )*C3) 

AU231K)=(2.*< (0+1. )*YU-2.*0)**2. )/((D+l.»*C2) 

3J1=( D«0 + 'J-1. )*U3-( YU-2. )-2.*D*D*(0 + 2. ) 

DJ2=ABS(BJ1) 

A03l(K) = (i5j2**2. )/(0*C3) 

70 CUHTlfjUu- 

00 7i J*2.30 
0 = J 

Ul = ( YU*( YU-A. )+<V.*U**2. >**.5 
U3 = ( YU*( YU-^.. )+A.*(D + l.)**2.)**.5 

C1*D*(U+1. )*YU*( YU-A. »+2.*(2.*0+l. )*(0-l.)*0*(0+l. » 

C2=YU*( VU-i. )+A.*D*(U+l. ) 

Ci=(D-l. )*(0+2. )*YU*(YJ-A. l+2.*(2.*0+l. )«0*(0+L. l*(U+2.) 

2Ir> = YU*( Y0-A.‘)+A.*U*U 
21 T*ABS( 21R) 

ORiJ=< YU-2. )* ( ( 2lT)**.a»-YU*( YU-A. )+A.*0 

Cl=Cl+CRU 

:3=C3+CKU 

BlL = Li*<4. 

ARlll J)»( ieiL-i.)*( (0+i:)*Ul-YU+2.*BIU**2. )/(Cl*«2.*0-l.»» 

M.i22(J) = (B.*(U-l. )*(D + 1. )**3. )/(U*C2) 

A«33( J) = ( (O+l. )*< (0*0-1. )*U3+(D-l.)*(YU-2. ) +2. *0*0* ( 0+2 . ) 1**2. 1/(0 
1*(2.*0-1. )*C3) 

iR12( J)*( (li*0-l. )*( (0+1. )*(YU-2.)-Ul)**2.)/(0*CU 

AK1.3( J) = ( ( (0-1. )**2. )*( 0 + 1. )*( (13+i. )*Ul-( YU-2. )-2.*0*( D+1.U**2. )/ 
1(D*(2.*0-1.)*C1) 

AR2K J)*(2.*(D*0-1. )*YU*YU)/((2.*0-l.)*C2) 

AR32( J»*( ( 0*0-1. ) *(U3+(D+1.)*( YU-2. n**2.)/(D*C3) 

AR23( J)=(2.*( D+1. ) *(D*( YU-A. 1+2. )**2. ) / ( 0* ( 2 . *0-1 . )*C2) 

AR3K J) = ( (U*U-1. )*( (0+1. )*U3+( YU-2. } -2 . *0* ( 13 + 2 . ) ) * *2 . > / ( ( 2. *0- 1 . ) * 
1C3) 

71 COnUNUE 

D:J 72 L = 2,30 
U = L-2 

Ul = ( YUM YU-A. )+<».*U**2.)**.5 
J3=(YU*( YU-A. )+A.*(U+l. 1**2. )**.5 

01=0* (0+1. )*YU*( YU-A. )+2.*(2.*U+l.)*(0-l.»*U*(U+l. I 
C2 = YU*( YU-4. )+A.*D*(0+l. ) 

C3 = (0-l. ) = (0 + 2. )»YU*(YU-A. ) + 2.*(2.*D+l. » *0* ( U+ 1 . ) • ( 0+2 . ) 

/li' = YU-M YU-‘t. )+A; -i>0*u 
2 1 I = \0:>(211') 

V. Y ..- 1 r ) . s ) -VU5- ( Yu-4. )4<t. *u 
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317 
32 0 
j21 

322 
3;* 3 

324 

325 

:J26 
32 7 
330 
351 

332 

333 
3 54 

335 

336 

337 
340 

342 

343 

344 
545 
^46 
:>47 

J 30 

351 

352 

353 

354 

355 

356 
:>57 

361 

362 
36 3 
jj64 

365 

366 

367 
370 
->71 
i/2 

373 

374 
3/5 

376 

377 
400 
421 
4f 2 
4. 5 
4 ' 4 


Cl«Cl+CftU 

C3*C3^CRU 

APlI(t)^(04rC04(04>2« |«014CD42.)*(7U-2* ) )4*2 

!• )/( (0>1.)4C2*«0^3. )«CIt 
AP22<U = (6.*0*D*D*i0+2-. n/l(0+l.)*C2) 

BK1=LUMJ3-C YU-2. )*2.*D^ (042.) 
bK^=ABS(LKl) 

AP33(U = (0<‘(Of2. )*(BK2**2.n/( C2.4U43. )*C3) 

AP12(U = (04(D+2. )*(D*(YU-2.)4U1)*42. )/( (D4l.)4^Cl) 

AP13(L)=(0*CO+2. )*(O^UI4YU-2.*040)**2. )/(( 2.*043.)4Cn 

AP2l(U = (2.*0*( (U4i. )4YU-2.4(2.*043.) )**2. )/( (04l. ) * ( 2. «Q43. ) *02 I 

BLl*04( YU-2. )-U3 

BL2^A6S(0L1) 

AP32(L) = (i)<‘{042.)*(BL2^42«) )/( (D4i. )*C3) 

AP23(L) = (2.4U*(042. )♦YU*YU)/C C2.^U43.)*C2) 

BMi=D4U3-( YU-2. )-2.*D4CU4l.) 

0M2=AB5(oMl) 

AP31 ( L) = Cl)<*C (042. )**2. )*( BM2**2.) )/( (0+l.)»(2.*043. )*C3) 

72 COrjTINUir 

CU>«PUTe WAVc NUMBERS FOR 27 BRANCHES FOR 29 ROT. SUQSTATES 
OU 21 L=2,29 

I*L4l 

N = L-l 
M=L41 
K = L-l 

Q22(L)=0WN2tFU2(L)-FL2CU 

XC^22( L )=R*AU22(L) *cXP(-FG(L)*E/T) 

R35(L )=UWN34FU3(M)-FL3(t) 

XK33( L)=R*AR33{L)4EXP(-FG(I)*E/n 
02KL )=UWN2 4FU2CL)-FLI (U 
Xr3cl(L)=R4Aw21(L)*. XP(-FG(L)4«£/T) 

R3: (L )=DWN34FU3(M)-FL2 CD 

XP32C L)=Ri*AR32(L)*EXP(-FGm*e/T) 

02 3(L )-bWN2*FU2(L)-FL3(U 
Xu?.i(L)=R4'AC23(L) ♦cXP(-FGCt)*£/T) 

R31U )=0WN34FU3(M)-FLi (U 

XR31 ( L)=R*AR3l (U*tXP(-FC(l )*E/T) 

aiUL )=UwNi4FUl (U-FLUU 

X011(L)=R*ACll(U«£XP(-FGCU*fc/T) 

R1HL)=0WN1^FU1CM)-FL1 (L) 

XKIK L >=R4AR11(L) ♦EXPC-FGC I )♦E/T| 

R22CL )=0WN24FU?(M)-FL2(U 
XR22(L)=R4AR22 (U ♦HXP(-FG(I )♦E/T) 

;*12(L )=DW\ll^FUl (L)-FL2(L) 

X012( L)=R*A ,^i2(U ♦cXP(-FG(U^E/n 
^12(L)=0WNl4FUl (M)-FL2(L) 

XR12(L) = R4AR12(L)4fcXP(-FG(n*E/T) 

U13(L)=UWNUFU1 {L>-FL3(D 

XQ13( U=R>«'A013(L)'PtXP(-FG(LM‘e/n 

R13(L )=0WN14FUI(M)-FL3<L) 

XR13(L)=R«AR1 iiD^'tXPC-FGC n*E/D 
R21( L)^0WN?«-FU2(M)-FL1 (L) 

I )=p<.ak 21 (U -t^LXPi-FGi I ) ♦E/n 
R^ i(L) =JwN2^FU2(M)-FL3(L) 

KlC* 3 ( L) ^ R ^ A ^ / 3 U ) 4c XP ( -F G (I ) ♦£ / r ) 
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405 

407 

4i0 

4.U 

412 

414 

415 

416 
41? 

420 

421 

422 

423 

424 

425 
4^6 
427 
4J.U 

431 

432 
43:i 
4i4 
435 
*♦ ->6 
4'*0 

441 

442 

443 

444 

*t4i> 

44 6 
447 
450 
•t5i 

452 

453 

455 

456 
-r57 

460 

461 

462 

463 

464 

465 

466 
46 ? 
470 

472 

473 
-♦?4 
47 5 

‘ilh 


21 continue 

DO 22 L=3»29 

N«L-l 
M = L+-l 
K = L-l 

P22(U=DWN2<'FU2(K)-FL2a) 
XP22(L)=R*AP22(L)*£XP(-FG<N)*E/T) 
03T(L )=0fc»N3^'FU3(U-FL3<L) 

X033( L)=R4'A0 33(L>’^‘tXP(-FGlL)*E/n 

P2KL )=uWN24FU2(K)«FLKU 

XP2l<U=R*AP2X(L)*EXP(-FGCN)*E/n 

a32(L)=UWN3^FU3(U-FL2(L) 

X032(L)=R»A032a)*cXPI-FGa)^6/T) 

P23(L)=DWM24FU2(K)-FL3(L) 

XP23(L) = R*AP2 3aM‘tXP(-FGS4^)*E/n 
03KL )^Dw.^ 5 + FU3(U-FL1(U 
X^iI(L)=R<‘A03I(L)«cXP(-FG(L)*E/n 
Pli(t)=DW 11^FU1(K)-FU (L) 
XPlUL)=R<=APmL)*LXP(-FG<N)*E/n 
P12i L )^jWMlfFUl(K)-FL2(L) 
XPl2(L)=R«AP12(L)^EXP(-FG(N)*e/n 
PI 3(U =UW.\|U1U1 (KI-FL3(U 
XPl3(U=R«APl3a)*c:XPi-FG(N)*E/T) 

22 CONTINUE 

DQ 23 L=4,29 
I=L4l 
N = L-1 

K=L-1 

P33(U=0WN3 + FU3(M-FL3(L) 

XP3i( L)=R«AP33tL)*hXP(-FG(N)^E/T) 
P32( L )=UWN3+FU3(K)-FL2(L) 

XP32( t )=R<‘AP32(U '>tXP(-FG(N)*E/TI 
P3i (L) =0WN/4FU3(K)-FLHL) 

XP3K L)=R*AP31 (L)*EXP(-FG(N)*E/TI 

23 CONTINUE 

Cl) >1 = 1 .0E + C8 
CONVERT TO ANGSTROMS 
DO 30 K=l,29 
ail(K)=CaN/QllCK) 

R1I(K)=C0:m/R11(K) 

R22(K)=i:nN/R22(K) 

012(K)=CON/012(K) 

R12(K)=CON/R12(K) 

K ) t:CDN/Q13(K) 

^13(K)=CON/R13(K) 

R2l(K)=CON/R21{K) 

R2 MK)=C0N/R23<K) 

30 CONTINUE 

DO 26 K=2,29 
Pil(K)=CON/Pll(K) 

J22(K)^CO J/022IK) 

R3^( K ) =^LfJN/R ■) U K) 
l>12(K > .0 l/Pi.?(K) 
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P13IK)*C0M/P13|K) 

500 


921(K)»C0N/Q21IKI 

301 


;<32(K) »CO:>l/R32(K> 

002 


923(KI*C0N/G23CKI 

3u3 


R3KKXON/R31CK) 

504 

26 

CONTI NUn 

506 


DO 27 K=3,29 

307 


P2?nO*Cn*4/P22(K) 

3l0 


Ci35(K)*C0 I/Q33(KI 

511 


P2MK1*C0'J/P21CK) 

512 


032IKMCON/032(KI 

513 


P23(K)^COM/P23CK) 

514 


03HKI*CU.<i/031(KI 

515 

27 

CUITTmUE 

517 


DO 28 K=4,29 

:>l\j 


P33(K)^C0.M/P33IK) 

521^ 


P32(K)=^C0N/P32(K| 

jZk, 


P31(K>-CO'l/P3UKI 

3? 3 

r’tt 

CUM ! 

325 


RtrAUl 5»15 > ) SW 1 1 ROT »N£ND« A 

350 
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P0.^MAT(2X,F6.1,12Xf 1 L92X»l3f2XfF6.0) 

531 


IF( IROT.Eu.O) GO TO 310 

5-4 


O'J 311 J»2#30 

335 


K»J-1 

5i6 


L*J^1 

5i7 


XPIK J)«R«APll ( JMcXPC-FOlIKI^F/ri 

540 


XOIK J)=R*AC,11C J)*eXP(-FUlCJ)*E/TI 

5s i 


XRIUJ)=R4AR11 (J)*cXPI-FUl CU«C/T> 

342 


XP22C Jl=R*AP22(J)*eXP(-FU2CK)*e/TI 

543 


X022I J )*R<‘AU22(J)^tXP(-FU2UME/T) 

3^4 


XR ?2( J)^R^AR22U)*LXP(-Fu2CU*E/TI 

545 


XP3 H J)=R«AP33CJ) ♦EXP(^FU3<KI#e/T) 

546 


Xw33( JI=R*Aij3 3(J)^£XP(-FU3( J)*E/T) 

547 


XR33U)*K*AR33( J)-^cXP(-FU3a)*r/T) 

550 


XP12I J)*R*AP12I JI*£XP(-FUlCKI*c/T) 

33 1 


XUl2( J)*R<^Awl2( JI^eXPl-FDlUl^e/T) 

552 


XR12( Jl=R*AR12(J)*cXP(-FUlCU*F/TI 

553 


XP13U1*R44P13( J)^«tXP<-FUl(K)*E/TI 

554 


X013< J)=R<‘AGl 3U)«tXP(-FUl (Jl*e/T) 

555 


XR13( JI*R4AR13CJH'tXP<-FDlCU«E/T) 

556 


XP23( J)-R«AP23U)4EXP|-FU2(K)*fc/T) 

557 


X023U)*R*AU23U)^£XP<-FU2( J>*E/T) 

560 


XR23CJ)=R<‘AR23( J)»cXPC-FU2IU*F/T| 

561 


XP3l(J)=R4AP3i ( J) ♦CXP<-FU3IK)*F/T) 

562 


XG3K J)*R*A03UJ)*£XP<-FU3CJ)»F/TI 

563 


XKjK JI=R*AR3lUH'cXP(-FU3(LI ♦6/TI 

564 


XP72( J)=R*AP32U)»eXP(-FU3(K)+r:/T> 

565 


Xu'32( J)=R*AU32< J) ‘•‘cXPJ -FUM J)*F/T) 

56to 


XR:>2( J)=R*AR32<J)<»tXP(-FU3<L)^E/T) 

567 


XP2U J)=R^AP21U)*tXP<-FU2CK)*F/T) 

370 


Xw2l( J)=R<'AG21U)*EXP(“FJ2lJ)*e/TI 

371 


XK2U Jl=R«AR21U)<'FXP(-FU2a>*F/T) 

3 72 

311 

CUNTINUi 

S74 

310 

CO iT I ^4Ut 


C Lisr “tSULTS 



I 


I 


III I III 


I Hill 


iiiiiniii in 


Ill II 


1>75 
57o 
y 77 
600 
oO 1 
o ‘.;2 
60 3 
fcJLO 
oil 
6lo 
6‘ 7 
620 
621 
t 

□ 33 
0*^0 
OM 5 
6*^6 
O'* 7 

Ot)<t 
oc. 1 
6f»6 
•:.73 
67^t 
o76 

702 
707 
7jOV 

721 

722 
7/3 
7^U 
7. >6 
7^2 
7-^7 

760 

761 
76b 

703 
770 

775 

776 

777 

i v .'0 4 

Isyll 
i. v/ i 6 

lv.^3 
t C 2 5 

iu?2 

ib37 

I 6 i 
lv 62 
’ :6 S 


31 FORMAT (3Xtl5F7.0) 

32 FORMATOX, 15F7.0) 

62 F0RMAT(10X,14F7*01 
53 F13-:MAT(17X,13F7.0» 

5^ FlJKMAT(23X,12F7.0) 

>5 FJHMAK ^X. 14r7.0) 

W.^:l TF(6,1^M J,J = lrl5) 

. 1<» F0KMAT(3X, 1517) 

WRITr(fe»l6M J,J = l6f29) 

15 F0KMAT(3X,14I7) 

WKlTe(6,l3l) 
lOJ FJHMATC^M J^ll) 

WRITt(o,5-') (Ril (K) ,K«3,15> 

WR I TE < 6, 5 3 M XP I U L ) , L*3 1 1 5) 
WRITE (6,55 J{ PI 1(K) tK^lb, 291 
WRITE (fc,56) (XPll (L) tL=l6f 29) 
WRlTE(o,lt)2) 

102 FURMATIVH Oil) 

WRITE (6,52 )( Oil (K) ,K=2,15) 
WRITE (6, 52 MXgll (L),L=2,15) 
WRITE (0,55 Moll (K) ,K=^16, 29) 
WRITE(6,56)(X0ll(L) ,L=16,29) 
WKlTE(o,103) 
i:)3 FOr<MAT(^H RID 

WRlTE(6r52) (KIDK) ,K«2,15) 
WRITE (6,52) (XRll (L) ,L=2,1 5) 
WRITE (6,55) (Kll (K) ,K=16,29) 
WRITc(o,55nXRll(U ,L = 1S,29) 
WRITE(6,10) 
iOh FGRMAT(^H l»?.2) 

WRITE(6,5j)(P22(K) ,K=3,15) 
WRITE(6,53 ) (XP22(L) ,L=3,15) 
WRI T£ (6,53) (P22(K) ,K^1 6, 2 9) 
WRITE (6,55 )(XP22(U ,L=1 6,29) 
WRITE (6, 105) 

105 FjRMAT( 4H 022) 
WRITE<6,52)(g22(K) ,K=2,15) 
WRITE(6,52 j(XU22(U ,L = 2,15) 
WRITE(6,55 )(u 22(K) ,K=16,29) 
WRIT6(6,55) (X022(L) ,L=16,29) 
WRITE (6, lOo) 

106 F0RMAT(4H R22) 

WRITE (6, 52 )(R22(K) ,K*2,15) 
WR1TE(6,52 ) (XR22(L) ,L=2,15) 
WRITE(6,56)(R22(K) ,K=l6f29) 
WRI TE(6, 55MXR22(U ,L=16,29) 
WRirc(6,l07) 

107 FURMAT(^H P33) 

WRITE (r,5^) (P33(K) ,K=^,15) 

WRI rE(o,5 O (XP33(L) ,L = ',,15) 

WR f re lo,5./MP33(K) ,K = l6f29) 
WRI TE (6,55 ) ( XP33(U ,L = 15,29) 
Wi< I TE ( 6, 1 >f5 ) 

ICa r J i,3 3) 

« D i iM >, 5 . ) { .,3MK ) ,K^3 , 1 5 J 
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1060 
lw65 
i r. 11 

L/ii 

i too 

I io(y 
IU3 
U20 

1125 

1126 
1127 
1130 
ii35 
11^2 
1167 
i i 56 
li55' 
1156 
lie>3 
1 170 
! 175 
1^:02 
1 3 
i 2' 
i^ll 
1216 
1223 
i 230 
4 2 31 
1232 

I 23 ^ 
i 266 
12 M 

1256 

1257 
1^60 
1265 
1272 
1277 
1306 
13C5 
1306 
1313 
1320 
1325 
1332 
13 33 
1 336 
1361 
1 366 
1353 
1 360 
i 3M 
\ 362 


WRITE (6» 53 HXQ33(U »L>3tl5l 
WRITE(6,55HQ33CK) »K«I6»29I 
WRITE(6f 55)(Xu33<U»Lxl6f29) 
WRITf=l6,109) 

109 FJRMAT(6H X33) 

^^lUTc (6,52 MR33(K| »K=2f 151 
WRITE(7>,52 MXR33(L>,L:=2tl5l 
WRI Tc(6, 55Mk33(K) rX*l6,29) 
WKI TE (6,55) (XK33CU ,L = 16,29| 
WRITE (6, 200) 

WRITE(6,11C) 

110 FORMAT(*»H P12> 

WRITE (6,5 i ) (P12(K) ,K = 3,15) 
WRIT£(6,53)(XP12(U ,L*3,15) 
WR1T6(6,55)(P12(K) ,K*16,29) 
WRITE (6,55) (XP12(U ,L*16, 29) 
WkirE(6,lll) 

111 F0RMAr(6H 012) 

WRITE (6, 52 ) (Q12(K) ,K=2,15) 
WRI re(6,52)(Xgi2(L) ,L=*2,15) 
WRITE ( 0,50 m 012(K) ,K = 16,29) 
WaiTE(6,5:>)(X(Jl2(L) ,L=16,29) 
WKITE(6,112) 

112 F0RMAT(6H R12) 

WiUTE(o,52) (K12(KI ,K«2,15) 

WR ITb (6,52) (XR12(L) ,L»2,15) 
WTITE(6,55) (R12(K) ,K=16,29) 
WRITE (6,55) (XR12(U ,L = 16,2 9) 
WRITE(6,H3) 

113 F3RHAT(6H M13) 

•WKlTH(t»,5 3)(P13(K) ,K = 3,15) 

WR irE(6,5i)(XP13(U ,L=^3,15) 
WRITc(o,55) (P13(K) ,K=16,29) 
WRI TE (6,55 )(XP13(L) ,L-16,29) 
WRITE (6,116) 

116 FOKMAT(^H oI3) 

WRITE(6,52)(U13(K) ,K=2,15) 
WRITE(6,52 i(X0l3(L) ,L=2,15) 
WRITE(6,55)(U13(K) ,K=16,29) 
WR I TE (6,55) (Xul3(L) ,L» 16,29) 
WRITE(o,U5) 

115 F0RMAT(6H R13) 
WRITE(o,52)(Rl3(K) ,K=2,15) 
WRITc(6,52 ) (XR13(L) ,L=2,15) 
WRirE(6,5D)(R13(K),K=l6,29) 
WRITE (6,55 )(XR13(L) ,L=1 6,29) 
wRiTE (6,116) 

116 FQRNiAT(6H P21) 

WRI TE (6,53) (P2KK) ,K = 3,15) 
WRI re(6,5A)(XP21(U ,L = 3,15) 
WRITE (6,55) (P2KK) ,K = 16,29) 
WRITE (6,5.>)(XP2l(L) ,L = 16,29) 
WRITl(6,1; 7) 

117 F^J•^M/4T^'♦^! 021) 

WR llh ( 6,52 ) ( 021 (K ) ,K = 2 , I''*) 


i it? 
1374 
1401 
1406 
^407 
1410 
1415 
I4k2 
14;.7 
.434 

1435 

1436 

1437 
1444 
1451 
145 6 

1463 
i 464 
i 466 
14/2 
1477 
1504 
i > i 1 
1p12 
I >13 
1520 
1515 
1632 
i 537 
i 54(J 
1541 
1 546 
1556 
1560 
1565 
I ?06 
I .>6 7 
: -»74 

16.’ A 
). oC*6 

1613 

1614 
lol5 
lo22 
lo27 
lo34 
lo4i 
1642 
i o4 3 
1650 
1655 
A t>62 
i 06 ? 
lr.70 
i .»7i 


MR I TE < 6, *^2 ) ( XQ21 ( U « L-2 » 1 5 > 
MRITE(6,5SMQ21(KI,K*16,29I 
WRIT£(6,55)JXu2UU»L*l6,29> 
MRITe«6,ll8) 

1X8 F:jRMAT<<*H R2l) 

MRITf.(6,b2)(R2l(K>fK*2f 15J 
WRIT£|6,52)»XR2HU.L»2,I5I 
wRirE(6,55)(«21(K)fK=»l6,29l 
WR 1 Tfc ( 6t 55 MXR2I (LI »L“l6t 291 
H-<irE<bt200) 

WRITL-(6,X19) 

1X9 FDHMAH'tH P32) 

WRITE(e>,5<.»(P32CKI tK=>^tX5) 
WRITE(6f bHUXP32(U tL*<>tX5» 
WRITE(6,5-3»(P32(M t<*16t29l 
wairS(t),5^){XP32(U tL=*X6,29l 
WR ITE (6, X20» 

120 FllRMAT(-»H Ci32) 

WK ITE ( Of 53 ) 1032(X) f K = 3.X5» 
WRITE (Of 1.3 MXU32IUfL = 3fX5> 
WRITE(6,f>o)»Q3?(X)fX»X6,29» 
WRITE (6,53 Mxa32tU tL = XSt 291 
W!UT6(0,X2X) 

X21 F0RHAT(AH R32) 

WK.IT£(c,f52MR32(K» ,K=2,i5) 
WRITE (6,52 )(XR32(L» fL=2,l 5 » 
WR! TE (6,55) (R'32(K» ,K = 16,29) 
WRITE. (6,55) (XR32(L),L=X6,29» 
WRlTf;(o,X22) 

X22 FORMAT(hH P23) 

WR1TE(6,53) (P23(X) ,K=3,X5) 

Wk I IE ( 6,53 ) t XP23(L) ,L»3, X5) 
WR1TE(6,55)(P23(K) ,K=16,29) 
WRITE (6,55) (XP23(L) ,L=X6, 29) 
WRITE) 6, 123) 

X23 F0RMAT(4H 023) 

WRIT£(6,52 )(U23(K) ,K=2,X5) 
WRITF(6,52i(XU23(L),l-«2,X5) 
WRITc(6,55)(a23(K),X=16,29) 
WRIT£(6,55)(XQ23(L) ,L=X5,29I 
WR(T£(6,X2A) 

X2A F3RMAT(-vH R23) 

WRITE(6,32)(k23(K),X»2,15) 
WRIT£(6,52)(XR23(L) ,L=2,X5) 
WRITE(6,55)(R23(K) ,K»X6,29) 
WR1TE(6,55)(XR23(L) ,L=X6,29) 
WRITE (6, X25) 

125 FORMAT (AH P31) 

HRITE(6,5A)(P3X(K),X=A,15) 
WRITE(6,5A) (XP31(L),L»A,X5) 
WRIT£(6,55)(P3l(K) ,K«X6,29) 
WRI TE(6,5':.)(XP3X(L) ,L = X6,29) 
rtRITE(6,l26) 
i2t> F !,'MAT(A|. .jU) 

«<l It (6,-.'l(:'TI(M,K = 3,l5) 



lo76 

1703 

1710 

i71t> 

1716 

1717 
17::4 
ir3i 
i 736 

1743 

1744 
I /46 


WRire(6,53MXg31(U »L»3fl5l 
WRITE (6« 53) (Q31 (K) »K*16»29} 
WR1TE(6»53MXQ31(U tL«16t29) 
WRITE(6,127) 

127 F0KMAT(4H R31) 

WRlTE<&«32)(K31(K)«K»2tl5) 
WRlTE(b,52)(XR31IL)#L*2tl5) 
WRlTk (Ut53 Mr<31(K) fK>16f29) 
WRITE (6»53HXR31(L) tL»16»29l 
GO TO 

99 CO'iTI.iUE 
E'^D 
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